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within  the  framework  of  applicability  of  its  basic  concepts. 
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Introducción 


Es  este  un  libro  sobre  Termodinámica  destinado  a  los  estudiantes  de  la  Licenciatura  en 
Física  y  otras  carreras  afines  en  esta  Universidad.  Durante  mi  desempeño  como  profesor  he 
tenido  oportunidad  de  dictar  numerosos  cursos  de  Termodinámica  y  varios  de  Mecánica 
Estadística  -base  microscópica  de  la  termodinámica-  para  las  licenciaturas  en  Física  y  en 
Astronomía.  A  partir  de  mis  notas  de  clase  he  ido  elaborando  los  capítulos  que  constituyen  el 
libro,  tratando  de  resumir  en  ellos  lo  esencial  de  cada  clase. 

En  el  breve  texto  presente  se  pretende  integrar  los  conocimientos  sobre  termodinámica 
previamente  adquiridos,  en  una  formulación  basada  en  postulados  que  se  constituye  en  una 
teoría  general  de  la  materia  y  la  energía.  En  el  desarrollo  lógico  posterior  a  partir  de  los 
postulados,  encuentran  su  lugar  las  diversas  variables  macroscópicas  conocidas  de  la  física 
general.  Analizados  los  criterios  de  estabilidad  se  introducen  las  transiciones  de  fase  asociadas 
a  la  violación  de  dichos  criterios  e  ilustradas  con  una  rica  fenomenología  en  diversos  campos. 
Luego  se  discuten  los  fenómenos  críticos  concluyendo  en  la  necesidad  de  un  postulado 
adicional:  la  hipótesis  de  escala.  En  los  cursos  habituales  se  aborda  finalmente  el  dominio  de  la 
termodinámica  fuera  de  equilibrio  y  las  transiciones  de  fase  disipativas  tendiendo  un  puente 
hacia  el  estudio  de  los  sistemas  biológicos  y  otros  sistemas  complejos.  Estos  temas  han 
quedado  fuera  del  presente  texto  y  tal  vez  sean  objeto  de  una  próxima  obra. 

El  objetivo  principal  de  este  libro  es  cubrir  adecuadamente  los  contenidos  de  una  materia 
fundamental  para  cualquier  licenciado  en  física  y  otras  disciplinas  afines.  Con  ello  pretende 
dotar  a  los  estudiantes  de  una  herramienta  universal  sencilla  para  encarar  el  estudio  de 
cualquier  sistema  macroscópico,  sus  propiedades  de  equilibrio  y  la  posible  existencia  de 
diversas  transiciones  de  fase.  Como  objetivos  específicos,  se  pueden  mencionar: 

-  ofrecer  una  versión  unificada  de  las  propiedades  macroscópicas  aprendidas  en  física  general; 

-  brindar  un  método  deductivo  sencillo  para  predecir  el  comportamiento  de  sistemas 
macroscópicos; 

-  capacitar  para  el  análisis  de  muy  diversas  transiciones  de  fase  posibles  y  reconocer  la 
universalidad  de  algunos  comportamientos. 


Luis  Mendoza  Zélis 
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CAPÍTULO  1 

Repaso  y  formalización  de  la  termodinámica 


En  este  capítulo  recordaremos  las  leyes  de  la  termodinámica  de  equilibrio,  aprendidas  en 
los  cursos  de  física  general.  A  partir  de  ellas  elaboraremos  una  formulación  en  base  a 
Postulados  que  facilitará  el  desarrollo  posterior. 


Las  leyes  de  la  Termodinámica 

La  termodinámica  es  la  ciencia  de  lo  macroscópico.  En  nuestro  estudio  físico  de  la 
naturaleza  nos  enfrentamos  a  sistemas  constituidos  por  un  número  enorme  de  partículas  (es 
en  este  sentido  que  los  llamamos  macroscópicos).  A  la  descripción  microscópica  provista  por 
las  leyes  fundamentales  de  la  física  debemos  unir  las  consecuencias  del  gran  número  de 
partículas  constituyentes.  De  aquí  surge  una  serle  de  propiedades  generales  que  todo  sistema, 
más  allá  de  sus  detalles  específicos,  debe  cumplir:  las  leyes  de  la  Termodinámica.  Se 
constituye  así  ésta  en  una  especie  de  teoría  general  de  la  materia.  Las  leyes  de  la 
termodinámica  tienen  una  fuerte  base  empírica  y  fueron  enunciadas  hace  ya  más  de  cien  años. 
A  continuación  recordamos  sus  aspectos  fundamentales. 


Ley  0 

Existe  una  propiedad  de  los  estados  de  equilibrio  de  los  cuerpos,  llamada  temperatura. 

Dos  cuerpos  puestos  en  contacto  mutuo  pueden  cambiar  sus  propiedades  hasta  alcanzar 
un  estado  de  equilibrio  térmico  entre  ambos,  a  partir  del  cual  éstas  no  cambian  más:  decimos 
entonces  que  los  cuerpos  se  encuentran  a  la  misma  temperatura.  Esta  relación  tiene  carácter 
transitivo:  dos  cuerpos  en  equilibrio  térmico  con  un  tercero  estarán  en  equilibrio  térmico  entre 
sí.  Esto  nos  permite  construir  termómetros  y  definir  acabadamente  esta  nueva  propiedad:  la 
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temperatura  o  estado  térmico  de  los  cuerpos.  A  partir  de  esta  ley  podemos  además  distinguir 
entre  paredes  diatérmicas  (las  que  permiten  que  se  alcance  el  equilibrio  térmico  mediante  un 
proceso  que  lleva  a  que  ambos  cuerpos  tengan  la  misma  temperatura)  y  paredes  adiabáticas 
(las  que  impiden  este  proceso). 

Antes  de  continuar,  es  conveniente  en  este  contexto  definir  con  cierto  rigor  matemático  la 
noción  de  propiedad  o  función  de  estado. 

F  es  función  de  estado  <=>  £dF  =  0  <=>  dF  es  un  diferencial  exacto. 

En  dos  variable  independientes  x  e  y,  la  forma  diferencial  X(x,y)dx  +  Y(x,y)dy  es  un  diferencial 
exacto  si  y  sólo  si  3X18 y  =  8Y I8x.  Para  la  propiedad  F  =  F(x,y), 

dF  =  (3F/3x)dx  +  (3F/5y)dy  y  entonces: 

c>X/dy  =  d2FI3ydx  =  32FI3xdy  =  dYIdx 

si  las  derivadas  segundas  existen  y  son  continuas. 

Además: 


fdF  =  f[(5F/5x)dx+(5F/5y)dy]  =  fdF[(5F/5x)i+(5F/5y)j].(dxi+dyj) 
e  introduciendo  Z  =  (5F/5x)i+(5F/5y)j  y  ds  =  (dx)i+(dy)j 
fdF  =  fz.ds  =  JJrotZ  dxdy  =  0 
pues  rotZ  =  dYIdx  -  dXIdy  =  52F/5x5y  -  02F/Sy5x  =  0. 


Ley  1 

La  energía  de  un  sistema  aislado  se  conserva. 

Durante  el  proceso  que  lleva  al  equilibrio  térmico  entre  dos  cuerpos,  éstos  intercambian 
energía  (llamamos  calor  Q  a  esta  forma  de  intercambio).  Además  la  energía  de  un  sistema 
puede  modificarse  realizando  trabajo  W  sobre  él,  que  es  una  forma  controlable  y  medible  de 
intercambiar  energía.  De  esta  manera  el  cambio  en  la  energía  interna  U  será  en  general:  AU  = 
Q  +  W  y  en  un  proceso  adiabático  (Q=0):  AU  =  Wad.  Esto  nos  permite  reconocer  a  la  energía 
interna  U  como  una  función  de  estado,  es  decir,  una  propiedad  característica  del  dado  estado 
de  equilibrio  del  sistema. 

Claramente,  en  un  sistema  aislado  (Q=W=0),  AU  =  0.  Microscópicamente,  esta  ley  es 
consecuencia  de  que  todas  las  interacciones  entre  partículas  constituyentes  responden  a 
fuerzas  conservativas. 
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Ley  2 


La  entropía  de  un  sistema  aislado  nunca  disminuye. 

No  todos  los  procesos,  aún  cuando  conserven  la  energía,  ocurren  espontáneamente.  En 
general,  si  naturalmente  ocurre  el  proceso  A^-B,  entonces  no  se  observará  el  proceso  inverso 
B^A.  Hay  una  dirección  natural  para  la  ocurrencia  espontánea  de  los  procesos.  Esta 
irreversibilidad  macroscópica  se  contrapone  a  la  reversibilidad  de  las  leyes  microscópicas  que 
describen  las  interacciones  entre  partículas  constituyentes. 

En  términos  del  rendimiento  q  de  una  máquina  térmica,  que  toma  Q-i  de  un  foco  a  T,  y 
entrega  -Q2  a  otro  a  T2,  el  teorema  de  Carnot  establece  que:  q  =  |W|/Qi  =  1-|Q2|/Q1  <  1  — T 2/T i . 
De  aquí  surge  que  debe  cumplirse:  0-,/T-i  <  |Q2|/T2  y  teniendo  en  cuenta  que:  |Q2|  =  -Q2, 
podemos  concluir  que:  Q^T-,  +  Q2/T2  <  0.  Resultado  que  Clauslus  generalizó  para  cualquier 

proceso  cíclico:  fdQ/T  <  0 

Se  define  entonces  una  nueva  función  de  estado  llamada  entropía  S,  a  través  de: 

AS  =  jrev(dQ/T)  donde  la  integral  debe  hacerse  a  lo  largo  de  un  proceso  cuaslestátlco 

(reversible).  Para  los  procesos  espontáneos  (irreversibles)  vales  la  desigualdad  AS  >  j(dQ/T). 

Claramente,  para  un  sistema  aislado  (dQ=0),  AS  >  0,  dependiendo  de  las  características  del 
proceso.  Ello  permite  caracterizar  a  los  estados  de  equilibrio  de  un  sistema  aislado  como  los  de 
máxima  entropía  S  para  el  dado  valor  de  la  energía  interna  U. 

Microscópicamente,  cada  estado  de  equilibrio  macroscópico  ( macroestado )  es  compatible 
con  un  gran  número  de  configuraciones  microscópicas  ( microestados ).  Boltzmann  estableció 
una  relación  entre  la  entropía  del  macroestado  y  el  número  W  de  microestados  consistentes 
con  el  mismo,  en  un  sistema  aislado: 

S  =  k.ln(W) 


Ley  3 

La  entropía  de  cualquier  sistema  se  anula  en  el  estado  T-0. 

Se  observa  que  los  cambios  en  la  entropía  en  procesos  Isotérmicos  se  van  reduciendo  al 
disminuir  la  temperatura,  de  modo  que  AST  -»  0  cuando  T-»  0.  Ello  llevó  a  formular  esta  ley, 
que  sólo  pudo  justificarse  con  el  posterior  desarrollo  de  la  mecánica  cuántica. 
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Formulación  en  base  a  postulados 


Estas  cuatro  leyes: 

Existe  una  propiedad  de  los  estados  de  equilibrio  de  los  cuerpos,  llamada  temperatura. 

La  energía  de  un  sistema  aislado  se  conserva. 

La  entropía  de  un  sistema  aislado  nunca  disminuye. 

La  entropía  de  cualquier  sistema  se  anula  en  el  estado  T=0. 

son  la  base  de  la  Termodinámica  y  de  ellas  puede  deducirse  casi  todo  cuanto  está  contenido 
en  este  libro.  Esta  tarea  se  simplifica  considerablemente  reformulando  la  teoría  en  base  a  una 
serie  de  postulados  (afirmaciones  con  verificación  a  posterlori).  Basados  en  el  texto  de  H. 
Callen  formularemos  entonces  tres  postulados  sobre  los  que  puede  fundarse  la 
Termodinámica. 

El  primero  se  refiere  a  la  posibilidad  de  describir  los  estados  de  equilibrio  de  los  sistemas 
simples  en  términos  de  un  pequeño  número  de  variables.  Denominamos  sistema  simple  a  todo 
sistema  macroscópicamente  homogéneo,  Isótropo,  sin  carga,  químicamente  inerte, 
suficientemente  grande  para  poder  despreciar  efectos  de  superficie,  no  afectado  por  campos 
externos  (eléctricos,  magnéticos  o  gravltatorios). 


Postulado  I.  En  los  sistemas  simples  existen  estados  de  equilibrio  que  macroscópicamente 
quedan  completamente  caracterizados  por  la  energía  interna  U,  el  volumen  V  y 
los  números  molares  Ni,  N2,...,Nr. 


Como  se  ve,  el  mismo  postulado  admite  la  existencia  de  estados  de  equilibrio,  en  el  sentido 
habitual  de  permanecer  Invariables  en  el  tiempo,  y  establece  las  propiedades  que  los 
caracterizan:  U,V,N1,N2,...,Nr.  Notemos  que  todas  ellas  son  propiedades  extensivas 
(proporcionales  al  tamaño  del  sistema).  Para  describir  sistemas  de  mayor  complejidad 
deberemos  incluir  nuevas  variables  (momento  magnético,  eléctrico  u  otras  propiedades). 

Se  sobreentiende  que  el  sistema  así  definido  tiene  límites  (paredes)  que  lo  aíslan  del  resto 
del  universo.  Esto  significa  que  estas  paredes  deben  ser  restrictivas  a  la  modificación  de 
cualquiera  de  los  parámetros  que  fijan  el  estado  de  equilibrio.  Deben  ser  rígidas  para  Impedir  la 
modificación  de  V  y  de  U  (a  través  del  trabajo  asociado  al  cambio  de  volumen),  deben  ser 
Impermeables  para  Impedir  la  modificación  de  cualquiera  de  los  N¡  y  de  U  (a  través  del  trabajo 
asociado  al  intercambio  de  partículas)  y  deben  ser  adiabáticas  para  Impedir  la  modificación  de 
U  a  través  del  intercambio  de  calor. 


9 


universo 


pared  rígida, 
impermeable 
y  adiabática 


Figura  1.1 


En  general  nos  interesa  estudiar  sistemas  abiertos  total  o  parcialmente  al  resto  del  universo 
(que  se  considera  en  sí  mismo  un  sistema  aislado)  mediante  paredes  dlatérmicas  y/o  móviles 
y/o  permeables.  Esto  es  equivalente  a  estudiar  sistemas  compuestos  de  dos  sistemas  simples, 
separados  entre  sí  por  paredes  de  un  dado  tipo  y  aislados  del  resto  del  universo.  Las  paredes 
internas  permitirán  o  no  el  intercambio  de  una  o  varias  de  las  propiedades  extensivas  (U,V,N¡) 
entre  los  sistemas  constituyentes,  manteniendo  fijo  el  total  de  cada  una  de  ellas  mediante  la 
pared  externa.  En  el  caso  trivial  en  que  ambos  subsistemas  estén  aislados  entre  sí  (pared 
completamente  restrictiva)  el  estado  de  equilibrio  del  sistema  quedará  determinado  por  el 
equilibrio  Individual  de  los  constituyentes.  El  problema  que  se  nos  presenta  es  el  de  predecir 
cuál  será  el  estado  de  equilibrio  hacia  el  que  evolucionará  el  sistema  cuando  se  remueva 
alguno  de  los  vínculos  o  restricciones  internas,  es  decir  se  modifique  total  o  parcialmente  el 
carácter  de  la  pared.  Para  ello  formulamos  el  siguiente  postulado  (dos  postulados  en  la  versión 
original  de  H.  Callen): 


Postulado  II:  Existe  una  función  continua  y  diferenciable  de  los  parámetros  extensivos  de 
cualquier  sistema,  monótonamente  creciente  con  U,  llamada  entropía  S  y 
definida  para  todos  los  estados  de  equilibrio.  Para  un  sistema  compuesto  la 
entropía  es  aditiva  sobre  los  subsistemas  constituyentes  y  al  remover  un  vínculo 
interno,  el  estado  de  equilibrio  que  el  sistema  alcance,  será  aquél  que  la 
maximice  sobre  todos  los  estados  compatibles  con  dicho  vínculo. 


Vemos  que  el  postulado  introduce  una  nueva  función  de  estado,  dependiente  de  las  otras  a 
través  de  la  llamada  relación  fundamental :  S  =  S^.V.NL,...),  a  la  que  además  de  ciertas 
propiedades  le  agrega  un  principio  extrema!  claramente  Inspirado  en  la  Ley  2.  Al  remover  un 
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vínculo  interno,  el  sistema  compuesto  elegirá  un  estado  entre  muchos,  cada  uno  de  ellos 
siendo  un  estado  de  equilibrio  realizable  con  el  vínculo  apropiado. 

Veamos  esto  a  través  de  un  ejemplo.  Supongamos  que  en  el  sistema  de  la  figura  1.2,  la 
pared  interna  es  dlatérmlca,  rígida  e  Impermeable.  Entonces  U1  y  U2  podrán  variar, 
redistribuyendo  la  energía  entre  los  subsistemas,  sujetos  a  la  condición  de  cierre  U1  +  U2  =  U  = 
cte.  Todos  los  valores  de  U1  (y  los  correspondientes  de  U2)  son  posibles  de  realizar  usando  una 
pared  adiabática,  pero  el  valor  que  adopte  el  sistema,  al  permitírsele  variar  mediante  una  pared 
dlatérmlca,  será  el  que  maximlce  la  entropía  total:  S  =  S1(U1,V1,N11,...)  +  S2(U2,V2,N-|2,...). 

pared  rígida, 
impermeable 


y  díatérmíca  universo 


Entonces  para  variaciones  en  torno  al  equilibrio: 

5S  =  5S1  +  5S2  =  (5S1/5U1)v(i),Ni(i)-8U1  +  (5S2/5U2)v(2),n¡(2)-5U2  =  0 
5S  =  [(SSVsuVuNid)  -  (5S2/5U2)v(2),n¡(2)].SU1  =  0;  V5U1. 


y  la  condición  necesaria  de  equilibrio  será: 

(0S1/ SU 1  )v(1  ),Ni(1 )  =  (3S2/5U2)v(2),N¡(2)- 


Además,  si  el  sistema  se  encontrara  ligeramente  fuera  de  equilibrio,  por  ejemplo  con: 

(5S1/ SU 1  )V(1  ),Ni(1 )  >  (ciS2/0U2)v(2),N¡(2) 


el  proceso  que  restablezca  el  equilibrio  deberá  aumentar  S.  Entonces: 

5S  =  [(¿5S1/¿?U1  )v(i),N¡(i)  -  (ciS2/3U2)v(2),n¡(2)]-SU1  >  0 
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y  como  en  la  situación  supuesta  [  ]  >  0,  será  entonces  5U1  >  0  (5U2  <  0).  Es  decir  que  fluirá 
energía  en  forma  de  calor  desde  el  sistema  2  al  1 ,  desde  el  de  menor  al  de  mayor  (oS/oU)v,ní- 

Esta  cantidad  tiene  entonces  las  propiedades  de  1/T. 

De  manera  análoga  se  puede  analizar  la  redistribución  de  volumen  o  de  partículas  que 
conduce  al  equilibrio  en  el  caso  de  paredes  diatérmicas  móviles  o  permeables.  Así 
obtendríamos  respectivamente  las  condiciones  de  equilibrio  adicionales: 

(¿?s1/¿5V1  )u(i),N¡(i) =  (as2/av2)u(2),N¡(2)  y 

(c,S1/3N¡1)U(i)iv(1),Nj?í¡(1)  -  (c,S2/3N¡2)U(2)iv(2),Nj?í¡(2)- 

El  criterio  de  maximizar  S  para  los  estados  de  equilibrio  puede  extenderse  a  los  sistemas 
simples.  Estos  pueden  considerarse  idealmente  divididos  en  dos  o  más  partes  por  paredes 
imaginarias  no  restrictivas.  La  distribución  de  U,V,Ni,...  entre  los  subsistemas  es  tal  que 
maximiza  S  para  los  dados  valores  totales  de  U.V.Ni,...  . 

Notemos  además  que  de  la  aditividad  sobre  subsistemas  constituyentes: 

S  =  Ea  S“  con  S“  =  Sa(Ua,Va,N1°,...) 

surge  la  extensividad  de  la  entropía.  En  efecto,  si  consideramos  un  sistema  compuesto  de  X 
subsistemas  iguales: 


S(A.U,XV,XN¡)  =  Za  S“(U,V,N¡)  =  A.S(U,V,N¡). 


Ello  es  equivalente  a  decir  que  la  entropía  es  una  función  homogénea  de  grado  1  de  todas  sus 
variables  y  nos  permite  escalear  las  propiedades  de  un  sistema  de  N  moles  a  las  de  1  mol. 
Eligiendo  A,=1/N,  con  N=2N¡,  e  introduciendo  las  magnitudes  molares  s=S/N,  u=U/N,  v=V/N  y 
las  fracciones  molares  x¡=N¡/N: 

S(U,V,N! . Nr)  =  . . . 

Y  para  un  sistema  de  1  componente: 

S(U,V,N)  =  Ns(u,v)  y 
s  =  s(u,v). 

Finalmente  veamos  que,  por  las  propiedades  de  la  función  S,  la  relación  fundamental  S= 
S(U,V,N-i,...)  puede  invertirse,  obteniéndose  U  =  U(S,V,N-|,...)  que  también  será  función 
homogénea  de  grado  1  de  todas  sus  variables.  Pasamos  así  de  la  llamada  representación 
entrópica  a  otra  completamente  equivalente:  la  representación  energética.  Teniendo  en  cuenta 
que  (3S/3U)v,n¡  >  0,  el  principio  extremal  puede  reformularse  para  caracterizar  el  estado  de 
equilibrio  de  un  sistema  compuesto  aislado  como  el  de  mínima  energía  interna  U  para  el  dado 
valor  de  la  entropía  S. 
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Figura  1.3 

La  equivalencia  entre  ambas  formulaciones  puede  demostrarse  por  reducción  al  absurdo, 
como  muestra  el  esquema  de  la  figura  1.3.  Partiendo  de  un  estado  de  equilibrio  con  la  máxima 
S  para  el  dado  valor  de  U,  afirmamos  que  esa  es  la  mínima  U  para  el  dado  valor  de  S.  Si  ello 
no  fuera  así,  sería  posible  quitarle  energía  interna  al  sistema  en  forma  de  trabajo,  manteniendo 
constante  la  entropía  (U,S)  — ►  (U-AU,S).  Si  luego  le  devolviéramos  la  misma  cantidad  de 
energía  en  forma  de  calor,  restableciendo  el  valor  de  U  pero  aumentando  la  entropía  (U-AU,S) 
— ►  (U,S+AU/T),  no  se  cumpliría  el  principio  extremal  original  (máximo  S  para  el  dado  U).  Puede 
hacerse  la  analogía  con  la  definición  de  un  círculo  como  la  figura  de  mayor  área  para  el  dado 
perímetro  o  la  de  menor  perímetro  para  el  área  dada. 

Para  completar  este  esquema  formal  agregamos  un  tercer  postulado  que  coincide  con  la 
Ley  3  de  la  termodinámica. 


Postulado  III  La  entropía  de  cualquier  sistema  se  anula  en  el  estado  para  el  que 
(cU/3S)vm  =  0. 


Esquema  formal  básico 


Es  claro  que  lo  esencial  de  la  descripción  termodinámica  de  los  estados  de  equilibrio  de  un 
sistema,  contenido  tanto  en  la  Leyes  como  en  los  Postulados,  es  la  existencia  de  una  relación 
fundamental  que  describe  los  estados  de  equilibrio  y  un  principio  extremal  que  permite  predecir 
su  evolución.  En  este  sentido,  podemos  citar  a  Roger  Balian: 

Los  estados  de  equilibrio  de  un  sistema  están  caracterizados  a  escala  macroscópica  por 
un  conjunto  de  variables  extensivas  A¡  y  por  una  función  de  estas  variables,  la  entropía  S, 
que  es  continuamente  diferenciable,  positiva  y  aditiva-,  la  entropía  de  un  sistema 
compuesto  es  la  suma  de  las  de  sus  partes,  la  entropía  de  una  sustancia  homogénea  es 
extensiva.  En  un  sistema  compuesto  aislado  la  liberación  de  alguna  restricción  (interna) 
puede  permitir  intercambios  entre  los  subsistemas,  que  se  reflejan  en  cambios  en  los  A¡; 
el  dominio  A  permitido  para  las  variables  A¡  está  limitado  por  las  restantes  restricciones  y 
por  las  leyes  de  conservación.  En  el  estado  de  equilibrio  final  que  el  sistema  alcance  el 
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valor  de  las  variables  A¡  se  determina  buscando  el  máximo  de  la  entropía  en  el  dominio  A. 
(R.Balian,  1991, 244) 

Con  independencia  de  la  formulación  elegida  los  principios  básicos  pueden  resumirse  como 
sigue: 

•  En  los  sistemas  macroscópicos  existen  estados  de  equilibrio  homogéneos  que  pueden 
caracterizarse  completamente  por  un  conjunto  pequeño  de  variables  extensivas:  la  energía 
interna  U,  el  volumen  V  y  los  números  molares  N¡,  en  el  caso  de  sistemas  simples  y  algunas 
otras,  tales  como  momentos  eléctricos  o  magnéticos,  en  sistemas  más  complejos. 

•  Para  todos  los  estados  de  equilibrio  existe  una  función  de  estas  variables,  la  entropía  S, 
definida  a  través  de  la  ecuación  fundamental  del  sistema  (representación  entrópica): 

S  =  S(U,V,N¡)  que  contiene  toda  la  información  sobre  las  propiedades  macroscópicas  del 
sistema. 

•  En  un  sistema  compuesto  de  dos  o  más  subsistemas  simples  separados  entre  sí  por  paredes 
parcialmente  restrictivas  a  la  modificación  de  las  variable  mencionadas  (vínculos  internos)  y 
aislado  como  un  todo,  la  entropía  (aditiva  sobre  los  subsistemas  constituyentes)  toma  el  valor 
máximo  compatible  con  dichos  vínculos. 

•  La  relación  fundamental  puede  invertirse  (representación  energética)  U  =  U(S,V,N¡).  Los 
estados  de  equilibrio  de  un  sistema  aislado  son  los  que  minimizan  la  energía  para  el  dado  valor 
de  la  entropía. 


Referencias 

Balian,  R.  (1991).  From  microphysics  to  macrophysics:  Methods  and  Applications  of  Statistical 
Physics,  voll.  Springer. 

Callen,  H.  (1985).  Thermodynamics  and  an  i ntroduction  to  thermosatistics.  Wiley. 
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CAPITULO  2 

Propiedades  intensivas  y  relaciones  formales 


En  este  capítulo  veremos  cómo  aparecen  en  nuestro  formalismo  las  propiedades  intensivas, 
asociadas  a  las  derivadas  primeras  de  la  relación  fundamental.  Asimismo  desarrollaremos 
algunas  relaciones  formales  que  se  desprenden  fácilmente  de  los  postulados. 


Propiedades  intensivas  y  derivadas  primeras 


Para  la  descripción  de  los  estados  de  equilibrio  de  un  sistema  hemos  utilizado  sólo 
propiedades  extensivas.  Al  interesarnos  en  procesos  y  en  los  cambios  en  los  parámetros 
extensivos  asociados  con  ellos  necesitaremos  la  forma  diferencial  de  la  relación  fundamental 
que,  en  la  representación  energética,  puede  escribirse: 

dU  =  (aU/aS)v,N¡.dS  +  (aU/5V)s,N¡.dV  +  Z¡  (0U/0N¡)Siv,  Nj*i-dN¡ 


Daremos  a  estas  derivadas  parciales,  nombres  que  pueden  justificarse  fácilmente. 


(5U/aS)v,N¡  s  T 
(SU/dV)s,N¡  -  -P 
(0U/0N¡)s,v,Nj*¡  =  P¡ 


temperatura 

-presión 

potenciales  químicos 


Razonando  de  manera  análoga  a  la  de  la  clase  anterior,  podemos  analizar  el  equilibrio  entre 
dos  sistemas  separados  por  una  pared  dlatérmlca  (rígida  e  Impermeable)  en  la  representación 
energética,  y  asimilar  la  cantidad  (oU/oS)v,ní  con  la  temperatura  T  en  su  sentido  físico  habitual, 
justificando  la  anterior  denominación.  De  modo  similar  pueden  analizarse  las  condiciones  para 
el  equilibrio  entre  sistemas  separados  por  una  pared  móvil  y  dlatérmlca  en  cualquiera  de  las 

dos  representaciones.  En  la  representación  energética: 
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8U  -  [(¿?U  1/¿5S1  )v(i  ),n¡(i )-(¿?U2/ ¿5S2)v(2),n¡(2)] -5S1  + 

[(SU  VSV1  )S(i  ),ní(i  r(5U2/ 8V2)s(2),ní(2)]  -5V1  =0;  VSS'.VSV1. 


De  aquí  surge,  además  de  la  Igualdad  de  temperaturas,  la  Igualdad  de  la  cantidad  (5U/3V)s,n¡  a 
ambos  lados  de  la  pared  móvil.  Una  breve  reflexión  posterior  permite  asimilar  esta  derivada 
cambiada  de  signo  a  la  presión  P  con  su  sentido  físico  habitual.  Finalmente,  analizando  el 
equilibrio  entre  dos  sistemas  separados  por  una  pared  permeable  al  componente  1-éslmo  y 
dlatérmlca,  se  obtiene  como  condición  necesaria  la  Igualdad  de  la  cantidad  (3U/3N¡)s,v,nh  a 
ambos  lados  de  la  pared.  Esta  magnitud,  asociada  a  los  cambios  energéticos  producidos  por  el 
agregado  de  partículas  se  denomina  potencial  químico  del  dado  componente. 

Resulta  entonces  la  llamada  ecuación  de  Gibbs  (forma  diferencial  de  la  relación 
fundamental): 


dU  =  TdS  -  PdV  +  S¡  n¡dN¡. 


En  ella  reconocemos  el  término  de  trabajo  mecánico:  dWmec  =  -PdV  correspondiente  a  un 
cambio  cuaslestátlco  en  el  volumen.  Usamos  el  símbolo  dWmec,  en  lugar  de  dWmec,  para  Indicar 
que  no  se  trata  de  un  diferencial  exacto.  En  un  proceso  infinitesimal  a  N¡  =  cte, 
dU  =  TdS  -  PdV  =  TdS  +  dWmec,  con  lo  que  podemos  Identificar  el  término  de  calor:  dQ  =  TdS. 
Finalmente  en  un  proceso  infinitesimal  cualquiera,  cada  uno  de  los  términos  en  dN¡  puede 
asociarse  al  trabajo  químico  respectivo,  cambio  en  la  energía  producido  por  un  agregado 
cuaslestátlco  de  partículas:  dWq  ¡  =  p¡dN¡.  Entonces: 


dU  =  dQ  +  dWmec  +  S¡  dWqj  =  dQ  +  dW. 


Vemos  entonces  que  a  cada  variable  extensiva  se  le  asocia  una  variable  intensiva  como 
derivada  parcial  de  la  energía  respecto  de  aquella.  En  general,  Y¡  =  (oU/6X¡)Xj*¡  se  denomina 
variable  intensiva  conjugada  de  X¡.  En  un  sistema  simple  con  r  componentes  tendremos  r+2 
pares  de  variables  extensivas  e  Intensivas  conjugadas. 

Comprobamos  también  que  estas  propiedades  se  vinculan  directamente  con  las 
condiciones  de  equilibrio:  concretamente,  si  entre  dos  sistemas  se  permite  el  intercambio  de  X¡, 
el  equilibrio  entre  ambos  requiere  la  Igualdad  de  sus  respectivas  Y¡  (condición  necesaria). 
Trabajando  en  la  representación  entrópica  obtenemos  la  correspondiente  ecuación  de 
Gibbs: 
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dS  =  (1/T).dU  +  (P/T).dV  -  Z¡  (ja¡/T).dNi 


que  introduce  los  respectivos  parámetros  intensivos  entrópicos: 

(5S/aU)v,N¡  s  1/T 
(5S/5V)u,Ni  =  P/T 
(0S/5N¡)UiV,  Nj#¡  =  -P¡/T 


Ecuaciones  de  estado  y  otras  relaciones  formales 

Como  derivadas  primeras  de  la  relación  fundamental,  los  parámetros  intensivos  serán  a  su 
vez  funciones  de  las  variables  independientes  elegidas,  resultando  así  las  relaciones: 

T  =  T(S,V,N¡) 

P  =  P(S,V,N¡) 

M¡  =  (S,V,N¡) 

habitualmente  denominadas  ecuaciones  de  estado. 

Sugerencia:  buscar  entre  las  ecuaciones  conocidas  de  la  física  general,  de  origen 
fenomenológico,  agüellas  que  puedan  interpretarse  como  ecuaciones  de  estado. 

Cada  una  de  ellas  contiene  información  parcial  sobre  el  sistema  y  sólo  el  conjunto  de  las  r+2 
ecuaciones  posee  una  información  equivalente  a  la  de  la  ecuación  fundamental.  Estas 
relaciones  funcionales,  siendo  derivadas  primeras  de  la  relación  fundamental,  que  es 
homogénea  de  grado  1,  resultan  homogéneas  de  grado  0,  lo  cual  está  de  acuerdo  con  el 
carácter  intensivo  de  estas  propiedades.  Por  ejemplo: 


T(A.S,AV,A,N¡)  =  A,°.T(S,V,N¡ )  =  T(S,V,N¡ )  =  T(s,v,x¡). 


que  expresa  el  hecho  de  que  si  aumentamos  el  tamaño  del  sistema  simple  no  modificamos  su 
temperatura. 

Por  otra  parte,  a  partir  de  las  propiedades  de  homogeneidad  de  la  relación  fundamental: 


U(A.S,A.V,A.N¡)  =  XU(S,V,N¡), 


derivando  respecto  de  X: 

[aU(XS,XV,XN¡)/5(XS)].S  +  [5U(XS,XV,XN¡)/5(1V)].V  +  ...  =  U(S,V,N¡), 
T(1S,^V,^N¡)S  -  P(^S,XV,XN¡)V  +  E¡  =  U(S,V,N¡). 


T(S,V,N¡)S  -  P(S,V,N¡)V  +  2¡  ju¡(S, V, N¡)N¡  =  U(S,V,N¡). 
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Resulta  entonces  la  llamada  ecuación  de  Euler. 


U  =  TS  -  PV  +  Z¡  |i¡N|, 

donde  debemos  tener  presente  que  las  variables  independientes  son  S,V,N¡  y  deben  tomarse 
T,P,|u¡  como  funciones  de  aquellas. 

Finalmente,  diferenciando  la  ecuación  de  Euler: 
dU  =  TdS  +  SdT  -  PdV  -  VdP  +  Z¡  p¡dN¡  +  Z¡  N¡  dp¡ 


y  teniendo  en  cuenta  la  ecuación  de  Glbbs,  resulta: 


SdT  -  VdP  +  Z¡  NidjUj  =  0 


llamada  ecuación  de  Gibbs-Duhem.  Esta  pone  de  manifiesto  que  las  r+2  variables  Intensivas 
no  son  todas  Independientes,  sino  sólo  r+1  de  ellas.  Decimos  entonces  que  el  sistema  tiene  r+1 
grados  de  libertad. 

Para  sistemas  de  un  componente  con  U  =  U(S,V,N)  =  Nu(s,v),  la  ecuación  de  Glbbs  molar 
resulta: 

du  =  T.ds  -  P.dv 

y  es  habitual  escribir  la  ecuación  de  Gibbs-Duhem  como: 
dpi  =  -s.dT  +  v.dP 


En  la  Tabla  2.1  hacemos  un  resumen  de  estas  relaciones  formales  incluyendo  también  las 
correspondientes  a  la  representación  entrópica  y  los  casos  particulares  de  sistemas  de  1 
componente. 
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Tabla  2.1:  Formalismo  termodinámico  para  un  sistema  simple 


representación  energética 

representación  entrópica 

ecuación  fundamental 

U  =  U(S,V,Ni . Nr) 

S  =  S(U,  V,Ni,...,Nr) 

estados  de  equilibrio 
de  un  sistema  aislado 

mínimo  U 

para  los  dados  S,V,N¡ 

máximo  S 

para  los  dados  U,V,N¡ 

ecuación  de  Gibbs 

dU  =  T.dS-P.dV+Ei|i¡.dNi 

dS  =  (1/T)dU+(P/T)dV-S¡(M¡/T)dN¡ 

ecuaciones  de  estado 

T  =  T(S,V,Ni . Nr) 

P  =  P(S,V,Ni . Nr) 

Mi  =  p¡(S,V,Ni . Nr) 

1/T  =  1/T(U,V,Ni,...,Nr) 

P/T  =  P/T(U,V,Ni,...,Nr) 

Mí/T  =  m¡/T(U,V,Ni . Nr) 

ecuación  de  Euler 

U  =  TS-PV+I¡m¡N¡ 

S  =  U/T +PV/T-Ep¡N¡/T 

ecuación  de  Gibbs-Duhem 

S.dT-V.dP+XiNi.dpi  =  0 

U.d(1/T)+V.d(P/T)-ZiN¡.d(M¡/T)  =  0 

sistema  simple  de  un  componente 

ecuación  fundamental 

U  =  N.u(s,v) 

S  =  N.s(u,v) 

ecuación  de  Gibbs  molar 

du  =  Tds-Pdv 

ds  =  (1/T)du+(P/T)dv 

ec.  Gibbs-Duhem 

dp  =  -sdT+vdP 

d(p/T)  =  ud(1/T)+vd(P/T) 

En  un  sistema  genérico  tendremos: 

estados  de  equilibrio  descriptos  por  variables  extensivas  X¡, 

y  la  relación  fundamental  será  U  =  U({X¡}), 

habrá  variable  Intensivas  conjugadas:  Y¡  =  (0U/0X¡)Xj*¡ 

que  satisfacen  ecuaciones  de  estado:  Yj  =  Yj  ({X¡}) 

y  valdrán  la  ecuación  de  Gibbs:  dll  =  E¡  Y¡dX¡ , 

la  de  Euler:  U  =  S¡  Y¡X¡ , 

y  la  de  Gibbs-Duhem:  2¡X¡dY¡  =  0  . 

Ejemplo  1.1:  ecuación  fundamental  de  un  gas  ideal  monoatómico 

Como  aplicación  de  estos  conceptos  veamos  cómo  podemos  obtener  la  ecuación 
fundamental  de  un  sistema  simple  (gas  Ideal)  a  partir  del  conocimiento  parcial  que  dan  las 
ecuaciones  de  estado.  Por  un  lado  tenemos  la  conocida  ecuación  de  los  gases  Ideales:  PV  = 
NRT  y  por  otro  la  relación  U  =  NCVT  (con  Cv=  3R/2  para  un  gas  monoatómico).  Estas 
ecuaciones  pueden  reescrlblrse  para  poner  de  manifiesto  su  forma  de  ecuaciones  de  estado  en 
la  representación  entrópica,  ya  que  las  variables  extensivas  que  aparecen  son  U,V  y  N.  En 
efecto: 


1/T  =  3NR/2U  =  3R/2u 
P/T  =  NR/V  =  R/v 
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Un  procedimiento  posible  es  obtener  la  ecuación  de  estado  faltante,  la  que  daría  p/T  en 
términos  de  u  y  v,  integrando  la  ecuación  de  Gibbs-Duhem: 

d(p/T)  =  u.d(1/T)+v.d(P/T)  =  (3RT/2)d(1/T)+(RT/P)d(P/T) 

y  luego  reconstruir  la  ecuación  fundamental  usando  la  ecuación  de  Euler. 

Otra  alternativa  es  integrar  la  ecuación  de  Gibbs  molar:  ds  =  (1/T)du+(P/T)dv  y  luego 
obtener  S=Ns(u,v).  En  ambos  casos  la  constante  de  integración  quedará  a  determinar. 
Seguiremos  este  último  camino  que  aparece  como  más  simple: 
ds  =  3Rdu/2u  +  Rdv/v 

=>  s  =  s0  +  Rln(u/u0)3/2  +  Rln(v/v0)  -s0  +  Rln(u3/2v/u03/2v0), 
finalmente: 

S  =  Ns  =  (N/N0)S0  +  NR.In  [(U/U0)3/2  (V/V0)  (N/N0)_5/2]  . 

S  =  N.[cte.  +  R.ln(U3/2V/N5/2)] 

que  es  la  ecuación  fundamental  de  un  gas  ideal  monoatómico. 


Ejemplo2.1:  ecuación  fundamental  de  la  radiación  del  cuerpo  negro 

La  radiación  electromagnética  en  equilibrio  a  la  temperatura  T  en  el  interior  de  un  cuerpo 
negro  satisface  las  ecuaciones:  U/V  =  rT4  (con  r  =  4o/c,  siendo  a  la  constante  de  Stefan- 
Boltzmann)  y  P  =  U/3V,  donde  debe  notarse  la  independencia  con  N  (p=0).  Podemos  obtener 
de  ellas  las  ecuaciones  de  estado  en  la  representación  entrópica  (variables  U  y  V): 

1/T  =  1/T  (U,V)  =  r1/4(V/U)1/4 
P/T  =  P/T  (U,V)  =  (r1/4/3)(U/V)3/4 

De  allí,  por  sustitución  en  la  correspondiente  ecuación  de  Euler,  resulta  la  ecuación 
fundamental: 

S  =  (1/T)U  +  (P/T)V  =  (4/3).r1/4V1/4U3/4 

Notar  que  puede  escribirse  la  expresión  mas  sencilla:  S  =  (4/3). rVT3  que  no  es  una  ecuación 
fundamental.1 


Más  adelante  veremos  que  si  quieren  tomarse  T  y  V  como  variables,  la  relación  fundamental  debe  expresarse  en 
términos  de  la  función  de  Massieu,  transformada  de  Legendre  de  la  entropía,  S[1/T]  =  S  -  1/T.U  =  1/3  rVT3.  A  partir  de 
ella  podría  obtenerse,  por  ejemplo:  U  =  -  (í)S[1/T]/5(1/T))v  =  rVT4.  Por  otra  parte  es  interesante  notar  que:  G  =  U-TS+PV 
=  rVT4  -  (4/3)rVT4  -  (1/3)rVT4  =  0  como  era  de  esperar  a  partir  de  que  G  =  Np. 
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CAPÍTULO  3 

Procesos,  energía  utilizable  o  trabajo  mínimo 


En  este  capítulo  discutiremos  cómo  aplicar  el  principio  extremal  a  procesos  reales  en  los 
cuales  el  sistema  se  aparta  de  la  situación  de  equilibrio. 


Procesos 

Hasta  aquí  hemos  definido  estados  de  equilibrio  y  sólo  concebimos  los  procesos  como 
procesos  cuas¡estát¡cos\  sucesión  Ideal  de  estados  de  equilibrio.  Los  procesos  reales  son 
esencialmente  fuera  de  equilibrio  y  su  estudio  detallado  será  tema  de  un  texto  complementario. 
Como  todo  proceso  comienza  y  termina  en  un  estado  de  equilibrio,  algunas  conclusiones 
pueden  obtenerse  a  partir  del  cambio  en  las  funciones  de  estado,  que  tiene  lugar  como 
consecuencia  del  mismo. 

La  ecuación  fundamental  de  un  sistema  simple  S=S(U,V,Ni...Nr)  define  una  superficie  en  el 
espacio  de  configuración  termodlnámlco  (S,U,V,Ni,...,Nr).  Cada  punto  de  la  superficie 
corresponde  a  un  estado  de  equilibrio  y  un  proceso  cuaslestátlco  se  representará  como  una 
curva  o  trayectoria  sobre  la  superficie.  La  descripción  de  estados  y  procesos  en  sistemas  más 
complejos  o  en  sistemas  compuestos,  se  hace  de  manera  análoga  incorporando  las  variables 
adicionales  necesarias. 

Representar  un  estado  fuera  de  equilibrio  requerirá  de  un  espacio  de  mayor  dimensión  (más 
variables)  y  es  Imposible  mostrarlo  en  ese  gráfico.  En  un  proceso  real  el  sistema  desaparece 
de  la  superficie  para  volver  a  aparecer  cuando  recupera  el  equilibrio.  Ello  ocurrirá  más  o  menos 
rápidamente  dependiendo  de  la  velocidad  con  que  se  conduzca  el  proceso  (tiempos 
característicos  t)  y  de  los  tiempos  de  relajación  propios  del  sistema  (x).  Los  diferenciales  dS, 
dU,  ...  que  hemos  considerado  y  que  satisfacen  la  ecuación  de  Glbbs  se  refieren  a  procesos 
infinitesimales  cuaslestátlcos.  Podremos  aplicarlos  a  procesos  reales  en  tanto  ellos  se 
aproximen  a  procesos  cuaslestátlcos  (x«t). 

Si  un  sistema  simple  se  encuentra  en  equilibrio,  no  evolucionará  espontáneamente.  Para 
que  tenga  lugar  un  proceso  necesitamos  perturbarlo,  sacarlo  del  equilibrio.  Ello  Implica  abrirlo  a 
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la  interacción  con  un  agente  externo  o  al  medio  en  general.  Concretamente  debemos 
considerar  un  sistema  compuesto,  una  o  varias  de  cuyas  ligaduras  internas  se  remueve  para 
provocar  la  evolución  al  equilibrio.  Estos  son  los  términos  en  que  hemos  planteado 
anteriormente  el  problema  central  de  la  termodinámica,  al  que  el  postulado  II  propone  una 
solución:  "...  al  remover  un  vínculo  interno  el  estado  de  equilibrio  que  el  sistema  alcance  será 
aquel  que  maxlmlce  S  sobre  todos  los  estados  compatibles  con  dicho  vínculo.” 

Veamos,  a  través  de  un  ejemplo,  cómo  se  aplican  estos  principios  a  un  proceso  concreto. 


Ejemplo  3.1 

Sea  un  cilindro  rígido,  cerrado  en  sus  extremos  y  aislado  adiabáticamente  que  lleva  en  su 
interior  un  pistón  móvil  y  dlatérmlco  que  lo  divide  en  dos  partes.  Iniclalmente  el  pistón  está  fijo 
en  el  centro  con  1  dm3  de  gas  Ideal  a  300K  y  2atm  de  presión  del  lado  Izquierdo  y  Idm3  de  gas 
Ideal  a  300K  y  1  atm  del  lado  derecho.  El  pistón  se  libera  alcanzando  su  estado  de  equilibrio  en 
una  nueva  posición.  Calcular  la  temperatura  y  presión  finales  y  el  aumento  de  entropía.  Discutir 
la  ¡rreversibilidad  del  proceso. 

El  sistema  compuesto  se  encuentra  aislado,  a  U,V,N  constantes  (no  hay  aportes  de  calor  ni 
de  trabajo  que  modifiquen  su  energía  interna).  El  nuevo  estado  de  equilibrio  será  el  que 
maxlmlce  S^+S2  con  U1+U2=cte  y  V1+V2=cte,  N1=cte  y  N2=cte.  Ese  equilibrio,  térmico  y 
mecánico,  queda  asegurado  por  la  Igualdad  de  temperaturas  y  presiones  entre  ambas  partes, 
como  se  desprende  de  las  condiciones  para  S  máximo:  1/T1  =  1/T2  y  pVt1  =  P2/T2.  Estas  dos 
Igualdades  son  ecuaciones  para  U1  y  V1  y  permiten  en  principio  resolver  el  problema. 

Con  la  condición  adicional  de  que  se  trata  de  un  gas  Ideal,  a  partir  de  la  ecuación 
fundamental  y/o  de  las  consecuentes  ecuaciones  de  estado  (U  =  3NRT/2  y  PV=  NRT)  resulta 
inmediatamente: 

T1  =  T2  =  T0 
P1  =  P2  =  3P0/2 
V1  =  2  V2  =  4V0/3. 


Se  sugiere  representar  en  un  diagrama  (P,V)  los  estados  inicial  y  final  de  ambos 
subsistemas  y  analizar  los  posibles  procesos  que  los  unen.  Internamente  ha  habido 
intercambios  de  energía  en  forma  de  trabajo  y  calor  entre  ambas  partes.  El  gas  de  la  Izquierda 
se  expande  (tiende  a  enfriarse)  realizando  trabajo  sobre  el  de  la  derecha  que  se  comprime 
(tiende  a  calentarse).  Habrá  entonces  un  flujo  de  calor  de  derecha  a  Izquierda  a  través  del 
pistón  dlatérmlco  que  compensa  la  transferencia  de  energía  en  forma  de  trabajo.  En  el  estado 
final,  de  equilibrio,  ambos  gases  tendrán  la  misma  energía  interna  que  al  comienzo  ya  que  su 
temperatura  no  habrá  cambiado. 

Calculemos  finalmente  el  aumento  de  entropía  asociado  al  proceso.  Según  la  ecuación 
fundamental  el  cambio  en  la  entropía  de  cada  parte  viene  dado  por:  AS1  =  NR.In(V'/V0)  y 
entonces: 
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AS  =  AS^aS2  =  N'R  ln(4/3)+  N2R  ln(2/3)  =  (P0V0/T0)  ln(32/27) 


La  entropía  del  sistema  Izquierdo  aumenta  y  la  del  derecho  disminuye,  resultando  un 
aumento  neto.  U,V  y  N  permanecieron  constante  pero  S  aumentó  como  esperamos  de  un 
proceso  irreversible  en  un  sistema  aislado.  Notemos  que  la  entropía  ha  aumentado  sin  que 
hayamos  podido  entregar  calor  al  sistema,  resultando  así  un  proceso  ineficiente  desde  un 
punto  de  vista  práctico. 

SI  el  proceso  del  ejemplo  se  realizara  cuaslestátlcamente,  realizando  trabajo  sobre  el  medio, 
entonces  AU<0  y  los  gases  se  habrían  enfriado:  AT<0.  En  tanto  S  no  habría  aumentado, 
llegando  así  a  un  estado  final  diferente.  Para  llegar  al  mismo  estado  final  que  antes  deberíamos 
entregar  esa  misma  cantidad  de  energía  en  forma  de  calor  para  aumentar  S  (ver  figura  3.1).  Es 
decir  que  de  esta  manera  (proceso  cuaslestátlco)  aprovecharíamos  el  aumento  de  entropía 
para  deshacernos  de  la  mayor  cantidad  posible  de  calor:  Q  <  T(AS)UiVN  (donde  el  signo  Igual 
corresponde  al  proceso  cuaslestátlco).  Cualquier  otro  proceso,  será,  menos  ventajoso  desde 
este  punto  de  vista.  Reencontramos  así  la  desigualdad. 

dQ  <  TdS 
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Potencial  de  trabajo  o  energía  libre 


Analicemos  ahora  un  proceso  a  S,V  y  N¡  constantes  para  ver  que  en  ese  caso  U(S,V,N¡) 
puede  considerarse  como  un  potencial  de  trabajo  o  energía  utilizable :  AL)  mide  el  trabajo 
externo  mínimo  necesario  para  realizar  el  proceso  o  el  máximo  que  puede  aprovecharse  en  el 
proceso  inverso  (-AU),  es  decir 


W  >  (AU)s,v,n. 


(Hemos  simplificado  la  notación  usando  N  para  representar  al  conjunto  de  los  N¡). 

Consideremos  un  sistema  compuesto  en  el  cual,  por  remoción  de  un  vínculo  interno,  ocurre 
un  proceso  isentrópico  que  lo  lleva  a  un  nuevo  estado  de  equilibrio  con  el  consiguiente  cambio 
en  la  energía  interna  AU<0.  Para  ello  acoplamos  nuestro  sistema  a  un  foco  reversible  de 
trabajo.2 

Supongamos  para  fijar  ideas,  un  sistema  como  el  del  ejemplo  1  en  el  que  S  será  constante 
si  los  procesos  son  cuaslestátlcos.  Para  ello  podemos  hacer  que  el  sistema  realice  trabajo 
sobre  un  agente  externo  (o  que  éste  haga  trabajo  sobre  aquél)  mediante  una  cuerda  sujeta  a  la 
pared  móvil  y  atravesando  el  extremo  del  recinto  de  la  Izquierda  y  una  polea,  con  una  masa 
sometida  al  campo  gravltatorlo  (foco  reversible  de  trabajo).  Para  realizar  el  proceso 
cuaslestátlcamente  podemos  suponer  la  masa  dividida  en  porciones  infinitesimales  que  pueden 
Irse  quitando  o  agregando. 

En  las  condiciones  dadas  la  reducción  de  energía  Interna  del  sistema  será  íntegramente 
aprovechada  para  realizar  trabajo  sobre  el  foco  y  en  el  proceso  inverso  todo  el  trabajo  realizado 
por  el  foco  será  íntegramente  convertido  en  energía  interna  del  sistema.  Para  un  proceso  finito 
tendríamos:  W  =  (AU)SiV,n 

En  un  proceso  real  (espontáneo)  no  se  podrá  aprovechar  todo  el  cambio  en  la  energía  como 
trabajo.  Habrá  una  agitación  que  calentará  el  gas  y  si  queremos  mantener  constante  S  será 
necesario  permitir  la  liberación  de  Q  al  exterior  aumentando  la  entropía  del  medio.  Entonces  no 
toda  la  reducción  en  energía  interna  será  aprovechable  para  realizar  trabajo  sobre  el  medio,  ya 
que  en  parte  será  disipada  en  forma  de  calor.  En  el  proceso  inverso,  no  todo  el  trabajo  externo 
podrá  convertirse  en  energía  interna  y  entonces  W  >  (AU)s,v,n,  o  en  términos  del  trabajo  externo 
-W  <  (-AU)s,v,n-  En  un  proceso  infinitesimal: 


dW  >  dll 

que  es  otra  forma  de  la  desigualdad  Implícita  en  el  principio  extremal. 


2  sistema  ideal  confinado  por  paredes  adiabáticas  e  impermeables  con  tiempos  de  relajación  tales  que  todos  los 
procesos  que  en  él  ocurren  son  cuasiestáticos  y  que  actuará  como  foco  o  sumidero  de  trabajo  cuasiestático 
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La  energía  interna  juega  entonces  el  papel  de  potencial  de  trabajo  en  procesos  a  S,V,N 
constantes,  determinando  el  máximo  trabajo  que  puede  aprovecharse  (energía  libre  o 
utillzable)  o  el  trabajo  mínimo  que  hay  que  invertir  para  llevar  a  cabo  el  proceso  inverso. 
Notemos  que  en  el  caso  trivial  en  que  no  se  intercambia  trabajo  con  ningún  foco  (W=0)  resulta: 
(AU)s,v,n  <  0  es  decir,  en  los  procesos  espontáneos  la  energía  disminuye.  Ello  es  compatible 
con  nuestra  descripción  de  los  estados  de  equilibrio  (estables)  como  estados  de  mínima 
energía  interna.  Si  estos  son  estados  de  mínima  energía  el  sistema  no  podrá  abandonarlos 
espontáneamente  y  si  ello  ocurriera  por  alguna  causa  (perturbación  externa  o  fluctuaciones)  el 
sistema  regresará  espontáneamente  a  ellos. 

En  el  siguiente  Capitulo,  luego  de  introducir  los  potenciales  termodinámlcos,  veremos  cómo 
cada  uno  de  ellos  resulta  un  potencial  de  trabajo  para  procesos  realizados  en  diferentes 
condiciones  específicas. 

Para  finalizar,  como  aplicación  de  estos  conceptos,  consideraremos  el  siguiente  proceso 
cuasiestático  y  analizaremos  los  cambios  que  ocurren  en  el  sistema  y  en  sus  partes. 


Ejemplo  3.2 

Sea  un  cilindro  rígido,  cerrado  en  sus  extremos  y  aislado  adiabáticamente  que  lleva  en  su 
interior  un  pistón  móvil  y  adiabático  que  lo  divide  en  dos  partes.  Inicialmente  ambas  partes 
tienen  volumen  V0  y  se  encuentran  llenas  con  un  gas  ideal  (Cv  =  2R)  a  T0  y  P0.  Por  medio  de 
una  resistencia  eléctrica  inserta  en  la  parte  izquierda  se  calienta  lentamente  el  gas  de  ese  lado, 
haciendo  circular  corriente,  hasta  que  se  alcanza  una  presión  27P0/8.  Hallar: 
el  estado  final  de  ambas  partes: 

-  del  lado  derecho:  AS2  =  0  (proc.  cuasiestático  adiabático);  V2=  4Vq/9;  T2=  3T J2 

-  del  lado  izquierdo:  V1  =  14V0/9;  T1  =  21T0/4 

los  aportes  de  W  y  Q  sobre  ambas  partes: 

-  trabajo  sobre  el  lado  derecho  W  =  AU2  =  NRT0 

-  trabajo  sobre  el  lado  izquierdo  \N’  =  -\N  -  -  NRT0 

-  calor  transferido  al  lado  izquierdo  Q  =  AU1  -  W’  =  AU1  +  AU2  =  19NRTo/2 

el  cambio  en  la  entropía: 


AS  =  AS1  pues  AS2  =  0 

AS1  =  N'R  ln[(T/T0)2(V/V0)]  -  (P0V0/T0)  ln(7/2)3 

Sugerencia:  discutir  los  cambios  energéticos  y  entrópicos  en  cada  parte,  en  el  sistema 
compuesto,  sin  y  con  calefactor,  y  en  el  universo. 
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CAPÍTULO  4 

Potenciales  termodinámicos 


En  este  Capítulo  introduciremos  los  potenciales  termodinámicos  mediante  cambios  en  las 
variables  independientes  para  adecuarnos  a  las  diversas  condiciones  en  que  podemos 
encontrar  al  sistema.  Para  conservar  la  información  contenida  en  la  ecuación  fundamental 
utilizaremos  la  transformada  de  Legendre. 


Formulaciones  alternativas 

Tanto  en  la  representación  entrópica  como  en  la  energética,  empleadas  hasta  ahora, 
describimos  los  estados  de  equilibrio  exclusivamente  en  términos  de  variables  extensivas,  en 
general  difíciles  de  controlar  o  medir.  La  experiencia  indica  que  son  las  variables  intensivas: 
temperatura,  presión,  etc,  más  fáciles  de  controlar  o  medir  y  más  apropiadas  entonces  para 
hacer  una  descripción  de  los  estados  de  equilibrio. 

Por  otra  parte,  en  las  representaciones  empleadas  el  principio  extremal  resulta 
adecuado  para  sistemas  aislados  con  los  que  no  es  frecuente  encontrarse.  Más  bien  nos 
interesan  sistemas  abiertos,  total  o  parcialmente,  a  la  interacción  con  el  medio.  En  esos 
casos  uno  debe  aplicar  los  principios  enunciados  al  gran  sistema,  aislado  en  su  conjunto, 
constituido  por  el  sistema  de  interés  y  el  medio  que  lo  rodea  y  de  ahí  obtener 
conclusiones  para  el  sistema  en  cuestión. 

Nos  interesa  entonces,  considerar  descripciones  alternativas  que  empleen  variables 
intensivas  para  describir  los  estados  de  equilibrio  y  que,  a  su  vez,  se  adecúen  naturalmente  a 
la  descripción  de  sistemas  abiertos.  Se  trata,  pues,  de  intentar  reemplazar  una  o  varias  de  las 
variables  extensivas  por  sus  conjugadas  intensivas  es  decir  por  las  derivadas  primeras  de  la 
relación  fundamental.  Esto  debe  hacerse  con  el  rigor  adecuado,  a  través  de  las 
transformaciones  de  Legendre,  de  modo  de  evitar  la  pérdida  de  parte  de  la  información 
contenida  en  aquella. 
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Sugerencia:  analizar  el  caso  del  gas  ideal  y  la  dependencia  de  S  con  U  o  con  T;  verificar 
que  una  vez  introducida  T  como  variable  no  puede  obtenerse  U  sino  a  menos  de  una 
constante  de  integración. 

Consideremos  el  caso  sencillo  de  una  función  de  una  variable  X,  U=U(X),  con  su  variable 
conjugada  Y=dll/dX,  tal  que  dll=YdX.  Intentamos  obtener  una  función  <t>=cj>(Y)  con  Y  como 
variable  independiente  y  el  mismo  contenido  en  información.  SI  Y=Y(X),  es  decir  si  no  se  anula 
la  derivada  segunda:  d2U/dX2  =dY/dX,  podemos  invertir  la  relación  entre  Y  y  X  y  obtener  una 
U=U(X(Y))  donde  Y  sea  la  variable  independiente.  En  el  proceso  habremos  perdido  parte  de  la 
información  ya  que  esta  relación  describe  toda  una  familia  de  curvas  con  el  mismo  valor  de  la 
pendiente  para  el  dado  valor  de  U  (figura  4.1).  Nos  falta  Incorporar  entonces,  la  información 
sobre  la  ordenada  al  origen  4>  para  poder  especificar  por  completo  la  curva  de  partida. 


pendiente  Y 


La  misma  curva  descripta  por  los  pares  de  puntos  (X,U)  que  satisfacen  U=U(X)  puede 
describirse  por  la  envolvente  de  una  familia  de  rectas  con  pendientes  Y  y  ordenadas  al  origen 
<t>,  que  satisfagan  entre  ellas  una  relación  ®=<I>(Y)  apropiada.  Gráficamente  es  fácil  comprobar 
que  Y=(U-<t>)/X,  de  modo  que  tal  relación  viene  dada  por  ®=U-XY,  que  constituye  la  llamada 
transformada  de  Legendre  de  U.  Esta  debe  entenderse  como:  ®=®(Y)=U(X(Y))-X(Y)Y  y  podrá 
obtenerse  siempre  que  d2U/dX2^0.  Además  notamos  que  d®=dll-XdY-YdX=  -X.dY  de  donde 
resulta  que  puede  recuperarse  la  variable  original  X  como  conjugada  de  la  nueva  variable  Y, 
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derivando  la  función  obtenida,  es  decir:  X=-d®/dY.  También  puede  pensarse  en  una 
transformación  inversa:  U=®+YX. 

Por  ejemplo,  para  la  función  f(x)=ax2,  la  variable  conjugada  de  x  es  y=df/dx=2ax  y  entonces 
x=y/2a.  La  transformada  de  Legendre  será:  ®(y)=f-xy=  ax2-2ax2=-  ax2=-y2/4a  (y  no  y2/4a  como 
se  obtendría  por  substitución  directa).  Verificamos  que  puede  obtenerse  x  =  -d®/dy  =  y/2a. 

La  generalización  al  caso  de  funciones  de  varias  variables  U(Xi,  ...,  Xn)  es  inmediata, 
pudiendo  reemplazarse  una  o  varias  de  dichas  variables.  En  general  se  utiliza  la  notación: 
U[Yi,...,Yk]  =  ®(Y1,...,Yk,Xk+1,...,Xn)  =  U-X1Y1-...-XkYk  para  indicar  la  transformada  de  Legendre 
de  U  en  la  que  se  han  sustituido  las  variables  extensivas  (Xi,  ...,  Xk)  por  sus  conjugadas 
intensivas  (Yi,  ...,  Xk). 


Transformadas  de  Legendre  de  la  energía  interna 


Apliquemos  las  anteriores  consideraciones  al  caso  de  la  energía  interna  de  un  sistema 
simple  U  =  U(S,V,N),  para  obtener  los  diferentes  potenciales  termodinámicos  y  sus  relaciones 
fundamentales  correspondientes.  Por  simplicidad  consideramos  el  caso  de  un  componente. 


Potencial  de  Helmholtz  F(T,V,N) 

Se  define  este  potencial  como  la  transformada  de  Legendre  de  la  energía  interna  en  la  que 
se  reemplaza  la  entropía  por  la  temperatura  como  variable  independiente,  es  decir:  F=U[T]. 
También  se  lo  denomina  energía  libre  de  Helmholtz.  En  el  cuadro  que  sigue  se  detallan  el 
cambio  de  variable,  la  ecuación  fundamental,  su  forma  diferencial  (ecuación  de  Gibbs)  y  las 
derivadas  primeras  del  potencial  en  cuestión.  Resaltamos  el  hecho  de  poder  obtener  S,  ahora 
como  variable  conjugada  de  T,  a  partir  de  la  derivada  primera  de  la  nueva  función. 


S  <-»  T=(aU/5S)v,N 
F  =  U-TS  =  F(T,V,N) 


-S  =  (3F/<9T)v,n 
-P  =  (5F/aV)T,N 


dF  =  -SdT-PdV+pdN 


|i  =  (0F/0N)TiV 


Entalpia  H(S,P,N) 

Se  define  este  potencial  como  la  transformada  de  Legendre  de  la  energía  interna  en  la  que 
se  reemplaza  el  volumen  por  la  presión  como  variable  independiente,  es  decir:  H=U[P]. 
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V  <->  -P=(5U/5V)s,n 
H  =  U+PV  =  H(S,P,N) 
dH  =  TdS+VdP+pdN 


T  =  (SH/dS)P,N 
V  =  (5H/5P)s,n 
|í  =  (5H/3N)s,p 


Potencial  de  Gibbs  G(T,P,N) 


Se  define  este  potencial  como  la  transformada  de  Legendre  de  la  energía  Interna  en  la  que 
simultáneamente  se  reemplazan  la  entropía  y  el  volumen  por  la  temperatura  y  la  presión  como 
variables  independientes,  es  decir:  G=U[T,P].  También  se  lo  denomina  energía  libre  de  Gibbs. 
Notemos  que  G  =  H-TS  razón  por  la  que  también  se  lo  denomina  entalpia  libre  y  que  en  virtud 
de  la  ecuación  de  Euler  (U=TS-PV-pN)  resulta:  G  =  Np  (S¡N¡pi¡  en  sistemas  multicomponentes). 


S  o  T=(5U/5S)v,n 
V  -P=(5U/5V)SiN 
G  =  U-TS+PV  =  G(T,P,N) 
dG  =  -SdT+VdP+pdN 


-S  s  (5G/5T)p,n 
V  =  (5G/5P)t,n 
p  =  (5G/3N)T,p 


Se  atribuye  a  Max  Born  un  esquema  mnemotécnlco  consistente  en  un  cuadrado  con  los 
potenciales  U(Energía),F,G,  H  en  los  lados  y  las  variables  S,V,T,  -P  en  los  vértices  (figura  4,2). 
Volveremos  a  usarlo  para  recordar  las  relaciones  de  Maxwell  en  el  capítulo  siguiente. 


Gran  Potencial  Q  (T,V,|u) 

Se  define  este  potencial  como  la  transformada  de  Legendre  de  la  energía  interna  en  la  que 
simultáneamente  se  reemplazan  la  entropía  y  el  número  de  moles  por  la  temperatura  y  el 
potencial  químico  como  variables  independientes,  es  decir:  Q=U[T,p].  En  virtud  de  la  ecuación 
de  Euler:  Q  =  -PV. 
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S  <-»  T=(au/5S)V,N 
N  <-»  p=(dU/dN)s,v 
Q  =  U-TS-pN  =  Q(T,V,  p) 
dQ  =  -SdT-PdV-Ndp 


-S  s=  (oQ/oT)v,  M 
-P  =  (oQ/oV)t,  „ 
-N  =  (cO/c>|i)T  ,v 


Es  importante  destacar  que  esta  necesidad  de  cambiar  variables  no  es  exclusiva  de  la 
Termodinámica.  En  Mecánica  Estadística  se  describe  un  sistema  aislado  con  el  conjunto 
estadístico  microcanónico  donde  vale  la  equiprobabilidad  de  los  microestados  y  cuya  propiedad 
macroscópica  asociada  es  S.  Para  describir  sistemas  abiertos,  de  energía  y  composición 
variable,  se  introducen  sucesivamente  los  conjuntos  estadísticos  canónico  y  gran  canónico, 
que  se  vinculan  con  las  propiedades  macroscópicas  FyQ,  respectivamente. 

Observemos  que  estos  potenciales,  derivados  de  la  energía,  son  propiedades  extensivas  y 
que,  al  depender  de  variables  independientes  tanto  extensivas  como  intensivas,  tienen 
propiedades  de  homogeneidad  híbridas.  Por  ejemplo,  para  F(T,V,N):  F(T,¡LV,XN)  =  ^F(T,V,N). 
Este  tipo  de  relaciones  se  denominan  de  homogeneidad  generalizada'.  f(Xpx,Xqy)  =  ?T(x,y). 

Veamos  finalmente  que  si  intentáramos  cambiar  todas  las  variables  obtendríamos  un 
potencial  trivial  idénticamente  nulo.  En  efecto:  Z(T,P,p)  =  U [T , P , pi]  =  U-TS+PV-pN  =0  en  virtud 
de  la  ecuación  de  Euler.  Ello  está  de  acuerdo  con  la  necesidad  de  incluir,  en  la  descripción  de 
un  sistema  macroscópico,  al  menos  una  variable  extensiva  para  dar  cuenta  de  la  dependencia 
con  el  tamaño  del  sistema. 


Transformadas  de  Legendre  de  la  entropía 


Se  las  denomina  funciones  de  Massieu  y  se  obtienen  reemplazando  una  o  varias  variables 
extensivas  por  sus  respectivas  conjugadas  intensivas.  En  general  pueden  relacionarse  con  los 
potenciales  ya  vistos: 

M-i  =  S[1/T]  =  S-(1/T)U  =  -(U-TS)/T  =  -F/T 

M2  =  S[1/T,P/T]  =  S-(1/T)U-(P/T)V  =  -(U-TS+PV)/T  =  -G/T 

M3  =  S[1/T,  p/T]  =  S-(1/T)U+(p/T)N  =  -(U-TS-pN)/T  =  -G/T. 

Las  fundones  de  Massieu,  al  igual  que  los  potenciales  derivados  de  la  energía,  son 
propiedades  extensivas  y  satisfacen  relaciones  de  homogeneidad  generalizada. 

Sugerencia:  analizar  la  aplicación  de  estos  conceptos  a  la  ecuación  fundamental  para  la 
radiación  del  cuerpo  negro. 
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El  principio  extremal  para  los  potenciales 


La  utilidad  de  estos  potenciales  para  su  aplicación  a  las  diferentes  condiciones  en  que 
puede  encontrarse  un  sistema  puede  concretarse  reformulando  el  principio  extremal.  Nos 
basaremos,  entonces,  en  el  principio  de  mínimo  para  un  sistema  aislado:  “El  valor  de  equilibrio 
de  cualquier  parámetro  interno  no  restringido  es  el  que  minimiza  la  energía  interna  para  los 
dados  valores  de  la  entropía,  el  volumen  y  los  números  de  moles  totales”. 


mínimo  Ut 

s+s0  =  St 
v+v0  =  vt 

N+N0  =  N, 


Consideremos  ahora  un  sistema  en  contacto  térmico  con  un  reservoño  de  calor  3  o  baño 
térmico  que  actúa  como  fuente  o  sumidero  de  calor  cuasiestático  sin  alterar  su  propia 
temperatura  T0  y  por  lo  tanto  la  del  sistema  (ver  figura).  Nos  proponemos  establecer  un 
principio  para  predecir  el  estado  de  equilibrio  que  alcanzará  el  sistema  al  remover  en  él  algún 
vínculo  interno.  Suponemos  que  el  conjunto  formado  por  el  sistema  y  el  reservorio  puede 
considerarse  aislado,  es  decir,  mantenido  con  valores  de  St,Vt  y  Nt  constantes.  En  esas 
condiciones  los  valores  de  equilibrio  de  los  parámetros  internos  serán  los  que  minimicen  Ut  = 
U+U0  de  modo  que: 

Ut  es  mínimo  =>  dUt  =  0  y  d2Ut  >  0 
con  las  condiciones  de  cierre: 

dS,  =  d(S+S0)  =  0  =>  dS0  =  -dS 


3  sistema  confinado  por  paredes  rígidas  e  impermeables  con  tiempos  de  relajación  tales  que  todos  los  procesos  que  en 
él  ocurren  son  esencialmente  cuasiestáticos  y  suficientemente  grande  para  que  el  intercambio  de  calor  no  afecte  su 
temperatura  (9T/9U)  ->  0 
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y  condiciones  semejantes  para  Vt  y  Nt.  Pero,  tanto  el  sistema  como  el  reservorio  están 
rodeados  de  paredes  rígidas  e  impermeables  y  entonces  todas  estas  variables  están  fijas: 
dV=dN=dV0=  dNo=0. 

Las  condiciones  de  equilibrio  pueden  transformarse  como  sigue: 

dUt  =  dU+dll0  =  dll+T0dS0  =  dU-T0dS  =  d(U-T0S)  =  0 
d2Ut  =  d2U+d2U0  =  d2U+(5To/aS0)(dS  0)2  =d2U  >  0 

pues  en  un  reservorio  las  variables  extensivas  V0  y  N0  son  constantes  de  modo  que  todos  los 
términos  que  incluyen  dV0  o  dN0  se  anulan  y  32U0/3S02  =  dTJdS0  =  0  para  un  reservorio.  Por 
otro  lado,  por  ser  T0  una  constante  y  S  una  variable  independiente: 

d2Ut  =  d2U-d2(T0S)  =  d2(U-T0S)  >  0 

Así  resulta  que  las  condiciones  de  equilibrio,  d(U-T0S)  -  0  y  d2(U-T0S)  >  0,  son  las 
requeridas  para  minimizar  (U-T0S)  que,  teniendo  en  cuenta  que  la  temperatura  del  sistema  está 
fija  en  T=T0,  coincide  con  el  potencial  de  Helmholtz  del  mismo.  En  otras  palabras,  minimizar  UT 
para  el  sistema  más  el  baño  térmico  equivale  a  minimizar  F  para  el  sistema  en  cuestión,  sobre 
el  conjunto  de  estados  con  T=T0. 

Entonces,  el  valor  de  equilibrio  de  los  parámetros  internos  libres  en  un  sistema  en  contacto 
térmico  con  un  reservorio  de  calor  es  aquel  que  minimiza  la  energía  libre  de  Helmholtz  para 
T=T0  y  los  dados  valores  de  V y  N. 

El  significado  intuitivo  del  razonamiento  anterior  es  el  siguiente.  Ya  que  d(-TS)  representa  el 
cambio  en  la  energía  del  reservorio:  -T0dS  =  T0dS0,  decir  que  es  mínima  la  energía  Interna  total 
es  equivalente  a  decir  que  es  mínima  la  energía  libre  F  del  sistema. 

Análogamente  podemos  decir  que  los  estados  de  equilibrio  de  un  sistema  son  los  de 

-  mínimo  H  para  los  dados  valores  de  S,  p  y  N. 

-  mínimo  G  para  los  dados  valores  de  T,  p  y  N. 

-  mínimo  Q  para  los  dados  valores  de  T,  V  y  p. 


Los  potenciales  termodinámicos  como  potenciales  de  trabajo 


Como  vimos,  la  energía  interna  U  juega  el  papel  de  potencial  de  trabajo  en  procesos  a 
S,V,N  constantes  determinando  el  trabajo  mínimo  que  hay  que  invertir  para  llevar  a  cabo  un 
proceso  o  el  máximo  trabajo  que  puede  aprovecharse  ( energía  libre  o  utlllzable)  en  el  proceso 
inverso  inverso  (-AU,  en  este  caso).  Notamos  también  que  en  el  caso  trivial  en  que  no  se 
intercambia  trabajo  con  ningún  foco  (W=0)  resulta  (AU)s,v,n  ^  0,  es  decir,  en  los  procesos 
espontáneos  la  energía  disminuye.  Ello  es  compatible  con  nuestra  descripción  de  los  estados 
de  equilibrio  (estables)  como  estados  de  mínima  energía  interna. 
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Mediante  razonamientos  análogos  para  los  potenciales  termodlnámicos,  podemos  ver  como 
cada  uno  de  ellos  resulta  un  potencial  de  trabajo  para  procesos  realizados  en  diferentes 
condiciones  específicas.  Su  cambio,  A®,  da  la  máxima  energía  aprovechable  en  forma  de 
trabajo  o  “energía  libre”en  procesos  cuaslestátlcos  conducidos  bajo  diferentes  condiciones.  Por 
ejemplo  puede  verse  que 


W  >  (AF)t,v,n  ■ 


Es  decir  que  el  trabajo  mínimo  necesario  para  un  proceso  a  T,V,N  constantes  viene  dado 
por  el  valor  de  AF  y  corresponde  a  un  proceso  cuasiestático;  en  cualquier  otro  caso  habrá  que 
invertir  una  cantidad  mayor  de  trabajo.  Viceversa,  la  máxima  energía  aprovechable  como 
trabajo,  energía  libre,  en  un  proceso  a  T,V,N  constantes  será  -AF.  Por  las  propiedades  de  las 
fundones  de  estado,  no  es  necesario  que  los  procesos  sean  estrictamente  a  T,V,N  constantes 
sino  que  es  suficiente  con  que  estas  variables  tengan  el  mismo  valor  al  comienzo  y  al  final  del 
proceso  en  cuestión. 

Asimismo  podemos  escribir: 

W  >  (AH)s,p,n  , 

W  >  (AG)t,p,n  . 

W  >  (AQ)T  V,p , 

con  significado  semejante. 


Ejemplo  4.1 

Sea  un  cilindro  rígido,  con  un  pistón  interno  a  cada  lado  del  cual  hay  1  mol  de  un  gas  ideal 
monoatómico.  Sus  paredes  son  diatérmicas  y  el  sistema  está  sumergido  en  un  baño  térmico  a 
0°C,  de  gran  capacidad.  Inicialmente  V1,  =  10  I  y  V2¡  =  1  I  y  luego  el  pistón  se  desplaza 
cuasiestáticamente  hasta  V\  =  6  I  y  V2f  =  5  I.  Hallar  el  trabajo  aprovechable  usando  el  potencial 
de  Helmholtz. 

Este  trabajo  vendrá  dado  por  -AF  del  sistema  compuesto. 

Obtenemos  F(T,V,N)  a  partir  de  M (1/T,V,N)  =  S[1/T]  =  -F(T,V,N)/T: 

como  vimos,  para  un  gas  ideal:  S  =  (N/N0)S0  +  NR.In[(U/U0)3/2(V/V0)(N/N0)_5/2] 
de  donde:  M-,  =  S  -  (1/T)U  =  (N/N0)S0  +  NR.In[(T/T0)3/2(V/V0)(N/N0)_1]  -  (3/2)NR 
y  entonces:  F  =  -  TM,  =  (NT/N0T0)F0  -  NRT.In[(T/T0)3/2(V/V0)(N/N0)-1]. 
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En  procesos  a  T  y  N  constantes,  AF  =  -NRT.  ln(Vf/V¡) 


y  entonces:  AF  =  AF(1)  +  AF(2)  =  -RT.In(6/10)  -RT.In(5/1)  =  -RT.In(3)  =  -2495  J. 

Todo  el  trabajo  obtenido  (-AF)  proviene  de  la  fuente  térmica  ya  que  la  energía  interna  del 
sistema  no  ha  cambiado  (AU  =  AT  =  0).  Notemos  que  en  general:  AF  =  AU-T.AS  y  la 
disminución  de  F  (energía  libre)  puede  tener  origen  en  una  disminución  de  la  energía  U  y/o  en 
un  aumento  de  la  entropía  S.  En  el  caso  presente  la  disminución  de  F  se  debe  exclusivamente 
al  aumento  de  entropía. 

Sugerencia:  analizar  el  proceso  espontáneo  que  tiene  lugar  cuando  se  libera  el  pistón  sin 
realizar  trabajo  sobre  un  agente  externo. 
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CAPÍTULO  5 

Funciones  respuesta 


En  este  capítulo  introduciremos  algunas  funciones  respuesta,  habitualmente  accesibles  al 
experimento,  y  las  relacionaremos  con  las  derivadas  segundas  de  los  potenciales 
termodlnámlcos.  Analizaremos  las  relaciones  entre  ellas  y  discutiremos  métodos  para 
reducirlas  a  un  conjunto  base. 


Funciones  respuesta  y  derivadas  segundas 

Ya  vimos  que  las  derivadas  primeras  de  U  (o  de  S)  tienen  sentido  físico.  De  hecho, 
introducen  las  propiedades  intensivas  T,P,|i  en  nuestro  formalismo  y  como  tales  aparecen  en 
las  condiciones  de  equilibrio.  Asimismo,  las  derivadas  primeras  de  los  potenciales 
termodlnámlcos  asocian  con  cada  variable  (intensiva  o  extensiva)  su  respectiva  variable 
conjugada  (extensiva  o  intensiva). 

También  las  derivadas  segundas  tienen  sentido  físico  pues  se  relacionan  con  magnitudes 
físicas  medibles,  características  de  las  sustancias:  las  funciones  respuesta.  Ellas  miden  cómo 
cambia  una  variable  de  estado  al  vahar  independientemente  otra,  bajo  condiciones  dadas; 
siendo  así  las  cantidades  termodinámicas  generalmente  más  accesibles  al  experimento.  Como 
derivadas  de  variables  intensivas  respecto  de  extensivas  (o  viceversa)  se  relacionan 
formalmente  con  derivadas  segundas  de  los  potenciales  termodlnámlcos  y  como  tales  jugarán 
un  papel  importante  en  las  condiciones  de  estabilidad  que  más  adelante  discutiremos. 

Recordemos  la  definición  de  algunas  funciones  respuesta  conocidas,  todas  supuestas  a  N 
constante: 

i)  coeficiente  de  dilatación  a 

definido  como  incremento  relativo  de  volumen  por  unidad  de  aumento  de  temperatura,  a 
presión  constante: 
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a  es  (1/V).(5V/5T)P  =  (1/V).(32G/3T3P)  =  (1/V).GTP 


donde  hemos  tenido  en  cuenta  que  V  =  (0G/0P)T  y  que  la  presencia  de  T  y  P  como  variables 
independientes  sugiere  trabajar  en  la  representación  de  Gibbs. 


ii)  capacidad  calorífica  molar  a  presión  constante  CP 

definida  como  cantidad  de  calor  que  es  necesario  agregar  cuasiestáticamente  por  mol  para  un 
incremento  unitario  de  temperatura,  a  presión  constante: 

CP  ss  (1/N).(dQ/dT)P  =  (T/N).(SS/5T)P  =  -(T/N).(52G/5T2)  =  -(T/N).Gtt 

donde  nuevamente  la  presencia  de  T  y  P  como  variables  independientes  sugiere  trabajar  en  la 
representación  de  Gibbs  y  hemos  tenido  en  cuenta  que  S  =  -(6G/0T)P. 


iii)  capacidad  calorífica  molar  a  volumen  constante  Cv 

definida  como  cantidad  de  calor  que  es  necesario  agregar  cuasiestáticamente  por  mol  para  un 
incremento  unitario  de  temperatura,  a  volumen  constante: 

Cv  ss  (1/N).(dQ/dT)v  =  (T/N).(5S/5T)V=  -(T/N).(52F/5T2)  =  -(T/N).Ftt 

donde  la  presencia  de  T  y  V  como  variables  independientes  sugiere  ahora  la  representación  de 
Helmholtz,  en  la  cual  S  =  -(0F/5T)V. 


iv)  compresibilidad  isotérmica  kj 

definida  como  disminución  relativa  de  volumen  por  unidad  de  aumento  de  presión,  a 
temperatura  constante: 

KT  ss  -(1/V).(0V/5P)t  =  -(1/V).(S2G/3P2)  =  -(1/V).GPP 

donde  nuevamente  la  presencia  de  T  y  P  como  variables  sugiere  la  representación  de  Gibbs. 
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v)  módulo  volumétrico  adiabático  Bs 


definido  como  incremento  de  presión  necesario  para  producir  una  disminución  unitaria  relativa 
de  volumen  a  entropía  constante: 


Bs  =  -V.(5P/5V)S  =  V.(a2U/5V2)  =  V.UW  =  1/ks 


donde  ahora  la  presencia  de  S  y  V  como  variables  sugiere  la  representación  energética,  en  la 
cual  P  =  -(SU/SV)s. 

Vemos  pues  como  estos  coeficientes  a,  CP,  Cv,  kt,  ks  y  otros  que  podrían  definirse  se 
relacionan  con  derivadas  segundas  de  potenciales  termodlnámicos.  Más  adelante  veremos 
cómo  relacionar  estas  derivadas  segundas  entre  sí  y  expresar  cualquier  función  respuesta  en 
términos  de  un  conjunto  base. 


Relaciones  de  reciprocidad  de  Maxwell 


Por  las  condiciones  de  continuidad  de  los  potenciales  y  sus  derivadas  primeras,  surge 
inmediatamente  la  Igualdad  de  sus  derivadas  segundas  mixtas: 

para  O  =  <t>(x1,x2)  tendremos:  d2®/5x- |3x2  =  32<t >/5x23x1 

que  en  términos  de  las  derivadas  primeras: 


yi  =  (c*t/c>xi)x2  =  yi(xi,x2)  e 


y2=  (0®/0x2)x1  =  y^.Xz) 


se  traduce  en  la  siguiente  relación: 

(5y2/5x1)x2  =  (3y1/3x2)x1 

Estas  propiedades  matemáticas  se  conocen  en  termodinámica  como  relaciones  de  Maxwell 
y  están  evidentemente  ligadas  con  el  hecho  de  que  los  potenciales  O  son  funciones  de  estado, 
es  decir  que  d®  es  un  diferencial  exacto,  o  bien  que  su  integral  sobre  un  camino  cerrado  en  el 
plano  (x1;x2)  se  anula.  Su  Importancia  radica  en  que  hacen  predicciones  físicas  no  triviales, 
experimentalmente  comprobables.  Veamos  algunas  de  las  más  conocidas.  Para  un  sistema 
simple  con  N  constante  (o  {N¡}  constantes)  sólo  quedan  dos  variables  Independientes  y  por 
tanto,  una  relación  de  Maxwell  o  bien  una  por  potencial,  a  saber: 
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en  U(S,V): 

(52U/5V5S)  =  (52U/5S5V) 

como  (5U/3S)V  =  T  y  (5U /<9V)S  =  -  P 

resulta:  (3TA9V)S  =  -(3P/3S)V 

en  F(T,V) 

(52F/5V5T)  =  (32F/3T3V) 

como  (3F/3T)V  =  -S  y  (<9F/3V)T  =  -P 

resulta:  (3S/3V)T  =  (3P/3T)V 

en  G(T,P) 

(<92G/3PST)  =  (32G/3T3P) 

como  (3G/3T)P  = -S  y  (3G/3P)T  s  V 

resulta:  (3S/3P)T  =  -(oV/oT)P 

en  H(S,P) 

(32H/3P3S)  =  (32H/3S3P) 

como  (3H/3S)p  =  T  y  (3H/3P)S  =  V 

resulta:  (3T/3P)S  =  (3V/3S)P 


Estas  cuatro  relaciones,  muy  útiles,  son  en  rigor  una  sola  ya  que  a  partir  de  cualquiera  de 
ellas  pueden  obtenerse  las  restantes  cambiando  de  representación.  Para  recordarlas  puede 
usarse  la  regla  mnemotécnlca  de  M.  Born  aplicada  al  cuadrado  que  definimos  anteriormente. 
Partiendo  de  cualquier  vértice  se  recorren  los  dos  siguientes  en  cierto  sentido,  generando  una 
derivada  parcial  que  incluye  las  3  variables  en  esos  vértices.  Esta  será  Igual  a  la  que  se 
obtenga  partiendo  del  otro  vecino  y  recorriendo  otros  dos  en  sentido  contrario,  con  el  signo 
cambiado  si  cambia  la  naturaleza  de  la  variable  de  partida  (extensiva  o  Intensiva). 


Es  interesante  analizar  el  tipo  de  predicciones  físicas  no  triviales  que  estas  relaciones 
encierran.  Consideremos,  por  ejemplo,  la  relación  de  Maxwell  en  la  representación  de  Gibbs: 
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(3S/5P)t  =  -(oV/oT)P  =  -  aV.  El  término  de  la  Izquierda  es  generalmente  negativo,  ya  que  debe 
liberarse  calor  al  comprimir  un  sistema  si  se  quiere  mantener  la  temperatura  constante.  Por 
tanto  deberá  ser  a  >  0  como  se  verifica  en  la  mayoría  de  las  sustancias.  La  excepción  más 
conocida  la  constituye  el  agua,  entre  0  y  4°C,  que  se  contrae  al  calentarla  (a  <  0).  En  este 
rango  entonces  deberá  ser  (3S/0P)T  >  0,  es  decir  que  será  necesario  agregar  calor  para  evitar 
que  se  enfríe  al  comprimirla,  tal  como  puede  verificarse. 

Para  N  variable,  en  un  sistema  simple  aparecen  un  par  de  relaciones  nuevas,  que  en  la 
representación  de  Gibbs  serían: 

(3p/aT)P,N  =  (52G/5T5N)p  =  (32G/3N3T)P  =  -(oS/oN)T,P 
(5|u/5P)t,n  =  (52G/5P5N)t  =  (32G/3N3P)T  =  (3V/3N)T,P 

Y  en  sistemas  multlcomponentes,  entre  otras  relaciones,  tendríamos  algunas  del  tipo: 

(3|ij/c>N¡)Nj  =  (32G/3N¡3Nj)  =  (32G/3Nj3N¡)  =  (3|i¡/3Nj)Ni 

todo  a  T,  P  y  los  restantes  Nk  constantes.  Estas  relaciones  muestran  la  reciprocidad  en  el 
cambio  resultante  del  potencial  químico  de  un  componente  al  agregar  moles  de  otro. 


Derivadas  segundas  independientes  y  reducción  de  derivadas 


En  general,  (en  cualquier  representación)  si  tenemos  n  variables  independientes,  habrá: 
n2  derivadas  segundas,  n  de  ellas  respecto  de  la  misma  variable, 
n2-n  =  n(n-1)  derivadas  segundas  mixtas, 
n(n-1)/2  relaciones  de  Maxwell  que  las  vinculan  de  a  pares  y 
n+n(n-1)/2  =  n(n+1)/2  derivadas  segundas  Independientes. 

Entonces  cualquier  derivada  segunda  podrá  reescriblrse  en  términos  de  un  conjunto  de 
n(n+1)/2  derivadas  elegidas  como  base,  mediante  cambios  de  coordenadas  apropiados. 

En  un  sistema  simple  monocomponente  tendremos  3  variables  independientes  (ej:  S,V,N)  y 
por  lo  tanto  9  derivadas  segundas,  de  las  cuales  sólo  6  serán  independientes  en  virtud  de  las  3 
relaciones  de  Maxwell  (ej:  Uss,  Uw>  Unn,  Usv~Uvs,  Usn=Uns  y  Uvn  =Unv)- 

SI  además  podemos  tomar  N  constante,  n=2  y  sólo  quedan  tres  derivadas  independientes 
(Uss,  Uw>  USv=UVs)-  En  la  práctica,  es  conveniente  elegir  como  conjunto  base  a  GTt,Gtp  y 
GPP  dada  su  relación  directa  con  CP,  a  y  kt,  magnitudes  accesibles  al  experimento  y 
habltualmente  tabuladas. 

A  continuación  describiremos  cómo  expresar  cualquier  otra  derivada  segunda  (función 
respuesta)  en  términos  de  éstas  o  de  cualquier  otro  conjunto  elegido  como  base.  Si  bien  este 
no  es  un  tema  fundamental  resulta  bastante  útil  en  las  aplicaciones  ya  que  es  habitual  tener 
la  necesidad  de  evaluar  alguna  función  respuesta,  la  que  puede  fácilmente  expresarse  como 
derivada  segunda  y  luego  reducirse  a  una  expresión  en  términos  de  un  conjunto  de 
magnitudes  conocidas. 

En  términos  generales  esta  reducción  Implica  un  pasaje  a  la  representación  de  Gibbs  con  el 
asociado  cambio  de  coordenadas.  En  principio,  cualquier  derivada  primera  que  involucre 
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parámetros  intensivos  (extensivos)  puede  expresarse  en  términos  de  derivadas  segundas  que 
impliquen  parámetros  extensivos  (intensivos)  reemplazando  aquellos  por  su  definición.  Si 
aparecen  diferenciales  de  otros  potenciales  termodinámicos  o  del  potencial  químico  p,  pueden 
reemplazarse  usando  respectivamente  las  ecuaciones  de  Gibbs  de  los  potenciales  en  cuestión 
o  la  ecuación  de  Gibbs-Duhem.  Para  esta  reducción  es  conveniente  utilizar  propiedades 
conocidas  de  las  derivadas  segundas  como  la  relación  cíclica  (Euler)  y  algunas  otras  que 
resulta  fácil  implementar  recurriendo  a  Jacobianos  y  sus  propiedades  (ver  Anexo).  En 
ocasiones  se  debe  recurrir  a  las  relaciones  de  Maxwell.  Veamos  cómo  proceder,  a  través  de 
algunos  ejemplos. 


Ejemplo  5.1 

Hallar  el  cambio  en  S  y  en  U  en  una  compresión  isotérmica. 
dS  =  (<9S/<9P)T,N.dP  =  -  (<9V/<9T)P,N.dP  =  -aV.dP, 
donde  hemos  usado  la  relación  de  Maxwell:  (dS/SP)T  =  -(c’V/c>T)P 
dU  =  (<9U/<9P)T,N.dP  =  [T(<9S/3P)t,n  -  P(3V/5P)T.N].dP  =  (-aT+  KTP)V.dP 
Para  gases  ideales:  a=1/T,  kt  =  1/P  y  entonces  dll=0 


Ejemplo  5.2 

Demostrar  que  para  cualquier  sistema  simple:  CP-CV  =  TVo2/Nkt. 
Evaluemos  Cv  =  (T7N)(dS/<9T)v,  en  términos  de  CP,  a  y  kt. 

(<9S/<9T)V  =  3(S,V)/3(T,V)  =  [5(S,  V)/S(T,  P)]/[S(T,  V)/S(T,  P)] 

donde  introducimos  las  variables  T  y  P  usando  la  regla  de  la  cadena. 

El  denominador  es  directamente  -Vkt  y  al  numerador  conviene  expandirlo: 
3(S,V)/3(T,P)  =  (<9S/«9T)p(<9V/3P)t  -  (<9S/<9P)t(<9V/3T)p  =  (NCp/T)(-Vkt)+(ciV)2 
Entonces: 

Cv  =  -(T/N)[-NCpVkt/T  +  (oV)2]/Vkt  =  CP  -  TVo2/Nkt  c.q.d. 

Para  gases  ideales:  a=1/T,  kt  =  1/P  y  entonces  CP-CV  =  PV/NT  =  R 


Ejemplo  5.3 

Hallar  el  cambio  en  T  en  una  expansión  isentálpica: 
dT  =  (<3T/dP)H,N-dP 

(<9T/<9P)h  =  3(T,H)/3(P,H)  =  [3(T,H)/3(T,P)]/[3(P,H)/3(T,P)]  =  -(3H/3P)t/(3H/3T)p, 
que  es  precisamente  la  relación  cíclica.  A  su  vez: 

(dH/dP)T  =  [T(3S/3P)t  +  V.1]  =  -T(<9V/5T)P  +  V  =  -aVT  +  V 
(<9H/<9T)P  =  [T(<9S/«9T)P  +  V.O]  =  NCP 
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Y  finalmente 

dT  =  -[(-aVT+V)/NCP].dP  =  (v/CP).(aT-1).dP 
Para  gases  ideales:  a=1/T  y  entonces  dT=0 

La  expresión  encontrada  es  útil  al  analizar  el  experimento  de  Joule  -  Thomson  de  expansión 
adiabática  través  de  un  tabique  poroso  o  estrangulamlento  ( throttling ),  que,  en  buena 
aproximación,  puede  considerarse  ¡sentálplca.  . 


Anexo:  Jacobianos  y  derivadas  segundas 


Siendo  (u,v)  funciones  de  (x,y)  se  define  el  jacobiano  J  como  el  determinante  de  la  matriz  de 
sus  derivadas  primeras: 

J  =  3(u,v)/3(x,y)  =  (3u/3x)(3v/3y)-(3u/3y)(3v/3x) 

Sus  propiedades  más  importantes  se  derivan  de  las  de  los  determinantes  de  matrices  y  son 
útiles  en  nuestra  discusión  por  su  relación  directa  con  las  derivadas  parciales: 

(Su  /3x)y  =  3(u,y)/3(x,y). 

Recordemos  algunas  propiedades: 

1 .  permutando  filas  o  columnas: 

3(u,v)/3(x,y)  =  -  3(v,u)/3(x,y)  =  -  3(u,v)/3(y,x) 

2.  usando  la  matriz  identidad: 

3(u,v)/3(x,y)  =  1/3(x,y)/3(u,v) 

3.  especie  de  regla  de  la  cadena: 

3(u,v)/3(x,y)  =  [S(u,v)/S(r,s)].[S(r,s)/S(x,y)] 

Como  aplicación,  demostraremos  la  propiedad  cíclica  (Euler)  de  las  derivadas  parciales: 
(3x/3y)z  (3z/3x)y  (Sy/Sz)x  =  -1 

(Sx/Sy)z  =  S(x,z)/3(y,z)  =  [S(x,z)/S(x,y)].[S(x,y)/S(y,z)]  =  -[3(x,z)/3(x,y)]/[3(y,z)/3(y,x)] 

=  -(Sz/Sy)x/(Sz/Sx)y, 

que  conduce  a  la  relación  buscada. 
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CAPITULO  6 

Criterios  de  estabilidad 


En  este  capítulo  analizaremos  por  completo  las  condiciones  de  extremo,  implícitas  en  el 
principio  extremal,  estableciendo  así  ciertos  criterios  de  estabilidad.  Analizaremos  también 
cómo  estas  condiciones  influyen  en  la  forma  de  los  potenciales. 


Criterios  de  estabilidad 

Las  derivadas  segundas  de  los  potenciales  termodinámicos  introducidas  en  el  capítulo 
anterior  aparecen  en  las  condiciones  de  estabilidad  de  los  estados  de  equilibrio.  Para  un 
sistema  aislado,  es  decir,  mantenido  con  valores  constantes  de  St,Vt  y  Nt,  los  valores  de 
equilibrio  de  los  parámetros  internos  serán  los  que  minimicen  la  energía  interna  total  Ut: 

Ut  es  mínimo  =>  dUt  =  0  y  d2Ut  >  0 

con  las  condiciones  de  cierre:  dSt  =  dVt  =  dNt  =  0. 

Consideremos  dentro  de  ese  sistema,  un  subsistema  pequeño  pero  macroscópico, 
constituido  por  N  moles  de  sustancia,  dejando  N'  en  el  resto,  con  N«N'=  Nt.  El  subsistema 
puede  entonces  intercambiar  energía  y  volumen  con  el  resto  y  como  hemos  visto,  ello 
determina  las  condiciones  de  equilibrio  que  surgen  de  dll,  =  0.  En  efecto: 


dUt  =  Ndu  +  N'du'  =  N(usds  +  uvdv)  +  N'(u's-ds'  +  u'vdv')  =  0 
como  N'ds'=  -Nds  y  N'dv'=  -Ndv,  resulta 
dUt  =  N(T-T')ds  -  N(P-P')dv  =  0;  V  ds,  V  dv 
=í>  T  =  T'  y  P  =  P' 

Dada  la  arbitrariedad  en  la  elección  del  subsistema  pequeño,  estas  igualdades  garantizan  la 
uniformidad  de  T  y  P  sobre  todo  el  sistema. 

Analicemos  ahora  las  desigualdades  que  surgen  de  la  condición  de  mínimo  d2Ut  >  0: 
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d2Ut  =  Nd2u  +  N'dV  >  O 

con  d2u  =  Uss(ds)2  +  2usv(ds)(dv)  +  uw(dv)2,  pues  usv  =  uvs  (Maxwell)  y  una  expresión  similar  para 
d2u'.  Pero,  teniendo  en  cuenta  que  |ds'|  =  (N/N’)|ds|  «|ds|  y  análogamente  |dv'|«|dv|,  puede 
despreciarse  N'd2u'  frente  a  Nd2u.  Así 

d2Ut  =  Nd2u  =  N[uss(ds)2  +  2usv(ds)(dv)  +  uw(dv)2]  >  0,  Vds,  Vdv  (excepto  ds=dv=0). 

Vemos  entonces  que  d2u  debe  ser  una  forma  cuadrática  definida  positiva  que  puede 
llevarse  a  una  forma  homogénea  mediante  un  cambio  de  variables  (s,v)— >(T,v).  En  efecto  como 
T=us 

dT  =  u3Sds  +  usvdv 

(dT)2/uss  =  uS3(ds)2  +  2usv(ds)(dv)  +  (usv2/uss)(dv)2 

Sumando  y  restando  (uSv2/uSs)(dv)2  a  la  expresión  para  d2u  resulta: 

d2u  =  (1/uss)(dT)2  +  (uvv-  usv2/uss)(dv)2  >  0  VdT,  Vdv 

Pueden  establecerse  así  dos  criterios  de  estabilidad  (intrínseca): 

uss  >  0 

Uss^w- Usv  '>  0 

que  también  pueden  expresarse  en  términos  de  derivadas  de  variables  intensivas  respecto  de 
sus  extensivas  conjugadas. 

(3T/5S)V  >  0 
-(3P/SV)t  >  0 

En  efecto: 

uss  =  (32u/3s2)v  =  (3T/3s)v  >  0 

UssUw-Usv2  =  3(us,uv)/3(s,v)  =  -S(T,P)/3(s,v)  =  -[3(T,P)/3(T,v)]/[3(s,v)/3(T,v)] 

=  -(3P/3v)t/(3s/3T)v  >  0 

Asimismo,  pueden  expresarse  en  términos  de  funciones  respuesta,  ya  que: 

Cv  =  (T/N)(oS/oT)v  =  T/(3T/3s)v  >  0 
kt  =  -(1/v)/(3P/3v)T  >  0 
es  decir: 

Cy  >  0 

Kj  >  0 

Así  formulados  estos  criterios  resultan  bastante  intuitivos.  Por  ejemplo,  Cv  >  0  implica  que 
para  que  un  sistema  sea  estable  debe  aumentar  su  temperatura  al  aportársele  calor.  Por  tanto 
si  el  sistema  sufre  una  desviación  local  de  su  estado  de  equilibrio  (T<T,  por  ejemplo)  el 
proceso  espontáneo  que  tiene  lugar,  flujo  de  calor  hacia  la  zona  más  fría,  tiende  a  restablecer 
la  igualdad  de  temperaturas  (equilibrio). 
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Si  el  cambio  de  variables  para  pasar  a  una  forma  homogénea  hubiera  sido:  (s,v)— >(s,-P),  se 
habrían  obtenido  un  par  de  criterios  equivalentes:  ks  >  0  y  CP  >  0  que  podrían  derivarse  de  los 
anteriores.  Asimismo  si  el  subsistema  se  hubiera  elegido,  por  ejemplo,  con  V  fijo  y  N  variable, 
se  habrían  obtenido  otros  criterios,  en  términos  de  derivadas  de  T  y  p  respecto  de  S  y  N,  que 
no  son  independientes  de  los  anteriores. 

En  general  para  un  sistema  con  n  variables  (n  =  r+2  en  un  sistema  simple  de  r 
componentes)  habrá  n-1  criterios  de  estabilidad  independientes,  que  pueden  escribirse  de 
manera  compacta: 


(dY¡/5X¡)cond  >  0  cond  =  Yi,  ...,  YM,  X¡+1> .... Xn  constantes; 

para  i  =  1,2,...,n-1.  Notemos  que  la  n-ésima  derivada  de  este  tipo,  (6Yn/6Xn),  resulta 
idénticamente  nula  al  tener  fijas  las  n-1  variables  intensivas  restantes. 


Principio  de  Le  Chatelier 


Nuestros  criterios  establecen  que  deben  ser  positivas  las  derivadas  del  tipo  (dY^SX-i),  tanto 
con  X2  constante  como  con  Y2  constante.  Es  fácil  comprobar  que  estas  derivadas,  a  su  vez, 
deben  satisfacer  una  desigualdad  entre  ellas: 

(3Y1/aXi)X2  >  (3Y1/aX1)Y2  >  0 

conocida  como  Principio  de  Le  Chatelier.  En  efecto: 


(5Y1/ax1)Y2  =  d(y1y2)id(x,y2)  =  [a(Y1,Y2)/3(x1,x2)].[a(x1,x2)/a(x1,Y2)] 
(5Y1/ax1)Y2  =  [(aY1/ax1)X2(aY2/ax2)xi  -  (aY1/ax2)Xi2]  /  (3Y2/ax2)x1 
(3Y m,)Y2  =  (8Y WÚY2  -  (dYm2)x:2l(dY2ldX2)x, 


La  derivada  que  no  tiene  signo  definido  aparece  al  cuadrado  y  todas  las  demás  son  positivas 
de  modo  que  se  comprueba  la  desigualdad  mencionada. 

Podemos  afirmar  entonces  que  al  producir  mediante  una  acción  externa  un  pequeño  cambio 
en  X-!,  sin  alterar  X2,  no  solo  se  modificará  Y^  sino  también  Y2;  y  si  las  condiciones  de  equilibrio 
exigen  la  constancia  de  Y2,  entonces  el  cambio  producido  en  Y-i  será  menor.  Por  ejemplo:  (AT)V 
>  (AT)P  para  un  mismo  cambio  en  S. 


Forma  de  los  potenciales 


La  positividad  de  uss  implica  que  u  debe  ser  una  función  convexa  en  s.  Quedándonos  con  la 
rama  de  pendiente  positiva  (temperatura  T  =  us  >  0),  u  resulta  entonces  una  función  creciente 
de  s,  que  exhibe  su  valor  más  bajo  cuando  s=0  y  T=0.  A  medida  que  s  crece,  crece  también  T  y 
así  u  resulta  una  función  monótonamente  creciente  de  T,  a  V  constante: 
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(dU)v,N  =  TdS  =  NCvdT 

Algo  similar  puede  decirse  de  la  entalpia  H,  a  P  constante: 

(dH)pN  =  TdS  =  NCpdT 

La  convexidad  de  u  como  función  de  s,  asegura  que  cualquier  subdivisión  del  sistema  en 
porciones  con  diferentes  s  (s'  y  s”)  y  por  lo  tanto  diferentes  T,  resulta  en  una  situación  con 
mayor  energía  que  el  sistema  homogéneo,  todo  con  la  misma  s  y  la  misma  T.  Puede 
comprobarse  que  en  ese  caso  (ver  figura  6.1),  el  valor  de  u  total  resultante  cae  sobre  la  recta 
secante  a  la  curva  de  equilibrio,  que  pasa  por  los  puntos  (u',s')  y  (u",s”) 

En  efecto  como  S=S'+S”, 

Ns  =  N's'+N"s”  =  N's'+(N-N')s” 

N'/N  =  (s  -s)/(s  -s  ) 

N'VN  =  (s-s  )/(s  -s  ) 


relaciones  conocidas  como  regla  de  la  palanca  que  determinan  la  fracción  molar  en  cada 
porción.  Entonces  la  energía  molar  será: 

U  =  Nu  =  N'u'+  N''u" 

u  =  u'(s''-s)/(s''-s')  +  u"(s-s')/(s"-s')  =  u'  +  (s-s')(u''-u')/(s''-s') 


que  puede  Interpretarse  geométricamente  como  el  punto  sobre  la  secante  de  abscisa  s.  Que  la 
función  sea  convexa  asegura  que  ese  será  siempre  un  punto  de  mayor  energía  que  el 
correspondiente  al  estado  homogéneo. 
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Del  mismo  modo,  la  positividad  de  los  calores  específicos,  compresibilidades  y  otras 
funciones  respuesta,  determina  la  curvatura  de  los  potenciales  termodlnámlcos  en  sistemas 
estables.  Puede  verse  entonces  que: 

FTt  =  -NCV/T  <0  =>  F  es  cóncava  en  T  (rama  decreciente:  S>0) 

Fw  =  (1/Vkt)  >  0  =>  F  es  convexa  en  V  (rama  decreciente:  P>0) 

GTt  =  -NCp/T  <0  =>  G  es  cóncava  en  T  (rama  decreciente:  S>0) 

GPP  =  -  vkt  <  0  =>  G  es  cóncava  en  P  (rama  creciente:  V>0) 

Qtt  =  -[(5T/5S)m,v]_1  <  0  =>  Q  es  cóncava  en  T  (rama  decreciente:  S>0) 

QMM  =  _[(3p/3N)T  v]"1  <  0  =>  Q  es  cóncava,  en  p  (rama  decreciente:  N>0) 

Sugerencia:  graficar  los  diversos  potenciales  y  sus  derivadas  primeras,  en  función  de  sus 
variables  independientes. 
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CAPÍTULO  7 
Tercera  ley 


En  este  capítulo  completaremos  nuestro  formalismo  analizando  el  tercer  postulado  y  sus 
consecuencias. 


Postulado  III 

Como  dijimos,  se  puede  observar  que  los  cambios  en  la  entropía  en  procesos  Isotérmicos 
se  van  reduciendo  al  disminuir  la  temperatura  del  proceso,  de  modo  que  AST  0  cuando  T 
0  (ver  figura  7.1).  Ese  es  el  contenido  del  postulado  de  Nernst  (1906)  que  equivale  a  suponer 
que  la  entropía  de  un  sistema  dado  tiende  a  un  valor  fijo:  S  ^  S0  cuando  T  — ►  0.  Ello  llevó  a 
Planck  (1911)  a  formular  una  ley  más  general,  que  sólo  pudo  justificarse  con  el  posterior 
desarrollo  de  la  mecánica  cuántica  y  que  Identifica  S0  con  0  para  cualquier  sistema.  Es  decir: 

La  entropía  de  cualquier  sistema  se  anula  en  el  estado  para  el  que  T  =  (¿U/¿S)V  Nl  =  0, 
enunciado  coincldente  con  nuestro  Postulado  III. 
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Esta  Ley  no  tiene  consecuencias  tan  trascendentes  como  las  dos  primeras.  Su  importancia 
deriva  de  que  ella  permite  determinar  valores  absolutos  de  S,  no  solamente  el  cambio  AS,  y  por 
consecuencia  permite  estimar  cambios  en  las  energías  libres  AF,  AG,  AQ  ya  que 
AG  =  -S.AT  +  etcétera. 

Microscópicamente  esta  ley  implica  que  el  número  de  microestados  accesibles  al  sistema, 
W,  se  reduce  al  a  T=0,  ya  que  S  =  k.lnW.  Siendo  U  monótonamente  creciente  con  T,  para  T=0 
el  sistema  y  cada  una  de  sus  partes  se  encontrará  en  el  estado  de  menor  energía  posible 
(estado  fundamental).  A  T=0  cualquier  parte  del  sistema  se  encuentra  en  un  determinado 
estado  cuántico  discreto,  el  postulado  requiere  que  este  estado  fundamental  sea  único,  no 
degenerado.  En  general,  si  persiste  una  degeneración  a  T=0,  es  entre  estados  inaccesibles 
entre  sí  (separados  por  grandes  cambios  macroscópicos).  Concluyendo,  a  T=0  el  sistema  cae 
en  su  estado  cuántico  de  menor  energía  y,  en  ese  sentido,  se  vuelve  completamente  ordenado. 


Calores  específicos  y  otras  derivadas  a  bajas  temperaturas 

Veamos  cuales  son  las  consecuencias  de  esta  Ley  sobre  el  comportamiento  de  algunas 
funciones  respuesta  en  las  cercanías  del  cero  absoluto. 

i)  kt  =  ks  y  CP  =  Cv  cuando  T=0. 

El  postulado  identifica  la  isoterma  T=0  con  la  adiabática  S=0,  por  lo  tanto  deben  coincidir  en 
ese  punto  las  propiedades  isotérmicas  con  las  adiabáticas.  En  particular  kt  =  ks. 

Por  otro  lado,  con  las  técnicas  de  reducción  de  derivadas  del  capítulo  5  pudimos  demostrar  que 
para  cualquier  sistema  y  =  CP/CV  =  kt/ks,  e  igual  a  1  en  este  caso.  Es  decir  CP  =  Cv. 
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¡i)  Cp  Cv  — >0,  cuando  T— >0. 


Todos  los  calores  específicos  se  anulan  a  T=0.  En  efecto,  podemos  escribir: 

CA=  (T/N).(0S/5T)a=  (1/N).(0S/5lnT)A. 

Cuando  T->0,  InT— >- °°  y  dado  que  S  tiende  a  un  límite  constante  (0),  esta  derivada  tiende  a 
cero.  Ello  quiere  decir  que  S  — >0  de  un  modo  suave.  Podemos  comprobar  además  que 
(Cp-Cv)  — >0  más  rápido  que  CP  Cv  — >0. 

¡II)  Todas  las  derivadas  Isotérmicas  de  S,  se  anulan  a  T=0. 

Como  ASt  0  cuando  T  ^  0,  entonces: 
dST  =  X(c,S/c’Xk)T.dXk  ->  o  y  entonces 
( SS/cXk)T  ->  0  para  cualquier  variable  extensiva  Xk 
Análogamente 

(3S/cYk)T  -»  0  para  cualquier  variable  intensiva  Yk 

iv)  Todas  las  derivadas  respecto  de  T  se  anulan  a  T=0. 

En  efecto,  a  partir  de  la  propiedad  anterior  y  según  las  relaciones  de  Maxwell: 

(3Xk/3T)Yk  =  (3S/cYk)T  ->  0  cuando  T-»  0  para  cualquier  variable  extensiva  Xk. 

(cYk/3T)Xk  =  -(oS/oXk)T  — >  0  cuando  T->  0  para  cualquier  variable  Intensiva  Yk. 

Entre  ellas,  por  ejemplo,  el  coeficiente  de  dilatación: 
a  =  (1/V).(5V/5T)P  =  -(1/V).(0S/5P)t  -»  0  cuando  T->  0 


Inaccesibilidad  del  cero  absoluto 

El  tercer  postulado  tiene  como  consecuencia  inmediata  la  Imposibilidad  de  reducir  la 
temperatura  de  un  sistema  hasta  el  cero  absoluto  mediante  un  número  finito  de  operaciones  de 
cualquier  tipo  (Nernst,  1912).  En  otras  palabra  el  estado  T  =  0K  no  puede  alcanzarse  por 
ningún  proceso  físico  realizable. 

Esta  Imposibilidad  es  más  bien  de  índole  práctica  que  conceptual.  En  el  laboratorio  se  han 
alcanzado  temperaturas  inferiores  a  InK  (10'9K)  y  se  espera  poder  alcanzar  temperaturas 
menores  aún.  Además,  los  comportamientos  observados  y  predlchos  de  las  variables 
termodinámicas  al  acercarse  a  T=0  son  suaves,  sin  singularidades  ni  divergencias,  de  modo 
que  se  consideran  válidas  las  extrapolaciones  a  0K  de  las  magnitudes  observadas  a 
temperaturas  finitas. 

El  postulado  Identifica  la  Isoterma  T=0  con  la  adiabática  S=0,  aunque  toda  otra  adiabática  e 
isoterma  son  diferentes.  Como  dos  adiabáticas  no  pueden  intersecarse,  no  existe  ninguna 
adiabática,  diferente  de  la  S=0,  que  pueda  intersecar  a  la  Isoterma  T=0.  Entonces  ningún 
proceso  adiabático,  iniciado  a  TfO  podrá  conducir  a  T=0,  es  decir  no  existe  un  proceso 
adiabático  simple  que  conduzca  de  una  temperatura  finita  a  la  temperatura  cero. 
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En  general,  para  enfriar  un  sistema  que  se  encuentra  a  una  temperatura  finita,  primero  se 
reduce  isotérmicamente  su  entropía  (AST)  haciendo  variar  algún  parámetro  X:  XA  ^  XB  y  luego 
se  baja  la  temperatura  isentrópicamente  (proceso  adiabático  cuasiestático:  (ATS)  recuperando 
el  valor  inicial  de  X:  XB  ^  XA  y  así  sucesivamente  como  muestra  la  figura  7.2.  Al  respecto  es 
útil  verificar  que  el  cambio  en  X  al  reducir  la  temperatura  tiene  signo  opuesto  al  producido  al 
reducir  la  entropía.  En  efecto,  para  cualquier  variable  extensiva  X  o  intensiva  Y: 

(3T/5X)S  =  -(5T/5S)x.(aS/aX)T 

y  como  (ST/oS)x  >  0  por  estabilidad, 

las  derivadas  (oT/oX)s  y  (oS/oX)T  tienen  signos  opuestos. 

Podría  tratarse,  por  ejemplo,  de  una  sucesión  de  compresiones  isotérmicas  y  expansiones 
adiabáticas  en  cuyo  caso  X=V.  Por  otra  parte,  según  el  postulado  de  Nernst  AST  ^  0  de  modo 
que  ambas  curvas  S(T)  tienden  a  juntarse  a  T=0  y  pese  a  los  sucesivos  enfriamientos  ATS 
obtenidos,  el  sistema  siempre  encuentra  antes  la  curva  X=XA  que  el  eje  de  ordenadas  (T=0). 


Comportamientos  típicos 


De  ordinario  S^-0,  según  una  ley  de  potencias  ATn.  Es  fácil  comprobar  que  CP,  Cv  y 
derivadas  del  tipo  (SV/5T)P  y  (<3P/c>T)v  tienden  a  cero  según  la  misma  ley. 

Experimentalmente  se  encuentran,  por  ejemplo,  los  siguientes  comportamientos: 

-  sólidos  no  conductores:  C  =  AT3  (vibraciones  de  red  o  fonones) 

-  sólidos  metálicos:  C  =  AT3+BT  (fonones  y  electrones  de  conducción). 

A  los  fines  prácticos,  si  se  conoce  el  calor  específico  como  función  de  la  temperatura:  CP(T), 
pueden  hallarse  por  integración: 
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S(T)  =  f0T  dT  CP(T)/T 
H(T)  =  H0  +  fgT  dT  CP(T) 

G(T)  =  H  -  TS  =  H 0  +  fgT  dT  CP(T)  -  T.Jor  dT  CP(T)/T 
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CAPÍTULO  8 

Inestabilidad  y  transiciones  de  fase 


En  este  capítulo  analizaremos  la  posibilidad  de  Incluir  en  nuestro  formalismo  las  transiciones 
de  fase  y  toda  su  fenomenología  asociada.  Veremos  que  ellas  pueden  asociarse  con  la 
presencia  de  inestabilidades  en  la  relación  fundamental  del  sistema.  Veremos  también  la 
posible  existencia  de  estados  heterogéneos  con  más  de  una  fase  coexistiendo  en  equilibrio. 


Inestabilidad  y  transiciones  de  fase 


La  experiencia  nos  Indica  que  en  algunas  situaciones  el  estado  de  equilibrio  de  un  sistema 
es  heterogéneo  con  dos  o  más  fases  coexistiendo  (líquido  y  vapor  de  una  sustancia  pura,  por 
ejemplo).  También  nos  muestra  cómo  pequeños  cambios  en  las  condiciones  del  sistema 
pueden  producir  cambios  notables  en  sus  propiedades:  decimos  que  el  sistema  sufre  una 
transición  de  fase. 

A  su  vez,  nuestros  principios  afirman  la  posibilidad  de  describir  los  estados  de  equilibrio 
estable  de  un  sistema  simple  con  ciertas  variables  características  (S,V,N),  pero  no  aseguran  la 
existencia  de  un  estado  de  equilibrio  (homogéneo)  para  cualquier  valor  arbitrarlo  de  estas 
variables.  Veremos  cómo  asociar  la  presencia  de  inestabilidades  en  la  ecuación  fundamental 
(violación  de  los  criterios  de  estabilidad)  con  la  heterogeneidad  del  estado  de  equilibrio  y  la 
ocurrencia  de  una  transición  de  fase. 

Supongamos  que  la  convexidad  de  u  como  función  de  s,  requerida  por  la  estabilidad  tal 
como  discutimos  en  el  capítulo  6,  no  se  satisface  en  todo  el  rango  de  las  variables 
termodinámicas,  sino  que  existe  un  intervalo  donde  u  es  cóncava.  Habrá  entonces  una  reglón 
donde  la  subdivisión  del  sistema  en  porciones  con  diferentes  s  (s'  y  s”)  reduce  el  valor  de  su 
energía  interna.  En  esas  condiciones  al  sistema  le  conviene  separarse  en  dos  subsistemas 
homogéneos  (fases)  coexistiendo  en  equilibrio.  Para  ello  es  necesario  que  ambas  tengan  la 
misma  temperatura  y  por  lo  tanto  la  separación  será  en  dos  fases  con  valores  de  s'  y  s”  tales 
que  la  pendiente  a  la  curva  u(s)  sea  la  misma  en  ambos  puntos:  T(s')=  T(s")=T0.  Analizando  la 
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curva  de  la  Figura  8.1  podemos  ver  que  esos  puntos  quedarán  determinados  por  la  envolvente 
inferior  de  la  curva  u(s). 

La  separación  de  fases  se  producirá,  no  sólo  en  la  reglón  en  que  uss^  0,  sino  en  todo  el 
intervalo  s'<s<s”,  en  tanto  que  para  valores  de  s<s'  o  s>s”  el  sistema  tendrá  estados  de 
equilibrio  homogéneos  normales.  SI,  partiendo  de  un  estado  homogéneo  (s<s'),  vamos 
variando  la  temperatura  y  consecuentemente  s,  al  pasar  por  el  valor  de  coexistencia  T0,  el 


Figura  8.1 

sistema  pasará  de  un  estado  homogéneo  (s=s')  a  otro  también  homogéneo  pero  muy  diferente 
(s=s”).  Esta  transición  de  fase  Irá  acompañada  entonces  por  un  salto  discontinuo  en  la 
entropía  molar:  As  =  s''-s'.  Similarmente  habrá  saltos  en  otras  magnitudes  como  el  volumen 
molar  Av  =  v"-v',  tal  como  se  desprende  de  un  análisis  de  la  dependencia  de  u  con  v. 
Dependiendo  entonces  del  valor  de  s  (y  consecuentemente  de  T)  los  estados  de  equilibrio 
serán  homogéneos  o  heterogéneos,  como  se  Indica  en  la  Figura  8.2. 


Figura  8.2 


53 


Dado  que  estas  transformaciones  se  manifiestan  en  saltos  discontinuos  en  las  magnitudes 
S  o  V,  que  son  derivadas  primeras  de  la  energía  libre  de  Gibbs,  se  denominan  transiciones  de 
fase  discontinuas  o  de  primer  orden.  Si  se  repite  el  análisis  anterior  para  distintos  valores  de  P 
y  se  llevan  sobre  un  gráfico  (T,s)  los  valores  de  la  temperatura  de  coexistencia  y  de  las 
entropías  molares  de  las  fases  que  coexisten,  resulta  un  diagrama  de  fases  como  el  mostrado 
en  la  Figura  8.3.  Similares  conclusiones  se  obtienen  para  las  variables  (P,v)  a  distintas  T. 


Figura  8.3 


Notemos  que,  bajo  ciertas  condiciones,  los  estados  homogéneos  pueden  prolongarse  más 
allá  de  estos  límites  como  estados  metaestables.  Ellos  podrían  ocurrir  para  valores  de  s 
mayores  que  s'  o  menores  que  s”,  pero  fuera  del  intervalo  central  de  inestabilidad  estricta 
(uSs<0)  determinado  por  la  condición  uss=0  e  indicado  en  la  figura. 


Coexistencia  de  fases:  regla  de  Gibbs 


El  estado  de  equilibrio  de  un  sistema  simple  de  un  componente,  abierto  al  medio  (T  y  P 
constantes),  queda  determinado  por  el  mínimo  de  la  energía  libre  de  Gibbs  G=G(T,P,N). 
Cuando  este  estado  es  heterogéneo,  las  fases  que  coexisten  en  equilibrio  deben  tener  la 
misma  temperatura  T,  presión  P  y  potencial  químico  p,  ya  que  la  pared  imaginaria  que  las 
separa  no  es  restrictiva  a  ninguna  de  las  variables  termodinámicas  (S,V  o  N). 

Para  la  coexistencia  de  dos  fases  (por  ejemplo  líquido  y  vapor  de  una  sustancia  pura)  deben 
cumplirse  entonces: 


T'=  T"  =  T 
P'=  P"  =  P 
p'(T,P)  =  p”(T,  P). 
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Esta  última  determina  una  dependencia  entre  temperaturas  y  presiones  de  coexistencia 
P=P(T).  Es  la  llamada  curva  de  equilibrio  o  de  coexistencia  (ver  Figura  8.4)  sobre  la  cual  es 
posible  la  coexistencia  de  ambas  fases.  Fuera  de  esa  curva  los  estados  de  equilibrio  serán 
homogéneos  de  un  tipo  u  otro,  al  atravesar  la  curva  ocurrirá  la  transición  de  fase 
(discontinuidad  en  s  y  v)  y  sólo  sobre  la  curva  será  posible  la  coexistencia.  En  el  caso  de  la 
transición  líquido  -  vapor,  corresponde  a  la  curva  de  presión  de  vapor  en  función  de  la 
temperatura.  También  pueden  describirse  de  manera  similar,  curvas  de  coexistencia  sólido- 
líquido  (fusión)  y  sólido-vapor  (sublimación). 


Para  la  coexistencia  de  tres  fases  (por  ejemplo  sólido,  líquido  y  vapor  de  una  sustancia 
pura)  deben  cumplirse 
T'  =  T"  =  T"'  =  T 

p-  =  p-  =  p-"  =  P 

p'(T,P)  =  p”(T,  P)  =  m"'(T,P). 

Estas  últimas  determinan  una  doble  dependencia  entre  temperaturas  y  presiones  de 
coexistencia  que  una  vez  resueltas  conducen  a  un  único  valor  de  T  y  P,  punto  triple  (Tt,Pt) 
sobre  el  cual  es  posible  la  triple  coexistencia.  Como  se  muestra  en  la  Figura  8.5,  este  punto  de 
coexistencia  de  tres  fases  es  también  el  punto  donde  se  unen  las  tres  curvas  de  coexistencia 
de  dos  fases  (sólido-líquido,  sólido-vapor  y  líquido-vapor,  por  ejemplo). 

Con  un  razonamiento  análogo  podemos  ver  que  en  tal  sistema  sería  imposible  la 
coexistencia  de  4  fases  en  equilibrio  ya  que  resultarían  tres  ecuaciones  para  T  y  P,  quedando 
sobredeterminada  la  solución. 
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El  resultado  puede  generalizarse  al  caso  de  sistemas  de  r  componentes,  caracterizados  por 

G=G(T,P,  N-i . Nr),  coexistiendo  en  cp  fases.  En  cada  fase  tendremos  (r-1)  variables  referidas  a 

la  composición,  además  de  T  y  P  comunes  a  todas  las  fases,  totalizando  (r-1)cp+2  variables  a 
determinar.  Las  condiciones  de  coexistencia  (condiciones  de  equilibrio  entre  las  fases)  implican 
la  Igualdad  de  T,P  y  todos  los  potenciales  químicos  p¡  (1=1,..  r)  sobre  todas  las  fases.  Tomando 
a  éstos  como  funciones  de  T  y  P,  de  su  Igualdad  resultan  r(cp-1 )  ecuaciones.  Para  que  el 
sistema  de  ecuaciones  no  resulte  sobre-determinado  debe  cumplirse:  (r-1  )cp+2  >  r(cp-1 )  de 
donde  surge  la  regla  de  las  fases  de  Gibbs\ 

cp  <  (r+2). 

El  número  de  fases  que  pueden  coexistir  en  equilibrio  es  como  máximo  r+2.  Esta  regla 
suele  expresarse  en  términos  del  número  de  grados  de  libertad  v  (número  de  variables 
intensivas  que  pueden  variarse  libremente)  que  resulta: 

v  =  (r+2)-cp  >  0. 

Aplicada  al  caso  simple  discutido  al  comienzo  (r=1)  Implica  que: 

-  para  1  fase,  quedan  2  grados  de  libertad:  T  y  P  pueden  variarse  independientemente; 

-  para  2  fases,  queda  1  grado  de  libertad  (T,  que  determina  P  sobre  la  curva  de  coexistencia); 

-  para  3  fases,  no  quedan  grados  de  libertad  (T  y  P  están  fijos  en  el  punto  triple). 


Regla  de  la  palanca 


El  estado  de  equilibrio  de  un  sistema  abierto  al  medio  (T  y  P  constantes),  queda 
determinado  por  el  mínimo  de  la  energía  libre  de  Glbbs.  Cuando  este  estado  es  heterogéneo, 
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las  fases  que  coexisten  en  equilibrio,  además  de  tener  los  mismos  valores  de  T,  P  y  p,  deben 
tener  la  misma  energía  libre  de  Gibbs  molar  g'=  g"=  g.  Caso  contrario,  como  G  =  N'g'+N"g”, 
el  sistema  preferiría  estar  por  completo  (estado  homogéneo)  en  la  fase  de  menor  energía  libre 
de  Gibbs  molar. 

En  las  condiciones  de  coexistencia  G  =  (N'+N")g,  y  en  principio  quedan  indeterminadas  las 
fracciones  de  sistema  en  cada  fase:  N'  y  N'\  El  equilibrio  es  Indiferente  a  dichas  cantidades, 
pero  ellas  sí  determinan  el  valor  de  magnitudes  como  la  entropía  y  el  volumen  totales: 

S  =  N's'+N"s”,  V  =  N'v'+N"v". 

SI  se  controla  externamente  alguna  de  estas  variables  (V  o  S),  las  fracciones  molares  x'=  N'/N 
y  x"=N"/N  =  1-x',  pueden  determinarse.  Por  ejemplo: 

V  =  N'v'+N"v"  — >  v  =  x'v'+x"v"  =  x'(v'-v”)  +  v" 

resultando  x'  =  (v''-v)/(v''-v') 

y  análogamente:  x"  =  (v-v')/(v"-v'), 

expresiones  que  se  conocen  como  regla  de  la  palanca. 

Sugerencia:  analizar  las  evidencias  experimentales  sobre  la  condensación  de  un  fluido  y  su 
descripción  a  partir  de  la  ecuación  de  estado  de  van  der  Waals:  P  =  RT/(v-b)  -  a/v2. 


Ecuación  de  Clausius  -  Clapeyron 


Hemos  visto  cómo  una  transición  de  fase  Implica  discontinuidades  As  y  Av  entre  las  fases  y 
también  cómo  estas  fases  pueden  coexistir  sobre  la  curva  de  coexistencia  P(T).  Veremos  cómo 
puede  relacionarse  la  pendiente  de  esta  curva  con  dichas  dlscontinidades  a  través  de  la 
ecuación  de  Clausius  -  Clapeyron. 

En  efecto,  sobre  la  curva  de  coexistencia: 

9'  =  9" 
dg'  =  dg” 

utilizando  la  ecuación  de  Gibbs  molar,  dg=-s.dT+v.dP,  para  cada  fase: 


-s'dT+v'dP  =  -s''dT+v''dP 
(s"-s')dT  =  (v"-v')dP 
y  finalmente: 

(dP/dT)cc  =  As/Av 

que  es  la  ecuación  buscada. 

El  cambio  en  la  entropía  se  vincula  al  calor  latente  de  transición  (vaporización  o 
condensación,  fusión  o  solidificación,  etc):  As  =  L/T.  Además,  para  la  transformación  líquido 
vapor  Av  =  vgas  -  V|¡q  *  vgas  =  RT/P.  Así  la  ecuación  toma  la  forma  más  conocida: 
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(dP/dT)cc  =  LP/RT2 
o  bien:  dP/P  =  LdT/RT2 

que  puede  integrarse  para  dar :  InP  =  cte  -  L/RT 

En  sistemas  complejos,  donde  se  requieren  otras  variables  extensivas  X¡  para  la  descripción 
de  los  estados  de  equilibrio,  con  sus  variables  intensivas  conjugadas  Y¡,  cuya  igualdad 
determina  el  equilibrio  entre  fases  con  diferentes  X¡,  la  pendiente  de  la  curva  de  coexistencia  en 
el  plano  (Y¡,T)  viene  dada  por  la  ecuación  de  Clausius  Clapeyron  generalizada: 


(dYj/dT)cc  =  -  AS/AX¡ 


Sugerencia:  discutir  los  diagramas  de  fase  (P,T),  (P,v)  y  (T,s)  de  diversas  sustancias. 


Punto  crítico.  Transiciones  continuas 

En  algunos  casos  la  diferencia  entre  fases  (discontinuidades)  va  disminuyendo  al  aumentar 
la  temperatura,  haciéndose  las  fases  separadas  por  la  transición  cada  vez  mas  similares  entre 
sí.  Cuando  esta  diferencia  desaparece  (As  =  Av  =  0)  el  sistema  se  encuentra  en  un  punto 
crítico  (TC,PC)  por  encima  del  cual  ya  no  habrá  fases  diferentes  ni  transición  de  fase. 

Tenemos  entonces  que,  como  muestra  la  figura  8.3, 
para  T>TC  los  estados  de  equilibrio  son  homogéneos  (1  fase  única); 
para  T<TC  los  estados  de  equilibrio  son  heterogéneos  (2  fases  coexistiendo); 
para  T=TC  una  única  fase  (las  dos  fases  de  baja  temperatura  coalescen  en  una  única  al 
desaparecer  sus  diferencias). 

En  esos  casos  la  curva  de  coexistencia  termina  en  el  punto  crítico  (figuras  8.4  y  8.5).  Un 
ejemplo  típico  lo  constituye  el  punto  crítico  en  la  transición  líquido  gas  de  un  fluido.  La 
transición  de  fase  que  ocurre  cuando  se  pasa  por  el  punto  crítico  se  denomina  continua  o  de 
segundo  orden  y  será  analizada  en  los  capítulos  que  siguen. 

Además  del  punto  crítico  de  un  fluido  podemos  mencionar  como  ejemplos  el  punto  de 
Curie  de  un  ferromagneto  (Tc,  H=0),  el  punto  crítico  de  la  transición  superconductora  a 
campo  aplicado  nulo  y  los  puntos  críticos  asociados  a  composición  y  grado  de  orden  en 
una  mezcla  binaria. 
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CAPITULO  9 

Sistemas  magnéticos 


En  este  capítulo  Introduciremos  las  interacciones  electromagnéticas  en  nuestro  formalismo  y 
describiremos  la  fenomenología  de  algunos  sistemas  magnéticos,  en  particular  la  transición  de 
para-  a  ferromagneto  que  es  una  transición  de  fase  de  primer  orden  comparable  a  la  de 
condensación  de  un  fluido. 


Termodinámica  de  un  sistema  magnético 


Para  describir  las  propiedades  termodinámicas  asociadas  con  la  presencia  de  un  campo 
magnético  e  Identificar  las  variables  relevantes,  debemos  analizar  el  término  de  trabajo 
magnético  asociado  a  la  entrega  de  energía  al  sistema.  Supongamos  que  nuestro  sistema 
consiste  de  una  muestra  de  volumen  V  encerrada  por  una  superficie  A  en  presencia  de  campos 
electromagnéticos  E  y  H.  De  acuerdo  con  el  teorema  de  Poynting  la  energía  entrante  al  sistema 
en  dt,  sera: 

dWem  =  -(c/4tt)  jdA.(EAH)  dt 
dWem  =  -(c/4tt)  jdVdlv(EAH)  dt 
dWem  =  (c/4tt)  JdV[E.rotH  -  H.rotE])  dt 


Teniendo  en  cuenta  las  ecuaciones  de  Maxwell 
rotE  =  -(1/c)3B/3t  rotH  =  (1/c)3D/3t+(4TT/c).jl¡5re 

y  suponiendo  que  no  hay  corrientes  libres  dentro  del  volumen  V 

dWem  =  (1/4tt)  jdV[E.,9D/<9t  +  H.3B/3t])  dt 
dWem  =  (1/4tt)  jdV[E.dD  +  H.dB] 
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donde  reconocemos  claramente  los  términos  eléctrico  y  magnético : 


dWe  =  (1/4tt)  jdV(E.dD) 
dWm  =  (1/4tt)  jdV(H.dB) 


Suponiendo  que  en  nuestro  sistema  el  vector  desplazamiento  eléctrico  D  es  uniforme  sobre 
el  volumen  V,  sólo  tendremos  el  término  de  trabajo  magnético.  SI  además  asumimos  que  el 
campo  magnético  y  el  material  son  uniformes  podemos  integrar  sobre  el  volumen  y  obtener: 

dWm  =  (V/4-rr)  H.dB  =  Hd6 
con  S  =  (V/4tt)Bh 

En  caso  contrario,  el  volumen  dejará  de  ser  una  buena  variable  y  deberemos  trabajar  con 
magnitudes  locales  (por  unidad  de  volumen). 

Podemos  adoptar  6  como  variable  extensiva  y  H  como  su  intensiva  conjugada  y  elaborar  la 
termodinámica  de  un  sistema  magnético  a  partir  de  ellas.  Sin  embargo,  es  conveniente  separar 
del  término  de  trabajo,  el  correspondiente  al  establecimiento  del  propio  campo,  presente  aún  en 
ausencia  del  material  que  constituye  nuestro  sistema.  En  el  sistema  de  unidades  gausseano4 
que  estamos  usando:  B  =  H+4ttM  y  entonces 

dWm  =  (1/4tt)  jdV(H.dH)  +  JdV(H.dM) 
dWm  =  (V/8tt)  d(H2)  +  V(H.dM) 

Dejando  de  lado  el  trabajo  asociado  al  campo  H,  la  nueva  expresión  para  el  trabajo 
magnético  será: 

dWm  =  VH.dM  =  HdM 

pudlendo  Identificar  las  variables  termodinámicas  extensiva  M  =  VMH  e  intensiva  H.  De  manera 
similar  habríamos  obtenido  para  el  caso  eléctrico,  donde  D  =  E+4ttP: 

dWe  =  VE.dP  =  EdP 

con  variables  extensiva  P  =  VPE  e  Intensiva  E. 

Volviendo  al  caso  que  nos  ocupa,  la  nueva  expresión  para  la  ecuación  de  Glbbs  resulta: 
dU  =  TdS  -  PdV  +  HdM  +  Z¡  |i¡dN¡. 
correspondiente  a  una  ecuación  fundamental: 

4  B  =  |j0(H+M)  en  el  internacional  SI 
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U  =  U(S,V,M,  Ni,  Nr). 


Suponiendo  N=cte  y  despreciando  los  cambios  de  volumen  que  en  general  son  pequeños  en 
sólidos,  resulta 

U  =  U(S,M)  dll  =  TdS  +  HdM 

pudlendo  adecuar  al  caso  presente  todas  la  expresiones  vistas  hasta  aquí,  mediante  el  simple 
cambio  de  variables: 

(-P,V)-»(H,M) 

De  manera  similar  a  lo  hecho  para  sistemas  simples,  podemos  definir  nuevos  potenciales 
termodlnámlcos  cambiando  variables  mediante  transformaciones  de  Legendre.  Por  ejemplo: 

F  =  F(T,M)  =  U-TS  dF  =  -SdT  +  HdM 

con  S  =  -(3F/3T)m  H  =  (3F/3M)T 

y  la  relación  de  Maxwell:  -(3S/3M)t  =  (3H/3T)m. 

O  bien,  el  análogo  de  la  energía  libre  de  Glbbs5 
G  =  G(T,H)  =  U  -  TS  -  HM  dG  =  -SdT  -  MdH 

con  S  =  -(3G/3T)h,  M  =  -(3G/3H)T 

y  la  relación  de  Maxwell:  (3S/3H)T  =  (3M/3T)H 

Podemos  definir  también  algunas  funciones  respuesta  como  la  susceptibilidad  magnética: 

Xt  =  (3M/3H)t  =  -  Ghh  >  0 

y  los  capacidades  caloríficas: 

CH  =  T(3S/3T)h  =  -TGtt  >  0 
CM  -  T(3S/3T)m  =  -TFtt  >  0 

Con  los  métodos  habituales  de  reducción  de  derivadas  podemos  probar  la  relación: 
Ch-Cm=T[(3M/3T)h]2/Xt 

semejante  a  la  ya  demostrada  para  CP  -  Cv  en  el  Capítulo  5. 


5  En  rigor  G=G(T,P)  y  cuando  se  preste  a  confusión  deberíamos  usar  otra  notación  para  G(T,H) 
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Propiedades  magnéticas 


El  comportamiento  magnético  más  sencillo  es  el  paramagnetismo  que  exhiben  muchas 
sustancias  sólidas  las  que,  en  presencia  de  un  campo  magnético  H,  adquieren  una 
magnetización  M  en  la  misma  dirección.  Esta  es  proporcional  a  H  para  bajos  campos  y  tiende  a 
un  valor  de  saturación  para  campos  grandes.  La  constante  de  proporcionalidad  depende  de  la 
temperatura  según  la  ley  de  Curie  y  entonces:  M  =  CH/T  donde  C  es  la  llamada  constante  de 
Curie.  La  correspondiente  susceptibilidad  será:  Xt  =  C/T.  De  acuerdo  con  la  citada  relación  de 
Maxwell 

(3S/<9H)t  =  (<9M/<9T)h  =  -CH/T2  <  0 


indicando  una  reducción  de  la  entropía  al  aplicar  un  campo  en  condiciones  Isotérmicas. 


Figura  9.1 

Microscópicamente  el  comportamiento  paramagnético  de  un  sólido  puede  explicarse  a  partir 
de  la  existencia  de  momentos  magnéticos  individuales  p,  asociados  a  los  átomos  de  alguno  de 
los  constituyentes,  fijos  en  puntos  dados  de  la  red  cristalina.  En  ausencia  de  campo,  estos 
momentos  se  orientan  al  azar.  Un  campo  magnético  externo  tiende  a  alinearlos  en  su  misma 
dirección  compitiendo  con  la  agitación  térmica  y  alcanzando  la  magnetización  de  saturación 
cuando  logra  orientar  la  totalidad  de  los  momentos.  En  ese  caso,  si  hay  N  momentos 
individuales,  el  momento  magnético  total  será  Ms  =  Np  La  ley  de  Curie  pone  de  manifiesto  la 
mayor  dificultad  para  orientar  los  momentos  a  medida  que  aumenta  la  temperatura.  Así,  las 


62 


isotermas  H(M)  típicas  de  un  paramagneto  tendrán  la  forma  mostrada  en  la  Figura  9.1,  donde 
se  puede  apreciar  la  evolución  de  esta  dependencia  al  ir  variando  la  temperatura. 

Algunas  sustancias,  como  el  Fe  y  otros  metales  de  transición  y  sus  óxidos  presentan,  por 
debajo  de  una  temperatura  característica  Tc,  una  magnetización  no  nula  aún  en  ausencia  de 
campo,  fenómeno  conocido  como  ferromagnetismo.  A  altas  temperaturas  T>TC  las  Isotermas 
H(M)  son  las  típicas  de  un  paramagneto,  con  pendientes  iniciales  cada  vez  más  pequeñas  al 
disminuir  la  temperatura.  Esta  pendiente  inicial  se  anula  al  alcanzar  T=TC  donde  la  Isoterma 
correspondiente  ( isoterma  crítica)  presenta  un  punto  de  inflexión: 

(3H/aM)T  =  ( 32H/<9M  2)t  =  0  para  H=0 

Para  temperaturas  inferiores  T<TC,  aún  a  campo  nulo  aparece  una  magnetización  finita  que 
luego  crece  con  el  campo  hasta  alcanzar  el  valor  de  saturación.  Esta  magnetización 
espontánea  y  el  momento  magnético  asociado  Me  aumentan  al  disminuir  la  temperatura, 
tendiendo  éste  a  Ms  cuando  T->0. 

Estas  características  se  muestran  en  los  diagramas  de  fase  H  vs  T  y  M  vs  T  de  las  figuras 
9.2  y  9.3  respectivamente. 


Figura  9.2 
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Figura  9.3 


Las  propiedades  del  punto  crítico  (Tc,  H=0)  serán  analizadas  en  los  capítulos  siguientes  al 
discutir  las  transiciones  de  fase  continuas. 


Modelo  de  Curie-Weiss 

Con  los  métodos  habituales  de  la  mecánica  estadística  pueden  describirse  bien  las 
propiedades  de  un  sistema  de  N  momentos  magnéticos  p  no  interactuantes  en  presencia  de  un 
campo  externo  H  (paramagnetismo).  Cuánticamente  estos  momentos  magnéticos  pueden 
asociarse  al  momento  angular  total  J  de  los  atomos  involucrados,  según:  p  =  gpB<J>,  donde  g 

es  el  factor  giromagnético  y  pB  =  eh/4TTmc  es  el  magnetón  de  Bohr.  Así  resulta  una  ecuación  de 

estado  magnética: 

M(T,H,N)  =  NgpBJ.Bj(y) 

donde  y  =  gpBJH/kT  y  aparece  la  función  de  Brillouin6  Bj. 

Para  el  caso  particular  de  electrones  J=1/2  y  g=2  y  entonces  p=pB,  resultando 

M(T,H,N)  =  NpB.tgh(pBH/kT) 


6  Bj(y)  =  (1+1/2J).cotgh[(1  +  1/2J)y]  -  (1/2J).cotgh[(1/2J)y] 
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Esta  ecuación  de  estado  describe  adecuadamente  el  comportamiento  paramagnético 
observado,  que  es  el  de  un  sistema  ideal  sin  interacciones.  Puede  verse  por  ejemplo  que,  para 
campos  pequeños,  cuando  puede  aproximarse  tangh(y)  =  y, 

M  =  NpB2H/kT 

satisfaciendo  la  ley  de  Curie  con  C=  NpB2/k. 

Por  otro  lado,  para  poder  describir  la  transición  de  fase  es  necesario  incluir  las  interacciones 
entre  momentos.  Con  ese  fin,  Weiss  desarrolló  una  teoría  de  campo  medio  en  la  que  el  efecto 
de  estas  interacciones  se  refleja  en  la  existencia  de  un  campo  adicional  Hcm,  proporcional  a  la 
magnetización,  que  actúa  sobre  cada  momento  individual,  además  del  campo  externo.  Así  el 
campo  efectivo  será: 

Hef  =  H  +  Hcm  =  H  +  XM 
con  lo  que  resulta: 

M(T,H,N)  =  NpB.tanh[pB(H+XM)/kT] 

que  es  una  ecuación  trascendente  para  M. 

Vamos  a  examinar  si  este  modelo  es  capaz  de  predecir  una  magnetización  espontánea,  es 
decir  M  ^0  para  H=0.  En  ese  caso: 

M  =  NpB.tanh(pB^M/kT) 

y  la  ecuación  trascendente  puede  resolverse  graficando  fi=M  y  f2=  NpB.tanh(pBXM/kT).  La 
condición  para  que  í-¡  y  f2  tengan  una  intersección  no  trivial,  además  de  aquella  para  M=0,  es 
que  la  pendiente  al  origen  de  la  función  f2  sea  mayor  que  1  (ver  figura  9.4).  Teniendo  en  cuenta 
que  para  pequeños  valores  del  argumento  tangh(y)  «  y,  cerca  del  origen:  f2=NpB2XM/kT  y  tal 
pendiente  será  NpB2X/kT  >  1.  Esta  condición  define  entonces  una  temperatura  crítica 

Tc=  NpB2Wk  =  CX 

por  debajo  de  la  cual  aparecerán  soluciones  no  triviales  indicativas  del  ordenamiento 
espontáneo. 
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Figura  9.4 

Por  otra  parte,  para  campos  no  nulos  pero  pequeños  y  T>TC  podemos  aproximar: 

M  =  NpB2(H+XM)/kT 
de  donde  podemos  despejar: 

M(1-NpB2WkT)  =  NpB2H/kT 
M  =  NpB2H/k(T-Tc) 

Que  expresa  una  dependencia  similar  a  la  de  la  ley  de  Curie  con  C  =  NpB2/k  y  un 
deplazamlento  del  eje  de  temperaturas  en  Tc,  tal  como  se  observa  en  ferromagnetos. 


Fenómenos  combinados 

Cuando  no  puedan  despreciarse  los  cambios  de  volumen,  tendremos  que  considerar  la 
siguiente  expresión  para  la  ecuación  de  Glbbs,  a  N=cte: 

dll  =  TdS  -  PdV  +  HdM 

Podemos  describir  entonces  algunos  fenómenos  que  combinan  cambios  en  las  variables 
térmicas,  mecánicas  y  magnéticas  y  definir  las  funciones  respuesta  asociadas. 
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Efecto  magnetocalórico 


cambio  en  la  temperatura  de  un  magneto  por  un  cambio  en  el  campo  magnético  aplicado, 
en  condiciones  adiabáticas: 

dT  =  (<9T/<9H)S,P  dH. 

Esta  función  respuesta  puede  expresarse  en  términos  de  otras  como  sigue: 

(<9TA9H)s,p  =  -(<9TA9S)h,p.(<9S/<9H)t,p  =  -(T/Ch.p).(3M/3T)h,p  >  0 


Eefecto  magnetostrictivo 

cambio  en  el  volumen  de  un  magneto  por  un  cambio  en  el  campo  magnético  aplicado,  en 
condiciones  isotérmicas  o  adiabáticas: 

dV  =  (dV/dl H)x,p  dH.  X  =  T,S 

Puede  demostrarse  la  siguiente  relación  entre  ellas,  donde  aparecen  el  coeficiente  de 
dilatación  y  el  magnetocalórico: 

(<9V/<9H)s,p  =  (3V/3  H)t,p  +  (<9V/<9T)h,p.(<9TA9H)s,p 

En  ocasiones  se  la  define  a  partir  del  cambio  en  la  longitud  del  imán  en  lugar  de  su  volumen. 


Eefecto  magnetoelástico 

cambio  en  el  momento  dipolar  M  por  un  cambio  en  la  presión,  aplicado  en  condiciones 
isotérmicas  o  adiabáticas: 

dM  =  (<9M/<9P)x,h  dP.  X  =  T,S 

Es  el  simétrico  del  efecto  anterior,  intercambiando  los  roles  de  las  variables  mecánicas  y 
magnéticas  entre  sí.  Puede  demostrarse  la  siguiente  relación: 

(<9M/3P)s,h  =  (dM/dP)T,H  +  (<9M/3T)p,H.(dT/dP)s,H. 

En  ocasiones  se  la  define  a  partir  de  variaciones  de  la  tensión  (unidimensional)  sobre  el  imán. 
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Sugerencia:  analizar  las  posibles  relaciones  de  Maxwell  entre  estos  coeficientes. 
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CAPITULO  10 

Transiciones  de  fase  continuas 


En  este  capítulo  consideraremos  los  fenómenos  críticos,  es  decir,  el  comportamiento  de  los 
sistemas  en  las  vecindades  del  punto  crítico.  Analizando  la  condensación  de  un  fluido  y  la 
transición  fe rromagn ética,  encontraremos  comportamientos  comunes  reflejados  en  la  similitud 
de  las  ecuaciones  resultantes  al  escribirlas  en  variables  reducidas. 


Punto  crítico 

Hemos  visto  cómo  en  la  condensación  de  un  fluido,  al  acercarnos  al  punto  critico  por  la 
curva  de  coexistencia,  las  discontinuidades  As  y  Av,  características  de  una  transición  de  fase 
de  primer  orden,  se  van  reduciendo  hasta  desaparecer  en  el  punto  crítico  (As=Av=0).  Lo  mismo 
ocurre  con  la  discontinuidad  en  el  momento  magnético  (AM  =  2Me),  al  aproximarnos  a  campo 
nulo,  a  la  Tc  de  un  ferromagneto.  Las  fases  separadas  por  la  transición  se  hacen  cada  vez  mas 
similares  entre  sí  hasta  Igualarse  en  el  punto  crítico,  por  encima  del  cual  ya  no  habrá  fases 
diferentes  ni  transición  de  fase.  Tenemos  entonces  que: 

-  a  T>TC  los  estados  de  equilibrio  son  homogéneos  (1  fase  única) 

-  a  T<TC  los  estados  de  equilibrio  son  heterogéneos  (2  fases  coexistiendo)  o  son  múltiples 
(vahas  fases  posibles) 

-  a  T=TC  habrá  una  única  fase  (las  fases  de  baja  temperatura  coalescen  en  una  única  al 
desaparecer  sus  diferencias). 

Podemos  decir  entonces  que  al  bajar  la  temperatura  por  debajo  de  Tc  hay  una  ruptura  de 
simetría  en  el  sistema  y  se  requiere  de  una  variable  adicional  para  describir  sus  estados  de 
equilibrio.  Llamamos  parámetro  de  orden  q  a  esta  variable  y  le  atribuimos  las  siguientes 
propiedades: 
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-  a  T>TC,  n  =  O 

-  a  T<TC,  n  ^  O  y  toma  más  de  un  valor 

-  cuando  T— >TC,  n  — >  0  continuamente. 

En  los  sistemas  magnéticos,  fácilmente  puede  Identificarse  n  =  M.  Para  la  condensación  de  un 
fluido  se  proponen  r]  =  vg-V|  o  r|  =  5g-5|  (5=N/V=v"1),  que  satisfacen  las  propiedades  enunciadas. 
Por  su  parte,  tanto  para  la  condensación  de  un  fluido  como  para  la  transición  ferromagnétlca, 
las  derivadas  segundas  de  la  energía  libre  de  Glbbs,  capacidades  caloríficas  y 
compresibilidades  o  susceptibilidades  son  discontinuas  o  divergentes  en  el  punto  crítico. 

La  transición  de  fase  que  ocurre  cuando  el  sistema  pasa  por  el  punto  crítico  se  denomina 
entonces,  continua  o  de  segundo  orden  y  el  comportamiento  de  las  distintas  variables 
termodinámicas  en  la  vecindad  del  mismo,  se  conoce  como  fenómenos  críticos.  Además  de  los 
casos  mencionados  podemos  citar  como  ejemplos:  el  punto  crítico  composlclonal  y  el  de 
transición  orden-desorden  en  una  mezcla  binarla  y  el  de  transición  superconductora  a  campo 
aplicado  nulo,  que  serán  discutidos  más  adelante. 


Estabilidad  del  punto  crítico 


En  los  dos  casos  considerados  hemos  visto  que  las  Isotermas  críticas  (T=TC)  presentan  un 
punto  de  ensilladura  en  el  punto  crítico,  es  decir: 

(0P/5v)t  =  (02P/3v2)t  =  0,  para  T=TC,  v=vc 
(5H/5M)t  =  (52H/5M2)t  =  0,  para  T=TC,  M=0. 

Recordando  que  la  estabilidad  de  los  estados  de  equilibrio  de  un  fluido  requiere: 

-(oP/ov)T  =  1/vkt  >  0 

debemos  revisar  los  criterios  de  estabilidad  para  hallar  las  nuevas  condiciones  suficientes,  al 
menos  sobre  la  Isoterma  crítica. 

Como  vimos  cuando  discutimos  los  criterios  de  estabilidad  intrínseca,  el  equilibrio  de  un 
sistema  aislado  (energía  total  mínima)  requiere: 

6Ut  =  dUt  +  Nd2u/2!  +  Nd3u/3!  +  Nd4u/4!  +..  >  0 

y  en  condiciones  normales  es  suficiente  que 

d2u  =  uS3(ds)2  +  2usv(ds)(dv)  +  uw(dv)2  >  0,  Vds,  Vdv 

o  bien,  introduciendo  T: 
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d2u  =  (1/uss)(dT)2  +  (uvv-  usv2/uss)(dv)2  >  O  VdT,  Vdv 

Sobre  una  isoterma  dT=0  y  el  coeficiente  del  término  en  (dv)2,  que  puede  escribirse  -(oP/ov)T, 
justamente  se  anula  en  el  punto  crítico.  Yendo  entonces  a  órdenes  superiores,  los  siguientes 
términos  relevantes,  siempre  sobre  una  isoterma,  serán: 

-  (52P/av2)T.(dv)3/3!  -  (d3P/cV3)T.(dv)4/4!  -  ...  >  0  Vdv 

El  primer  término,  impar,  debe  anularse  y  ello  se  satisface  por  ser  el  punto  crítico  un  punto  de 
inflexión: 

(52P/5v2)t  =  0. 

Entonces  el  nuevo  criterio  de  estabilidad  para  el  punto  crítico  será: 

-  (53P/av3)T  >  0. 


Principio  de  estados  correspondientes 

Consideremos  un  fluido  de  van  der  Waals,  que  satisface 
P  =  RT/(v-b)  -  a/v2, 

y  en  el  punto  crítico  debe  cumplir: 

(<3P/3v)t  =  -RTc/(vc-b)2  +  2a/vc3  =  0 
(<32P/Sv2)t  =  2RTc/(vc-b)3  -  6a/vc4  =  0 
(03P/5v3)t  =  -6RTc/(vc-b)4  +  24a/vc5  <  0  . 

De  las  dos  primeras  y  de  la  ecuación  de  estado,  podemos  obtener  las  coordenadas  del  punto 
crítico  en  términos  de  las  constantes  a,  b  y  R: 


vc  =  3b 
Tc  =  8a/27bR 
Pc  =  a/27b2. 

Podemos  verificar  además  que  con  estos  valores  se  satisface  el  criterio  de  estabilidad,  ya  que 
(c>3P/Sv3)t  =  -a/81  b5  <  0  . 
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Introduciendo  las  coordenadas  del  punto  crítico  en  la  ecuación  de  estado,  para  lo  que  se 
sugiere  calcular  directamente  P/Pc,  resulta: 

P/Pc  =  8(T/Tc)/(3v/vc  -  1 )  -  3/(v/vc)2 

o  bien,  definiendo  coordenadas  reducidas : 

P=P/PC;  T=T/TC;  v=v/vc, 

P  =  877(3v-1)  -  3/v^  =  P(T,v) 

que  es  una  relación  universal,  válida  para  todos  los  fluidos  de  van  der  Waals,  cuando  se 
expresan  las  magnitudes  en  unidades  críticas.  Esta  predicción  viene  a  corroborar  la 
observación  empírica,  conocida  como  principio  de  estados  correspondientes  y  reflejada  en  los 
conocidos  gráficos  de  Guggenheim.  En  ellos  puede  apreciarse  cómo  sistemas  muy  diversos 
(Ne,  Ar,  Kr,  Xe,  N2,  02,  CO  y  CH4)  muestran,  en  la  transición  líquido-vapor,  un  único  diagrama 
de  fases  (T vs  v-1=  p/pc)  y  una  única  curva  de  equilibrio  (InP  vs  T1)  para  valores  de  T  entre  Tc  y 
0.55TC. 

Coordenadas  y  ecuaciones  de  estado  reducidas 

Para  analizar  en  detalle  el  comportamiento  de  las  propiedades  termodinámicas  en  la 
vecindad  del  punto  crítico  es  conveniente  definir  magnitudes  reducidas  que  miden  el 
apartamiento  de  la  dada  propiedad  respecto  de  su  valor  crítico,  en  unidades  críticas  (cuando  el 
valor  crítico  no  es  cero).  Así  definimos  la  temperatura  reducida: 

6  =  (T-Tc)/Tc  =  r-i 
y  otras  como: 

tt  =  (P-Pc)/Pc  =  P- 1  u  =  (v-vc)/vc  =  v  -1 . 

Veamos  como  queda  expresada  la  ecuación  de  estado  de  van  der  Waals  en  estas 
coordenadas: 

tt  =  8(1+0)/(2+3u)  -  3/(1 +u)2-  1  =  [80-3u3/(1+u)2]/(2+3u). 

Como  al  aproximarnos  al  punto  crítico  u-»0,  podemos  desarrollar  en  serie: 
tt  «  (40— 3u3/2+3u4-  ...).(1-3u/2+9u2/4-  ...) 
y  conservando  sólo  los  términos  dominantes: 
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tt  =  40  -  60u  -  3u3/2 

que  describe  las  isotermas  de  van  der  Waals  en  las  cercanías  del  punto  crítico.  Se  puede 
verificar  que  satisface  las  propiedades  requeridas. 


Exponentes  críticos 


Con  el  fin  de  cuantiflcar  el  comportamiento  de  las  diversas  magnitudes  termodinámicas  al 
aproximarse  al  punto  crítico,  es  conveniente  definir  los  exponentes  críticos-,  exponentes 
dominantes  en  el  comportamiento  aslmptótlco  al  tender  al  punto  crítico.  Por  ejemplo,  si  la 
propiedad  F(0)  se  comporta  como  F(0)  «  ¡01a,  cuando  0— >0,  diremos  que  tiene  un  exponente 
crítico  X.  En  general,  definimos: 

X  =  lim0_o[ln|F(0)|/ln|0|]. 

Valores  positivos  de  X  Implican  que  F  tiende  a  cero,  con  pendiente  nula  si  X>\  y  con  pendiente 
infinita  si  X<1.  Para  el  caso  en  que  se  presenta  una  dlscontinidad  en  los  valores  aslmptótlcos 
de  F,  según  se  tienda  al  punto  crítico  desde  temperaturas  superiores  o  inferiores  a  la  crítica,  se 
asigna  X  =  0.  Finalmente,  valores  negativos  de  X  Implican  divergencias. 

A  continuación  definiremos  los  seis  exponentes  críticos  (macroscópicos)  habitualmente 
utilizados: 

a,  a'  exponentes  de  la  capacidad  calorífica: 

C(0)  ~  (0)"°,  para  0^O+ 

C(0)  ~  (-0)"° ,  para  0^0" 

en  condiciones  críticas:  v=vc  (H=Hc=0): 

(3  exponente  del  parámetro  de  orden: 
q(0)  ~  (-6)p  para  0— >0" 

en  condiciones  críticas:  v=vc  (H=Hc=0): 

Y,  y’  exponentes  de  la  compresibilidad  (o  susceptibilidad): 
kt(0)  ~  (0)"v,  para  0— >0+ 

kt(0)  ~  (-0)'v ,  para  0— >0" 
en  condiciones  críticas:  v=vc  (H=Hc=0). 

5  grado  de  la  Isoterma  crítica  (0=0) 

-tt  ~  us,  para  u— >0 

(H  ~  M5,  para  M— >0) 


Exponentes  para  el  punto  crítico  de  un  fluido  de  van  der  Waals 


Veamos  que  exponentes  críticos  resultan  para  un  fluido  que  satisfaga  la  ecuación  de  estado 


de  van  der  Waals. 
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Recordando  la  ecuación  de  estado  en  coordenadas  reducidas 
tt  =  40  -  60u  -  3u3/2 

se  concluye  inmediatamente  que  la  isoterma  crítica  (0=0)  se  comporta  como 
tt  «  -  3u3/2 
es  decir,  con  5-3. 

B 


Tomando  como  parámetro  de  orden  q  =  uG-uL,  para  hallar  p  es  necesario  considerar  la 
separación  de  materia  en  dos  fases,  por  debajo  de  Tc.  Por  la  condición  de  coexistencia  de 
fases  pL=pG,  a  T=cte  debe  anularse  la  integral  de  L  a  G 
I  dp  =  |  v.dP  =  Pcvc  |  (1  +u).du  =  Pcvc  I  u.dTr  =  0 

pues  también  son  iguales  las  presiones  tt  entre  las  fases  que  coexisten. 

1  u.drr  =  1  u(3TT/du)e.du  =  -1  u(60  +  9u2/2).du  =  -1  (60u  +  9u3/2).du  =  0 
Como  el  integrando  es  impar,  debe  ser: 

uL  =  -uG 

Es  decir  que  cerca  del  punto  crítico  los  volúmenes  molares  reducidos  son  iguales  en  módulo: 
uG  se  aparta  tanto  de  vc  en  más,  como  vL  lo  hace  en  menos;  la  curva  de  coexistencia  es 
simétrica  respecto  de  vc.  Es  importante  notar  que  como  consecuencia  de  esta  relación,  cerca 
del  punto  crítico  la  materia  se  distribuye  en  fracciones  ¡guales  entre  ambas  fases.  En  efecto,  uL 
=  -uG  Implica  vG  =  2vc-vl  y  en  virtud  e  la  regla  de  la  palanca: 

xG  =  (v-vL)/(vG  -vL)  =  (v-vl)/2(vc  -vl)  — ►  !4,  cuando  v— >vc 


Consecuentemente  xL  — ►  !4. 

Considerando  ahora  la  condición  ttl— ttg 

40  -  60ul  -  3ul3/2  =  40  -  60ug  -  3uG3/2 

60ul  +  3ul3/2  =  60ug  +  3ug3/2 

120ug  +  3ug3  =  0 

uG2  =  -40 

uG  =  2(-0)1/2 

q  =  uG-uL  =  4(-0)1/2 

con  lo  que  p  =  !4 


Y.  Y' 

Es  conveniente  considerar  la  inversa  de  la  compresibilidad: 
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1/kt  =  -v(5P/av)T  =  -Pc(1+u)(c,TT/c,u)e  =  -Pc(1+u)(66  +  9u2/2) 

i)  para  0>O  y  en  condiciones  críticas  (u=0): 

1/kt  =  -6PC.0  =  01,  cuando  0— >0+ 

con  lo  que  y=1 

ii)  para  0<O  hay  que  considerar  la  separación  de  fases  en  proporciones  xG  =  xL  =  %.  En  ambas 
fases  u2  =  -40  y  entonces: 

1/kt  =  -Pc(1+u)(60-180)  =  12Pc[1+2(-0)1/2].0  *  01,  cuando  0^0' 
con  lo  que  también  y'=1 . 


a,  a' 

La  entropía  molar  de  un  gas  de  van  der  Waals  puede  escribirse: 
s  =  cte  +  Cv°.lnT  +  R.ln(v-b) 

i)  para  0>O 

Cv  =  T.(c>s/c>T)v  =  Cv°,  constante 

ii)  para  0<O,  hay  separación  de  fases  y  aún  manteniendo  constante  el  volumen  total  v,  el 
volumen  de  cada  una  de  ellas  varía  con  T. 

s  =  (sg+sl)/2  =  cte  +  Cv°.lnT  +  (R/2).ln(vG-b)  +  (R/2).ln(vL-b) 

s  =  cte'  +  Cv°.lnT  +  (R/2).[ln(1+3uG/2)  +  ln(1+3uL/2)] 

s  =  cte'  +  Cv°.lnT  +  (R/2).ln(1-9uG2/4)  =  cte'  +  Cv°  InT  +  (R/2).ln(1+90). 

Desarrollando  en  serie  de  potencias  de  0: 
s  =  cte'  +  Cv°  InT  +  R.(90/2  -  81 02/4  +  ...) 
y  entonces: 

Cv  =  Cv°  +  RT.(9/2  -  810/2  +  ...)/Tc  =  Cv°  +  R.(9/2  -  360  +  ...) 

Cv  — >  Cv°  +  9R/2  para  0— >0" 

Es  decir  que  Cv  presenta  una  discontinuidad  finita:  ACV=  9R/2  con  lo  que  resulta  a  =  a'  =  0. 

Resumiendo,  obtenemos  para  un  fluido  de  van  der  Waals  los  exponentes  críticos: 
a  =  a'  =  0 
(3  =  y2 
Y  =  Y  =  1 
5  =  3, 

valores  denominados  “clásicos". 


Exponentes  para  el  punto  crítico  magnético  en  el  modelo  de 
Curie-Weiss 

Veamos  ahora  qué  exponentes  críticos  resultan  para  un  sistema  magnético  en  el  modelo  de 
de  Curie-Weiss.  Como  vimos,  la  ecuación  de  estado  puede  escribirse: 

M(T,H,N)  =  Np.tanh[p(H+^M)/kT] 

o  bien,  teniendo  en  cuenta  que  Ms  =  Np  y  Tc  =  ^Np2/k  =  XpMs/k, 
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M  =  Ms.tanh(|jH/kT+TcM/TMs) 

Definiendo  T  =T/TC  y  m  =M/MS: 
m  =  tanh(|jH/kT+m/7"). 

Como  cerca  del  punto  crítico  el  argumento  es  pequeño,  podemos  desarrollar  en  serle: 
m  -  (gH/kT+m/T)  -  (pH/kT+m/T)3/3 


5 

Para  la  Isoterma  crítica  (T= 1) 

(pH/kTc)  -  (pH/kTc+m)3/3  =  0 

m  •  (3pH/kTc)1/3  -  (pH/kTc)  •  (3pH/kTc)1/3  para  H-»Hc=0 
y  entonces  5  =  3 

P 

Para  H=Hc=0 
m  =  (m IT)  -  (m/T)3/3 
(1-1/7)  « -m2/373 

m2  =  -3T2(T -1)  =  -30(1+0)2  =  -30  para  0~>O 
y  entonces  (3  =  !4 


Y,  Y' 

Al  primer  orden,  para  T>TC 
m  «  pH/kT+m  IT 
m  =  pH/k(T-Tc) 

Xt  =  Np.(am/aH)T  =  Np2/k(T-Tc)  =  1/70 
y  entonces  y  =  1. 

Vemos  entonces  que  resultan  exponentes  críticos  similares  a  los  de  la  condensación  de  un 
fluido  de  van  der  Waals. 


Teoría  de  Landau 

Landau  (1937)  propuso  una  teoría  fenomenológlca  de  las  transiciones  de  fase  continuas  a 
partir  del  cambio  de  simetría  en  el  punto  de  transición.  Como  vimos,  al  bajar  la  temperatura  por 
debajo  de  Tc  hay  una  ruptura  de  simetría  en  el  sistema  y  se  requiere  de  una  variable  adicional, 
parámetro  de  orden  q,  para  describir  sus  estados  de  equilibrio. 

-  A  T>TC  el  estado  de  equilibrio  consiste  de  una  fase  única  simétrica  (q  =  0) 

-  A  T<TC  existen  varias  fases  no  simétricas  posibles  (q  #  0) 
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-  Cuando  T— >T0,  n  — *-0  continuamente. 

Suponiendo  que  se  trata  de  un  sistema  magnético,  q  =  M,  intentaremos  reflejar  estas 
características  de  los  estados  de  equilibrio  a  través  de  un  potencial  termodinámico 

O  =  0(T,H;M) 

en  el  que  M  tiene  un  carácter  diferente  como  variable.  Mientras  T  y  H  pueden  darse 
arbitrariamente,  el  valor  de  M  queda  determinado  por  la  condición  de  equilibrio:  mínimo  O  para 
los  valores  de  T  y  H. 

Para  estos  valores  de  M:  0(T,H;Meq)  =  G(T,H).  Supondremos  que,  en  torno  al  punto  crítico, 
este  potencial  puede  desarrollarse  en  serle  de  potencias  de  M 

<D(T,H;M)  =  Oo  +  aM  +  AM2  +  BM3  +  CM4  +  ... 

donde  0O,  a,  A,  B,  C  son  funciones  de  T  y  H.  Las  principales  propiedades  estarán  contenidas 
en  el  término  cuadrático  y  su  dependencia  con  la  temperatura.  Como  dijimos,  los  valores  de  M 
resultan  de  las  condiciones  de  equilibrio: 

(<90A9M)th  =  a  +  2AM  +  3BM2  +  4CM3  =  0 
(<920/<9M2)th  =  2A  +  6BM  +  12CM2  >  0. 

En  condiciones  críticas  (H=0),  debe  haber  una  solución  con  M=0  para  T>TC  y  ello  sólo  es 
posible  si  a=0.  A  su  vez,  para  T<TC,  las  soluciones  con  M^O,  deben  ser  simétricas  y  entonces 
debe  ser  B=0  además  de  a=0.  Aunque  un  término  lineal  de  la  forma  -HM,  podría  satisfacer 
estos  requisitos  también.  Tendremos  entonces: 

0(T,H;M)  =  0o(T,H)  -  HM  +  A(T,H)M2  +  C(T,H)M4 

que  constituye  la  base  de  la  teoría  de  Landau.  De  aquí  se  obtienen,  para  H=Hc=0: 

(<90/<9  M)th  =  2AM  +  4CM3  =  0 
(<920/<9M2)th  =  2 A  +  12CM2  >  0. 

con  soluciones: 

M  =  0;  <920/<9M2  =  2A  =>  A  >  0. 

M  =  ±(-A/2C)1/2;  <920/<9M2  =  -4A  =>  A  <  0. 
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adecuadas  para  T>Tcy  T<TC,  respectivamente.  A  T=TC,  punto  crítico,  resulta  un  mínimo  chato 
( d2O/dM2=0 )  del  potencial,  donde  coalescen  las  dos  soluciones  de  bajas  temperaturas  y  debe 
ser  A=0.  En  ese  punto  el  equilibrio  requiere: 

(330A9M3)th  =  24CM  =  0,  que  se  satisface  pues  M  =  0  y 
(^cD/aM4)™  =  240  O  =>  C  >  0. 

Por  lo  tanto  una  expresión  como 
A  =  a(T-Tc)/Tc  =  a6 

con  a=cte.  asegurará  las  condiciones  anteriores.  Por  su  parte,  el  requerimiento  de  C  >  O  puede 
extenderse  a  las  temperaturas  cercanas  a  Tc.  De  hecho,  asumiremos  que  C  no  depende  de  la 
temperatura.  Tenemos  entonces: 


O  =  0o-HM  +  a.6M2  +  CM4 


con  a  y  C  constantes  positivas. 

Veamos  qué  exponentes  críticos  resultan  en  esta  teoría. 


P 

Para  H=Hc=0 

M  =  ±(-A/2C)1/2  =  ±(a/2C)1/2(-0)1/2  => 

es  decir  p  =  !4 


a,  a' 

La  entropía  puede  hallarse  como  -30/dT  para  los  valores  de  M  de  equilibrio.  La 
dependencia  explícita  de  M  con  T  puede  ignorarse  al  derivar,  ya  que  ya  30/3M  =  0. 

Luego  podemos  hallar  CH,  derivando  una  vez  más. 

S  =  -(00/0T)h  =  -  (30O/3T)  -  {8Aldl)M2  =  S0  -  (a/Tc)M2 

i)  para  0>O;  M  =  O 
S  =  S0 

CH  =  T(3S0/3T)h  =  CH,o 

ii)  para  0<O;  M2  =  -a0/2C 
S  =  S0  +  (a20/2CTc) 

CH  =  T(3S/3T)h  =  CH,o  +  (a2T/2CTc2)  -  CH,0  +  (a2/2CTc) 

Vemos  así  que  CH  presenta  una  discontinuidad  finita  (a2/2CTc)  por  lo  que  asignamos: 
a  =  a'  =  O 
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5 

Para  hallar  los  exponentes  5,  y  y  Y'  es  necesario  explicitar  la  dependencia  de  H  con  M 

(ecuación  de  estado).  Como 

(<90/<9M)th  =  -  H  +  2a6M  +  4CM3  =  0 

H  =  2a6M  +  4CM3 

Para  la  Isoterma  crítica  (0  =0) 

H  =  4CM3 
y  entonces  5-3. 


Y,  Y' 

La  Inversa  de  la  susceptibilidad  a  H=Hc=0  será: 

1/Xt  =  (5H/5M)t  =  2a0  +  12CM2 

i)  para  0>O;  M  =  0 

1/Xt  =  2a0  =>  y  =  1 

ii)  para  0<O;  M2  =  -a0/2C 
1/Xt  =  -4a0  =>  y'  =  1 

Así  podemos  ver  que  la  teoría  de  Landau,  basada  en  criterios  generales  sobre  la  simetría  de 
las  soluciones  en  torno  al  punto  crítico,  predice  el  mismo  conjunto  de  valores  clásicos  de  los 
exponentes  críticos. 
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CAPÍTULO  11 
Hipótesis  de  escala 


En  este  capítulo  confrontaremos  las  predicciones  sobre  el  comportamiento  crítico  obtenidas 
en  el  capítulo  anterior  con  los  resultados  experimentales,  comprobando  diferencias 
cuantitativas  sistemáticas.  Para  subsanar  estas  discrepancias  formularemos  un  postulado 
adicional  sobre  la  homogeneidad  de  las  funciones  termodinámicas  en  la  vecindad  del  punto 
crítico. 


Valores  experimentales  de  los  exponentes  críticos 


Hemos  visto  cómo  en  las  teorías  clásicas  (campo  medio)  resultan  exponentes  críticos 
similares  para  sistemas  tan  diversos  como  un  fluido  que  condensa  o  un  paramagneto  que  se 
ordena  magnéticamente.  Incluso,  una  teoría  general  como  la  de  Landau  predice  los  mismos 
exponentes  críticos.  Estos  valores,  sin  embargo,  no  están  de  acuerdo  con  los  observados 
experimentalmente. 

La  determinación  experimental  de  exponentes  críticos  es  en  general  dificultosa  ya  que  requiere 
estabilizar  al  sistema  en  estados  de  equilibrio,  de  temperatura  definida,  en  condiciones  en  que 
las  funciones  respuesta  del  mismo  divergen  y  se  producen  grandes  fluctuaciones  en  las 
magnitudes  macroscópicas.  Aún  así,  a  partir  de  los  años  50  se  fueron  reportando  valores 
precisos,  obtenidos  en  muy  diversos  sistemas,  que  mostraban  un  apartamiento  sistemático  de 
los  predicóos  y,  a  su  vez,  una  coincidencia  entre  ellos:  universalidad.  Por  ejemplo,  para  el 
exponente  (3  del  parámetro  de  orden,  en  la  condensación  de  diversos  fluidos  (Guggenheim, 
1945),  en  el  ordenamiento  antlferromagnétlco  del  MnF2  (Heller  y  Benedek,  1962)  y  en  la 
separación  de  fases  de  la  mezcla  binarla  líquida  CCI4+C7Fi4  (Thomson  y  Rice,  1964),  se 
determinaron  valores  cercanos  a  1/3  en  todos  los  casos  (ver  Tabla15.1).  En  lo  que  sigue 
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intentaremos  comprender  las  razones  de  esta  discrepancia  entre  teoría  y  experimento  y 
formular  las  hipótesis  adecuadas  para  corregirla. 

Un  aporte  importante  en  este  sentido  surge  de  cálculos  basados  en  modelos  que  describen 
las  interacciones,  no  a  través  de  un  campo  medio,  sino  a  partir  de  la  interacción  entre  pares  de 
primeros  vecinos  (modelo  generalizado  de  Heisemberg  y  derivados).  El  hamiltoniano  del 
modelo  puede  escribirse  como: 

H(d)  =  -J.Z<ij>Si(D).Sj(D) 

donde  los  espines  S¡(D)  son  vectores  unitarios  de  dimensión  D  y  -J  es  la  energía  de  un  par  de 
espines  paralelos  <¡j>,  vecinos  próximos  localizados  en  los  sitios  i  y  j  de  la  red.  Existe 
numerosa  bibliografía  sobre  este  tema,  propio  de  la  mecánica  estadística,  a  la  cual  el  lector 
puede  remitirse.  Aquí  nos  limitaremos  a  mencionar  algunos  de  los  principales  logros  de  estas 
teorías. 


Tabla  15.1:  Exponentes  críticos 


Valores 

a,  a' 

P 

Y.  Y' 

5 

clásicos 

0 

y2 

1 

3 

experimentales 

0-0.2 

0. 3-0.4 

1 

4-5 

modelos  nn 

0 

1/3 

4/3 

5 

En  primer  término,  como  se  aprecia  en  la  tabla,  las  predicciones  de  este  modelo  para  los 
sistemas  considerados  coinciden  con  la  evidencia  experimental.  Por  otra  parte,  modelos 
derivados,  donde  puede  variarse  la  dimensionalidad  n  del  parámetro  de  orden  (escalar,  vector 
en  2  o  3  dimensiones)  y  la  dimensión  d  del  sistema,  muestran  la  existencia  de  grupos  de 
fenómenos  con  un  mismo  conjunto  de  exponentes  críticos:  clases  de  universalidad,  lo  cual  se 
verifica  experimentalmente.  A  continuación  presentamos  estas  clases  con  los  modelos  y 
sistemas  asociados  en  una  tabla  adaptada  de  Kenneth  Wilson. 
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Tabla  15.2:  Clases  de  universalidad 


clase 

modelo  teórico 

sistema  físico 

parámetro  de  orden 

d 

n 

2 

i 

Ising  en  2d 

películas  adsorbidas 

densidad  superficial 

2 

2 

XY  en  2d 

películas  de  Helio-4 

amplitud  de  fase  superfluida 

2 

3 

Heisemberg  en  2d 

magnetización 

3 

0 

camino  al  azar  sin 

conformación  de 

densidad  de  finales  de  cadena 

intersecciones 

polímeros  cadena  larga 

3 

1 

Ising  en  3d 

ferroimán  uniaxial 

magnetización 

3 

1 

“ 

fluido 

diferencia  de  densidad 

3 

1 

“ 

mezcla  binaria 

diferencia  de  concentración 

3 

1 

" 

aleación 

orden  de  largo  alcance 

3 

2 

XY  en  3d 

ferroimán  plano 

magnetización 

3 

2 

11 

Helio-4 

3 

3 

Heisemberg  en  3d 

ferroimán  isótropo 

magnetización 

4 

32 

cromodinámica  cuántica 

quarks  ligados 

Relaciones  termodinámicas  entre  exponentes  críticos 


Al  avanzar  la  investigación  sobre  exponentes  críticos  se  comprobó  que  satisfacían  una  serle 
de  relaciones  entre  ellos.  La  primera  y  más  conocida  es  la  desigualdad  de  Rushbrooke  (1963) 
que  puede  probarse  mediante  argumentos  termodinámlcos  a  partir  de  los  criterios  de 
estabilidad. 

En  un  sistema  magnético  se  cumple,  como  vimos: 

CH  =  CM  +  T[(<9M/<9T)h]2/Xt 
y  como  debe  ser  C  >  0,  entonces 

M 

CH  >  T[(3M/3T)h]2/Xt 

En  condiciones  críticas:  H=0  y  T— *TC',  puede  escribirse: 

CH  ~  (-6)"°  ;  (3M/<9T)h  ~  (-0)3"1;  Xi  ~  (-6)  Y 


y  entonces: 

A(-6)"a  >  BTC  [(-e)2p  2+Y’](1+e). 
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Teniendo  en  cuenta  la  definición  rigurosa  de  exponente  crítico,  esta  desigualdad  implica: 


a'+2(3+Y'  £  2 

que  es  la  desigualdad  de  Rushbrooke.  Sorprendentemente,  ésta  parece  verificarse  como  una 
igualdad,  tanto  para  los  valores  experimentales  de  los  exponentes  críticos  como  para  aquellos 
obtenidos  con  teorías  de  campo  medio  o  las  más  exitosas  de  primeros  vecinos. 

Griffiths  (1965)  pudo  demostrar  varias  desigualdades  más,  entre  ellas, 

cf+p(1+5)>2, 

basado  en  la  convexidad  de  la  energía  libre  y 
y'>p(5  -1). 

suponiendo  que  la  curva  de  magnetización  M(H)  debe  ser  cóncava  (susceptibilidad  decreciente 
con  H). 

Todas  estas  desigualdades  parecen  verificarse  también  como  igualdades.  De  hecho,  esta 
última,  previamente  conjeturada  por  Widom  (1964)  como  igualdad,  podría  deducirse  a  partir  de 
las  dos  anteriores,  tomadas  como  igualdades.  Por  otro  lado,  parecen  ser  iguales  los 
exponentes  primados  y  no  primados,  es  decir: 

a=a';  y=Y  ■ 

Es  interesante  notar  que,  de  satisfacerse  estas  relaciones  como  Igualdades,  el  conjunto  de 
6  exponentes  críticos,  para  cada  clase  de  universalidad,  quedaría  reducido  a  2  exponentes 
independientes.  Podríamos  escribir,  por  ejemplo: 


a  =  a'  =  2-p(5+1) 

Y  =  y'  =  P(5-1). 

conocidas  como  relaciones  de  escala  de  los  exponentes  críticos. 


El  rol  de  las  fluctuaciones 

Antes  de  formular  la  hipótesis  de  escala,  es  conveniente  razonar  sobre  el  estado  físico  en 
que  se  encuentra  un  sistema  en  la  vecindad  del  punto  crítico.  Nos  preguntamos  porqué  falla  en 
esta  situación  la  termodinámica  clásica,  encarnada  en  la  teoría  general  de  Landau  o  en  teorías 
específicas  de  campo  medio.  La  razón  parece  estar  en  el  hecho  de  que  un  sistema,  cerca  del 
punto  crítico,  se  encuentra  dominado  por  las  fluctuaciones  termodinámicas. 
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Por  ejemplo  en  un  fluido,  la  compresibilidad  diverge,  dando  lugar  a  grandes  cambios  locales 
en  la  densidad,  que  demoran  en  recuperarse.  Podemos  Imaginar  al  sistema  atravesado  por 
variaciones  continuas  de  la  densidad,  desde  valores  propios  del  líquido  hasta  valores  típicos 
del  gas,  reflejando  la  inminencia  de  la  transición  de  fase.  Esta  situación  se  manifiesta 
claramente  en  el  fenómeno  de  opalescencia  crítica.  Un  tubo  conteniendo  un  fluido  (C02)  a  vc,  al 
disminuir  su  temperatura  justo  por  debajo  de  Tc,  dispersa  la  luz  de  cualquier  longitud  de  onda  y 
se  vuelve  opalescente,  como  consecuencia  de  las  variaciones  en  el  índice  de  refracción 
asociadas  a  las  fluctuaciones  en  la  densidad. 

Estas  fluctuaciones  termodinámicas,  de  todo  tamaño,  dominan  entonces  el  comportamiento 
del  sistema.  Las  teorías  de  campo  medio  Ignoran  por  completo  la  presencia  de  fluctuaciones, 
asumiendo  un  caréete  uniforme  para  el  sistema  y  describiendo  la  Interacción  entre  partículas 
como  un  efecto  promedio.  De  modo  que  no  es  sorprendente  que  estos  modelos  fallen  al 
describir  los  estados  cercanos  al  punto  crítico.  Más  aún,  podríamos  preguntarnos  si  la 
termodinámica,  tal  como  la  hemos  formulado,  es  capaz  de  describir  correctamente  tales 
estados.  Nuestro  postulado  II  identifica  el  estado  de  equilibrio  como  el  de  máxima  entropía.  En 
rigor,  el  valor  de  equilibrio  observado  de  una  dada  propiedad  corresponde  a  un  valor  promedio 
según  probabilidades  que  estudia  la  mecánica  estadística.  Como  la  función  de  distribución  es 
muy  aguda  el  promedio  coincide  aproximadamente  con  el  valor  más  probable.  Este  último  es 
más  fácil  de  calcular  pues  es  el  que  maxlmlza  la  probabilidad  exp[S/k]  y  por  tanto  maxlmlza  S. 
Nuestro  postulado  Identifica  incorrectamente  el  valor  mas  probable  con  el  valor  de  equilibrio  o 
promedio.  En  un  sistema  dominado  por  fluctuaciones  esto  no  será  necesariamente  cierto  ya 
que  la  distribución  se  vuelve  ancha  y  asimétrica  y  es  de  esperar  que  nuestras  predicciones 
fallen,  al  menos  cuantitativamente. 


Hipótesis  de  escala 


Para  subsanar  las  discrepancias  observadas  postularemos  propiedades  adicionales  para  los 
potenciales  termodlnámlcos  (Wldom,  1965).  Suponemos  que  cerca  del  punto  crítico  la  energía 
libre  de  Glbbs  consiste  de  una  parte  regular  G0  y  otra  singular  G  que  describe  la  física  relevante 
para  la  transición  de  fase.  La  hipótesis  de  escala  estática  consiste  en  suponer  que 

La  parte  singular  de  la  energía  libre  de  Gibbs  es  una  función  homogénea  generalizada  de  sus 
variables  intensivas. 

Es  decir,  para  un  sistema  magnético: 

G(Xa0,XbH)  =  «3(6, H) 


84 


con  0  =  (T-Tc)/Tc.  Los  exponentes  a  y  b  se  denominan  parámetros  de  escala.  Es  importante 
notar  que  la  hipótesis  no  especifica  el  valor  de  estos  parámetros,  los  que  podrán  expresarse, 
como  veremos,  en  término  de  los  exponentes  críticos  del  sistema. 

Derivando  respecto  de  H  obtenemos  para  M: 


XbM  (A.a0,XbH)  =  XM  (0,H) 


o  bien 

M(|ia0,  pbH)  =  pM(0,H) 

es  decir  que  M  (su  parte  singular)  es  también  una  función  homogénea  generalizada  con 
parámetros  de  escala:  a'=a/(1-b)  y  b'=b/(1-b).  De  hecho,  la  parte  singular  de  cualquier  derivada 
de  G  o  de  cualquiera  de  sus  transformadas  de  Legendre,  será  también  función  homogénea 
generalizada  de  sus  variables  intensivas. 

Veamos  ahora  cómo  expresar  los  exponentes  críticos  del  sistema  en  términos  de  los 
parámetros  a  y  b. 


P 

Para  H=Hc=0  y  0^0' 
ibM(^a0,O)  =  xM(0,O) 

=>  M(0,O)  =  ^b_1M(xa0,O) 

eligiendo  X  =  (-0)"1/a 

M(0,O)  =  (-0)<1_b)/aM(-1,O)  ~  (-0)p 

con  lo  que  (3  =  (1-b)/a 


5 

Para  0=0  y  H->0 
XbM(0,XbH)  =  XM(0,H) 

=>  M(0,H)  =  A.b_1M(0,>.bH) 
eligiendo  X  =  H"1/b 
M(0,H)  =  H(1_b)/bM(0,1)~  H1/5 
con  lo  que  5  =  b/(1-b) 


Y 

Derivando  G  dos  veces  respecto  de  M,  para  H=Hc=0  y  0— >0 
^XT^O.O)  =  ^Xi(e.o) 
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eligiendo  X  =  0  \ 

xt(0.o)  =  lel"(2b"1)/aXT(±i,o) ~  le  I™ 

con  lo  que  y  =  y'  =  (2b-1)/a 


a 

Derivando  G  dos  veces  respecto  de  0,  para  H=Hc=0  y  0— >0 

^2aCH(^a0,O)  =  XCH(0,O) 

eligiendo  X  =  |  0  |  "1/a 

CH(0,O)=  |  0  | -(2a-1)/aCH(±1 ,0)  ~  |0|_a'_a' 

con  lo  que  a  =  a'  =  (2a-1)/a 

Entonces  M  satisface  la  relación  de  homogeneidad  generalizada: 

M(|i1/p0,  |iSH)  =  pM(0,H) 

Además  podemos  identificar: 

b  =  5/(5+ 1 ) 
a  =  1/p(5+1) 

con  lo  que  la  hipótesis  de  partida  puede  reescribirse: 

G(a1/P(6+1)0,a6/(5+1)H)  =  XG(0,H) 

Ahora  que  hemos  expresado  todo  los  exponentes  críticos  en  términos  de  los  parámetros  de 
escala,  es  fácil  verificar,  que  las  desigualdades  termodinámicas  (Rushbrooke  y  Griffiths)  se 
satisfacen  como  igualdades,  en  el  marco  de  esta  hipótesis.  Por  ejemplo: 
a +2p+y’  =  (2a-1)/a  +  2(1-b)/a  +  (2b-1)/a  =  2a/a  =  2 


Ecuación  de  estado  reducida 

A  partir  de  la  hipótesis  de  escala  podemos  deducir  la  forma  de  la  ecuación  de  estado  para 
M(0,H)  cerca  del  punto  crítico.  Partiendo  de  la  relación  de  homogeneidad: 

nM(0,H)  =  M(|i1/p0,  p5H) 

si  elegimos  p  =  0  |  "p 


86 


Definiendo  magnitudes  escaleadas: 


m  =  M/ 1 6 1 p  h  =  H/ 1 6 1 p5 

resulta  una  ecuación  de  estado  universal  para  la  variación  de  M  con  H,  convenientemente 
escaleados,  a  diferentes  T. 
m  =  m(±1,h)  =  m+(h). 

Esta  dependencia  describe  una  curva,  donde  coalescen  medidas  de  la  relación  MvsHa 
distintas  temperaturas  (Stanley,  1999).  Alternativamente  eligiendo  p  =  H"1/6  puede  analizarse  la 
relación  de  M  con  T  para  diferentes  valores  de  H: 

M(0,H)/H1/5  =  M(0/H1/S|3,  1) 

Definiendo  magnitudes  escaleadas: 

m'  =  M/H1/6  t  =  0/H 1/5,3 

resulta  una  ecuación  de  estado  universal: 
m'  =  m'(t) 

que  también  se  verifica  notablemente  mostrando  la  coincidencia  de  curvas  M  vs  T,  de  distintos 
sistemas  magnéticos  para  diferente  valores  de  H  (Stanley,  1999). 
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CAPITULO  12 

Sistemas  superconductores 


En  este  capítulo  aplicaremos  nuestro  formalismo  a  sistemas  superconductores,  limitándonos 
a  describir  sus  propiedades  macroscópicas  y  a  analizar  las  características  de  la  transición  de 
fase  involucrada. 


Fenomenología  de  la  superconductividad 


Numerosos  metales,  aleaciones  y  compuestos  pasan,  al  bajar  la  temperatura,  a  un  estado 
cuyo  propiedad  más  espectacular,  descubierta  en  1911  por  Kamerlingh  Onnes,  es  la 
superconductividad,  es  decir  la  ausencia  de  resistencia  eléctrica  (p=0)  al  pasaje  de  una 
corriente  continua.  Cuando  estudiaba  la  resistencia  del  mercurio  a  muy  bajas  temperaturas, 
encontró  que  ésta  caía  a  valores  casi  nulos  por  debajo  de  4.2K. 

En  general  esperamos  que  la  resistividad  eléctrica  de  un  metal  disminuya  con  la 
temperatura  como  resultado  de  la  reducción  en  las  vibraciones  de  red  (fonones),  tendiendo  a 
OK,  a  un  valor  finito  de  resistividad  residual  originado  en  la  presencia  de  Impurezas  y  otros 
defectos  estáticos.  Sin  embargo  los  materiales  superconductores,  a  una  cierta  temperatura 
crítica  Tc,  sufren  una  transición  abrupta  a  un  estado  de  resistividad  nula.  Los  valores 
observados  de  Tc  a  partir  del  descubrimiento  de  la  superconductividad  han  alcanzado  valores 
cada  vez  más  altos,  por  encima  de  150K,  ampliando  las  posibilidades  concretas  de  aplicación 
de  este  fenómeno. 

Una  vez  establecida  una  corriente  en  un  anillo  superconductor,  ésta  podría  entonces 
persistir  indefinidamente:  conductor  perfecto  (o= oo).  No  obstante,  no  es  la  ausencia  de 
resistencia  eléctrica  la  propiedad  fundamental  de  los  superconductores.  Son  las  propiedades 
magnéticas  las  que  sufren  los  cambios  más  profundos  en  la  transición.  El  campo  magnético  no 
penetra  nunca  en  un  superconductor;  luego  B=0  en  su  interior:  diamagnetismo  perfecto.  SI 
enfriamos  un  superconductor  en  presencia  de  un  campo  aplicado,  en  el  momento  de  la 
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transición  las  líneas  son  expulsadas.  Propiedad  descubierta  por  Meissner  y  Ochsenfeld  en 
1933  y  que  hoy  se  conoce  como  efecto  Meissner. 

Este  efecto  permite  distinguir  entre  un  diamagneto  perfecto  y  un  conductor  perfecto, 
mostrando  que  la  propiedad  fundamental  de  un  superconductor  es  la  primera.  En  ambos  casos, 
si  enfriamos  el  cuerpo  por  debajo  de  Tc  y  luego  aplicamos  un  campo  externo  sus  líneas  no 
penetrarán  en  el  cuerpo.  Si  primero  aplicamos  el  campo  y  luego  enfriamos,  las  líneas  serán 
expulsadas,  en  el  primer  caso,  conduciendo  al  mismo  estado  final.  En  cambio,  en  el  segundo 
caso,  no  habrá  expulsión  pues: 

(oB/ot)  =  -  c.rot  E  y  E  =  j/o 
y  si  o— >oo,  rotE  — >  O  y  B  permanece  constante. 

En  rigor  B=0  no  se  cumple  en  una  delgada  capa  superficial  del  cuerpo  cuyo  espesor  se 
denomina  longitud  de  penetración  A,  lO-IOOnm.  Por  otro  lado  la  expulsión  es  completa  sólo  en 
superconductores  puros  del  tipo  I.  El  campo  puede  penetrar  parcialmente  formando  tubos  de 
flujo  cuantificado  en  los  de  tipo  II. 

El  análisis  microscópico  de  la  superconductividad  requiere  de  un  tratamiento  cuántico.  En 
primer  lugar  se  ha  comprobado  que  los  portadores  de  carga  no  son  electrones  simples  como 
en  el  estado  normal,  sino  pares  de  electrones  de  espines  opuestos,  acoplados  para  formar  los 
llamados  pares  de  Cooper,  de  espin  nulo.  El  acoplamiento  entre  ambos  se  da  a  través  de  las 
vibraciones  de  la  red  (fonones),  mediante  una  interacción  electrón-fonón-electrón  que  es  la 
base  de  la  teoría  BCS  (por  sus  autores  J.  Bardeen,  L.N.  Cooper  y  J.R.  Schrlefer)  cuyo 
tratamiento  excede  el  marco  de  este  curso. 


Corriente  y  campo  críticos 


Consideremos  la  superficie  límite  de  un  medio  superconductor.  En  la  frontera  entre  dos 
medios  la  componente  normal  del  campo  de  inducción  magnética  Bn  debe  ser  continua  pues 
divB  =  0.  En  este  caso  Bnext  =  Bnnt  =  0,  es  decir  que  el  campo  exterior  será  tangente  a  la 
superficie  del  cuerpo:  las  líneas  de  fuerza  envuelven  al  superconductor. 

Por  otra  parte,  rotB  =  0  Implica  que  la  densidad  media  de  corriente  en  el  superconductor  es 
nula,  es  decir  que  las  corrientes  macroscópicas  volumétricas  son  imposibles  en  él;  toda 
corriente  será  superficial.  La  densidad  de  corriente  g[A/m]  puede  calcularse  a  partir  de  la 
discontinuidad  en  B,  en  la  superficie  de  separación  entre  el  cuerpo  y  el  exterior: 

g  =  (c/4tt)  n/\[Bext  -  Bint]  =  (c/4tt)  nAH 

donde  H  es  el  campo  en  el  exterior  suponiendo  p=1 . 

Esta  corriente  superficial  tiene  un  efecto  de  apantallamiento  y  es  la  responsable  de  que  el 
campo  en  el  interior  se  anule.  Puesto  el  superconductor  en  presencia  de  un  campo  H,  se 
establece  g  superficial  de  modo  de  anular  B  en  su  interior.  Cuando  el  campo  aplicado  es 
suficientemente  grande  puede  destruir  la  superconductividad.  Para  H>HC  el  campo  penetra  en 
el  cuerpo  y  éste  se  vuelve  normal.  A  T=TC,  cualquier  campo  finito  destruye  la 
superconductividad,  Hc=0,  y  para  T<TC  cuanto  menor  sea  la  temperatura  mayor  será  el  campo 
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crítico.  Experimentalmente  se  ha  determinado  para  la  mayoría  de  los  superconductores,  la 
siguiente  dependencia,  mostrada  en  la  figura  12.1: 

Hc(T)  =  H0[1-(T/Tc)2] 
con  H0  =  Hc(0). 

Esta  curva  de  campos  y  temperaturas  críticos  interdependientes,  determina  en  el  plano 
(H,T)  las  regiones  en  las  que  el  cuerpo  es  superconductor  o  normal,  pudiendo  coexistir  ambos 
estados  sobre  la  misma.  Notemos  que  la  existencia  de  un  campo  crítico  determina  también  una 
corriente  crítica,  máxima  que  puede  circular  sin  destruir  la  superconductividad:  gc  =  cHc/4tt. 


Figura  12.1 


Propiedades  magnéticas 

Si  bien  en  el  interior  del  superconductor  no  tiene  sentido  definir  M  y  H,  podemos 
introducirlos  de  una  manera  formal  (B=H+4ttM),  para  calcular  el  momento  magnético  total  M  = 
VM  que  conserva  su  sentido  físico  habitual.  Dado  que  B=0,  resulta: 

M  =  -VH/4tt 

como  si  el  cuerpo,  en  el  estado  superconductor  respondiera  con  una  susceptibilidad  x--1/4tt. 
En  tanto,  en  el  estado  normal  podemos  asumir  x=0  y  entonces  B=H  (ver  figura  12.2). 

En  consecuencia,  podemos  hacer  nuestra  descripción  termodinámica  considerando  al 
cuerpo  superconductor  como  un  diamagneto  uniforme  que  en  presencia  de  un  campo  externo 
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H  desarrolla  un  momento  magnético  -VH/4tt.  La  curva  HC(T)  en  el  plano  (H,T)  lo  dividirá  en  dos 
reglones  donde  existen  las  fases  superconductora  y  normal,  respectivamente.  La  propia  curva 
HC(T)  será  la  curva  de  equilibrio  o  coexistencia  de  fases.  En  la  figura  12.3  se  muestran  los 
diagramas  B  vs  H  en  esta  descripción. 


Figura  12.2 


Figura  12.3 
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Energía  de  condensación 


Vamos  a  formular  un  modelo  macroscópico  del  estado  superconductor,  basado  en  su 
respuesta  magnética,  con  el  cual  podremos  predecir  algunas  de  sus  propiedades  más 
Importantes.  Consideremos  un  cuerpo  superconductor  de  forma  cilindrica  en  presencia  de  un 
campo  longitudinal  H.  En  el  estado  superconductor  tendrá  una  energía  magnética 
suplementaria  con  respecto  al  estado  normal  dada  por: 

-jH.dM  =  (V/4TT)jHdH  =  VH2/8tt. 


Hemos  supuesto  que  en  el  estado  normal  la  energía  no  cambia  por  aplicación  de  un  campo 
magnético,  es  decir  despreciamos  el  para  o  dlamagnetlsmo  que  pudiera  tener  el  cuerpo;  p=1. 
Es  por  esta  energía  adicional  que  al  aumentar  el  campo  el  estado  superconductor  podrá 
resultar  menos  favorable  que  el  normal  al  superar  Hc. 

Para  T<TC,  el  cuerpo  puede  ser  ,  normal  (n)  o  superconductor  (s),  según  el  valor  de  H. 
Podemos  escribir  su  energía  libre  de  Helmholtz,  en  cada  caso 

Fn(T,V)  =  F  n  0(T  ,V) 

FS(T,V)  -  FS0(T,V)  +  VH2/8tt 

donde  el  subíndice  (o)  indica  el  valor  en  ausencia  de  campo.  Derivando  obtenemos 
PS(T,V)  =  -(oF/oV)t  =  P0(T,V)  -  H2/8tt 

que  es  la  presión  que  actúa  sobre  el  superconductor.  Podemos  obtener  entonces  la  energía 
libre  de  Gibbs7 

GS(T,P)  =  FS(T,V)  +  PSV  =  Fs,0  +  VH2/8tt  +  (P0-H2/8tt)V 
GS(T,P)  =  Fso  +  P0V  =  GS0(T,P+H2/8tt) 
pues  P0  =  P+H2/8tt. 

Es  decir  que  la  energía  libre  de  Gibbs  del  superconductor  en  presencia  de  un  campo  H  es  la 
misma  que  en  ausencia  de  campo  pero  a  una  presión  superior  en  H2/8tt.  Lo  mismo  ocurrirá  con 
la  entropía  y  el  volumen  que  serán  entonces  menores.  Esto  último  es  consistente  con  la  ¡dea 
de  una  contracción  de  volumen  producida  por  las  corrientes  superficiales.  En  efecto,  derivando 
el  potencial  de  Gibbs: 

Ss(T,P)  =  Sso(T,P+H2/8tt) 

V,(T,P)  =  Vso(T,P+H2/8tt) 


7  Esta  es  la  energía  libre  de  Gibbs  convencional,  función  de  T  y  P  y  no  la  que  usamos  en  sistemas  magnéticos,  que  era 
función  de  T  y  H. 
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A  su  vez,  para  el  estado  normal,  la  presencia  del  campo  no  produce  ningún  cambio  y 
podemos  escribir: 

Gn(T,P)  =  Gno(T  ,P) 

Finalmente,  para  T  y  P  dados,  si  H=HC,  debe  cumplirse  Gs=Gn  pues  ambas  fases  coexisten,  es 
decir: 


Gs,o(T,P+Hc2/8tt)  =  Gn,0(T,P) 

Esta  es  una  ecuación  termodinámica  exacta.  En  general,  el  término  Hc2/8tt  es  sólo  una 
pequeña  corrección  a  la  presión,  de  modo  que  podemos  desarrollar  en  serle  y  escribir  la 
expresión  aproximada: 


GS,0(T,P)  +  (Hc2/8tt).Vs,0(T,P)  =  Gn,0(T,P) 


Es  decir  que  en  ausencia  de  campo,  el  potencial  por  unidad  de  volumen  en  el  estado 
superconductor  es  menor  en  Hc2/8tt  que  en  el  estado  normal. 

Esta  ganancia  de  energía,  al  pasar  del  estado  normal  al  superconductor,  se  denomina 
energía  de  condensación.  Al  aplicar  un  campo  H  al  superconductor,  esta  energía  Irá 
aumentando  en  H2/8tt  hasta  que  para  H=HC,  Gs=Gn. 


Algunas  propiedades 

Sobre  la  curva  de  coexistencia  Gs=Gn,  es  decir: 

GS(T,P)  =  GS0(T,P+Hc2/8tt)  =  Gn(T,P) 

Derivando  respecto  de  T,  a  P=cte 

-Sso(T,P+Hc2/8tt)  +  Vs,0(T,P+Hc2/8TT).d(Hc2/8rr)/dT  =  -Sn(T,P) 
entonces: 

Sn(T,P)  -  Ss(T,P)  =  -  (1/8rr)Vs(T,P).dHc2/dT 

Con  lo  que  se  obtiene  una  expresión  para  la  discontinuidad  en  la  entropía8 : 

Sn  -  Ss  =  -  (VsHc/4Ti).dHc/dT. 

Recordando  la  forma  de  la  curva  HC(T),  vemos  que  dHc/dT<0  en  todo  el  rango  0-Tc,  de  modo 
que  Sn<Ss:  el  estado  superconductor  es  más  ordenado.  En  particular: 


8  Es  interesante  notar  que  esta  ecuación  es  análoga  a  la  ecuación  de  Clausius  &  Clapeyron  para  la  curva  de  equilibrio 
en  las  variable  T  y  H,  involucrando  las  discontinuidades  en  S  y  M,  propias  de  una  transición  de  fase  de  primer  orden. 
En  efecto:  AM  =  Mn-Ms  =  VsHo/4n 
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a  T=0,  como  dHc/dT=0,  Sn  =  Ss  en  acuerdo  con  el  Postulado  III 
a  T=TC,  como  Hc=0,  Sn  =  Ss,  desapareciendo  la  discontinuidad. 

Es  decir  que  entre  0  y  Tc  habrá  calor  latente  y  la  transición  será  de  primer  orden.  En  tanto  que  a 
T=TC,  desaparecen  las  discontinuidades  y  la  transición  de  fase  será  continua. 

Volviendo  a  derivar  respecto  de  T,  a  P=cte,  sobre  la  curva  de  coexistencia,  podremos 
obtener  la  diferencia  de  calores  específicos  entre  ambas  fases.  La  diferencia  Sn(T,P)  -  Ss(T,P) 
puede  escribirse: 

Sn(T,P)  -  Ss  o(T,P+Hc2/8tt)  =  -Vs,0(T,P+Hc2/8rr).d(Hc2/8TT)/dT 
y  derivando: 

[3Sn(T,P)/ST]p  -  [0Ss(T,P)/0T]P  -  [aSs(T,P)/eP]T.d(Hc2/8Tí)/dT  = 

=  -[aVs(T,P)/5T]P.d(Hc2/8Tr)/dT  -  [SVs(T,P)/dP]T.[d(Hc2/8TT)/dT]2  -  Vs(T,P).d2(Hc2/8Tr)/dT2. 
Usando  la  relación  de  Maxwell  (3S/3P)T  =  -(3V/c>T)P 

(5Sn/5T)P  -  (aSs/5T)P  =  -  Vs.d2(Hc2/8TT)/dT2  -  2(aVs/5T)P.d(Hc2/8Tr)/dT  -  (5Vs/5P).[d(Hc2/8TT)/dT]2. 


En  general,  las  derivadas  dVs/5T  y  3VS/3P  pueden  despreciarse,  además  a  T=0  y  T=TC  sus 
coeficientes  se  anulan.  Entonces: 

(5Sn/5T)P  -  (3Ss/dT)P  «  -  Vs.d2(Hc2/8TT)/dT2 

o  bien: 

CP,s-Cp,n  =  (TVs/4TT).[(dHc/dT)2+  H0(d2Hc/dT2)]. 

El  primer  término  es  positivo  y  el  segundo  negativo,  en  el  intervalo  0  -  Tc. 

Para  T=TC.  Hc=0  y  resulta  la  expresión  debida  a  Rutgers  (1934): 

CP,s-Cp,n  =  (TVs/4TT)(dHc/dT)To2 

Habrá  entonces  una  discontinuidad  finita  en  la  capacidad  calorífica  como  corresponde  a  una 
transición  estrictamente  de  segundo  orden. 

A  T=TC  la  diferencia  se  achica  y  pasa  a  ser  negativa  pues  prevalece  el  segundo  término. 

De  manera  análoga  pueden  estudiarse  los  efectos  ligados  al  volumen.  Derivando  respecto 
de  P  la  expresión  exacta  para  la  energía  de  Glbbs: 


Vn- Vs  =  (VsHc/4Tr).dHc/dP. 


A  T=TC,  Hc=0  y  entonces  Vn  =  Vs,  desapareciendo  la  discontinuidad  en  el  volumen,  pero  con 
una  discontinuidad  finita  en  la  compresibilidad. 
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CAPÍTULO  13 
Sistemas  binarios 


En  este  capítulo  examinaremos  detalladamente  las  propiedades  asociadas  a  la  composición 
de  los  sistemas,  analizando  el  caso  más  sencillo  de  un  sistema  con  dos  clases  de  partículas. 


Sistemas  con  dos  clases  de  partículas 

Nos  interesa  estudiar  estas  mezclas  binarlas  en  particular  porque  pueden  presentar 
diversas  transiciones  de  fase  ligadas  a  su  composición.  Consideremos  un  sistema  compuesto 
por  Na  moles  de  la  especie  A  y  NB  moles  de  la  especie  B  tal  que  en  la  representación 
energética  sus  variables  Independientes  serán  S,V,NA,NB  y  en  la  representación  de  Glbbs,  que 
usaremos,  T,P,NA,NB.  En  particular  tendremos  como  variables  intensivas  conjugadas  de  los 
números  de  moles,  los  potenciales  químicos: 
pA- (c,G/3Na)t,p,nb  y  |íb  -  (c,G/5Nb)t,p,na 
y  una  relación  de  Maxwell  entre  ellas: 

(3|ía/5Nb)na  =  32G/3Na3Nb  =  (3|íb/3Na)nb 

De  acuerdo  a  lo  visto,  la  estabilidad  intrínseca  de  un  tal  sistema  queda  asegurada  por  tres 
criterios: 

(3T/3S)v,na,nb  >  0  - >  Cy  >  0, 

-(3P/3V)t,na,nb  >  0  — : >  Kj  >  0 

(3|ía/3Na)t  P  NB  =  52G/3Na2  >  0. 

La  ecuación  fundamental  puede  escribirse: 


G  =  G(T,P,Na,Nb)  =  Ng(T,P,xA,xB)  =  Ng(T,P,x); 
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con  N  =  Na+Nb  y  donde  llamamos  x  =  xB  =  NB/N  y  consecuentemente  xA  =  NA/N  =  1-x.  Es 
habitual  nombrar  esta  mezcla  como  A-^Bx.  De  ahora  en  adelante  nos  Interesaremos 
principalmente  por  la  dependencia  con  la  composición,  analizando  las  propiedades  de  g  como 
función  de  x,  suponiendo  constantes  T,P  y  N. 


Propiedades  de  g(T,P,x)9 

Usando  la  ecuación  de  Euler  resulta: 

G=U-TS+PV=  pANA+|iBNB 
g(x)  =  (1  -x)|ía(x)+X|íb(x) 

y  de  la  ecuación  de  Glbbs-Duhem: 

-SdT+VdP+NAdpA+NBdpB=  0 

obtenemos  para  T  y  P  constantes: 

(1-x)dpA+  xdpB=  0 

Derivando  g(x)  a  T,P,N  constantes 

dg  I8x  =  -pA(x)  +  (1-x)5|ía/5x  +  pB(x)  +  xc>|i  B/5x 

y  teniendo  en  cuenta  la  anterior: 

dg Idx  =  |iB(x)  -  pA(x). 

Esta  diferencia  sería  entonces  la  variable  Intensiva  conjugada  de  x  en  la  representación  de 
Glbbs.  De  las  expresiones  para  g  y  su  derivada  obtenemos  los  potenciales  químicos: 

pA(x)  =  g  -  x.3g/c>x 
pB(x)  =  g  +  (1-x).0g/9x 


En  lo  sucesivo  omitiremos  la  dependencia  con  T  y  P,  a  menos  que  sea  necesario  explicitarla. 

98 


Ha'''/ 


MbM 


x 


O 

A 


x 


1 

B 


Figura  13.1 


relaciones  ¡lustradas  en  la  figura  1.  Allí  puede  verse  que  la  tangente  a  g(x)  en  cualquier  punto  x 
corta  a  los  ejes  (x=0)  y  (x=1)  en  los  valores  pA(x)  y  pB(x)  respectivamente. 

Derivando,  a  su  vez,  los  potenciales  químicos,  obtenemos: 

5|ía/3x  =  -x.5g2/3x2  =  -x.gxx 
5|íb/3x  =  (1-x).3g2/3x2  =  (1-x).gxx, 

que  como  veremos  tienen  un  signo  definido. 

En  efecto,  podemos  analizar  cómo  se  traduce  el  nuevo  criterio  de  estabilidad  en  g(x), 
usando  las  relaciones  anteriores  y  teniendo  en  cuenta  que  en  el  cambio  de  variables  de  (NA,NB) 
a  (N,x)  tendremos: 

(5x/3Na)nb  = -x/N  y  (3x/3Nb)na  =  (1-x)/N 

Entonces: 

(32G/3Na2)nb  =  (3|ía/3Na)nb  =  (3|ía/3x)(3x/5Na)nb  =  x2gxx/N  >  0 
de  donde: 


Es  decir  que  g(x)  debe  tener  curvatura  definida  (cóncava  hacia  arriba)  como  requisito  de 
estabilidad  para  la  mezcla  de  composición  x.  Vemos  entonces  que  el  gráfico  de  la  figura  1 
corresponde  a  una  situación  en  que  la  mezcla  homogénea  es  estable  para  todas  las 
composiciones  0<x<1.  Es  el  caso  de  las  mezclas  gaseosas  (homogéneas  en  todo  el  rango  de 
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composiciones)  y  de  la  mayoría  de  las  mezclas  líquidas,  salvo  algunos  caso  de  inmiscibilidad. 
De  todos  modos,  durante  la  condensación,  transición  gas  líquido,  de  una  mezcla  puede  tener 
lugar  una  separación  de  fases,  como  resultado  de  la  competencia  entre  ambas  fases 
homogéneas.  Por  su  parte,  en  mezclas  sólidas  es  habitual  tener  una  solubilidad  limitada,  de  A 
en  B  y  de  B  en  A,  produciéndose  la  separación  en  dos  fases  que  además  de  diferir  en 
composición  pueden  diferir  en  estructura  cristalina. 


Separación  de  fases 

Imaginemos  una  situación  en  la  que  g(x)  tenga  dos  mínimos  con  un  máximo  entre  ellos, 
determinando  una  región  de  inestabilidad  (gxx  <  0)  entre  x'“  y  x'p  (figura  13.2).  Recordando 
nuestro  análisis  del  capítulo  8,  sobre  las  inestabilidades  en  las  curvas  u(s),  concluimos  que  en 
estas  condiciones  al  sistema  le  conviene  separarse  en  dos  fases  de  distinta  composición 
coexistiendo  en  equilibrio,  es  decir,  con  iguales  valores  de  T,  P  y  de  los  potenciales  químicos 

pa  y  pb- 

La  separación  de  fases  es  termodinámicamente  favorable  toda  vez  que  para  una  dada 
composición  x  puedan  encontrarse  estados  adyacentes,  de  composición  entre  0  y  x  y  entre  x  y 
1  respectivamente,  cuya  combinación  en  proporciones  adecuadas  resulte  en  una  energía  libre 
menor  (estabilidad  global).  Ello  ocurre  para  todos  los  sistemas  de  composición  comprendida 
entre  xa  y  xp,  valores  determinados  por  la  tangente  común,  que  se  separarán  en  dos  fases  de 
composiciones  x“  y  xp,  en  proporciones  que  enseguida  determinaremos.  En  virtud  de  las 
expresiones  para  los  potenciales  químicos  y  de  su  interpretación  gráfica  esta  construcción  de  la 
tangente  común  asegura  además  la  Igualdad  de  los  mismos  entre  las  fases  coexistentes: 

Pa(x°)  =  pA(xP)  y  Pb(xq)  =  pB(xP). 

Notemos  que  las  fases  en  que  se  separa  el  sistema  difieren  en  su  composición, 
determinando  un  salto  Ax=  xp-xa  entre  ellas.  Esta  discontinldad,  típica  de  una  transición  de  fase 
de  primer  orden,  juega  aquí  el  mismo  papel  que  As  y  Av  en  las  transiciones  ya  discutidas. 
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HB“=flBP 


Figura  13.2 

Hallemos  ahora  la  fracción  de  moles  (o  de  átomos)  de  cualquier  clase  (A  o  B)  en  cada  fase, 
ca  y  cp  =  1-ca,  para  una  mezcla  de  composición  x: 

a  a .  3  3  aa.B  a  B 

x  =  cx+cx  =  cx+x-cx 

de  donde  surgen  las  fracciones: 

ca  =  (xp-x)/(xp-xa)  en  la  fase  de  composición  xa  (fase  a) 
cp  =  1-ca  =  (x-xa)/(xp-xa)  en  la  fase  de  composición  xp  (fase  (3). 

Relaciones  que  constituyen  la  llamada  regla  de  la  palanca.  Como  aplicación  se  puede  verificar 
que  con  estos  valores  resulta  ca.g(xa)  +  cp.g(xp)  <  g(x)  para  cualquier  composición  en  el 
intervalo  xa<x<xp. 

A  su  vez  en  los  intervalos  0<x<xa  y  <xp<x<1  la  mezcla  será  homogénea  presentándose 
como  fase  a  o  fase  (3,  respectivamente.  Resulta  así  dividido  el  espacio  de  composiciones  en 
tres  regiones  con  estados  de  equilibrio  homogéneos  en  dos  de  ellas  y  heterogénos  en  la 
restante  (figura  13.3  arriba).  En  general,  la  forma  detallada  de  g(x)  Irá  variando  con  T  y  P, 
supuestos  constantes  hasta  ahora,  determinando  una  dependencia  de  los  límites  de  dichas 
reglones  con  las  variables  T  y  P.  Esta  suele  ilustrarse  en  diagramas  de  fase  bldimensionales  en 
los  planos  (T,x)  o  (P,x)  que  muestran  mediante  curvas  T(x)  a  P=cte  o  P(x)  a  T=cte  los  límites 
entre  las  reglones  de  1  y  2  fases  (figura  13.3  abajo). 
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estado  heterogéneo: 
fases  a  y  p 


A 


X 


xP 


B 


Figura  13.3 

Este  análisis  puede  extenderse  a  situaciones  donde  en  lugar  de  tener  una  curva  g(x)  con 
zonas  de  inestabilidad,  simplemente  tengamos  dos  fases  homogéneas  posibles  ga(x)  y  gp(x) 
compitiendo  en  estabilidad.  La  tangente  común  (envolvente  inferior)  determinará  las 
composiciones  x“  y  xp  que  dividen  el  espacio  de  la  forma  mencionada. 

De  esta  manera  pueden  entenderse  fácilmente  los  diagramas  de  fase  de  mezclas  binarias 
con  regiones  en  forma  de  cigarro  para  la  coexistencia  de  fases  gaseosa  y  líquida,  donde 
compiten  las  respectivas  energías  libres  g9as(x)  y  gliq(x).  Al  enfriar  una  mezcla  gaseosa,  el 
primer  condensado  (fase  líquida)  tendrá  en  general  una  composición  diferente  a  la  composición 
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global  del  sistema,  alterando  también  la  composición  de  la  fase  gaseosa  remanente  y  dando 
lugar  aunAx  =  X|¡q  -  xgas  entre  las  fases  coexlstentes.  Asimismo,  puede  describirse  así  el  caso 
de  mezclas  sólidas  con  solubilidad  limitada  de  B  en  A  (fase  a)  y  de  A  en  B  (fase  (3).  En  la  zona 
de  composiciones  intermedias,  entre  los  respectivos  límites  de  solubilidad,  el  sistema  se 
separará  en  dos  fases  que  coexisten  en  equilibrio. 

Recordemos  que  según  la  regla  de  las  fases  de  Gibbs  el  número  de  fases  que  pueden 
coexistir  en  equilibrio  es  como  máximo  r+2,  es  decir  4  en  este  caso.  Esta  regla  también  puede 
expresarse  en  términos  del  número  de  grados  de  libertad  v  (número  de  variables  Intensivas 
que  pueden  variarse  libremente)  que  resulta  v  =  4-cp  >  0  en  el  caso  presente  (r=2).  SI  se 
representaran  las  fase  de  equilibrio  en  un  sistema  de  coordenadas  T,P,pA  (o  pB),  el  equilibrio 
entre  dos  fases  quedaría  determinado  por  una  superficie  (T,P  y  pA  deben  ser  los  mismos  en 
ambas  fases).  La  coexistencia  de  tres  fases  se  dará  sobre  líneas  triples  que  resultan  de  la 
Intersección  de  dos  de  tales  superficies.  Análogamente  la  coexistencia  de  cuatro  fases  se  dará 
en  puntos  cuádruples  que  resultan  de  la  intersección  de  dos  de  tales  líneas. 


Solución  ideal 

Analicemos  ahora  algunos  modelos  simples  capaces  de  describir  adecuadamente  los 
comportamientos  observados:  solución  Ideal  y  solución  regular. 

Una  mezcla  o  solución  de  varios  componentes  se  denomina  ideal  si  el  mezclarlos  no 
produce  efectos  energéticos,  solo  entrópicos.  El  caso  más  simple  es  el  de  una  mezcla  de 
gases  Ideales,  es  decir  una  mezcla  de  partículas  no  interactuantes  (ver  Anexo).  En  casos  más 
realistas,  puede  admitirse  que  haya  interacción  entre  las  partículas  pero  que  ésta  no  resulte 
afectada  por  la  mezcla.  En  cualquier  caso,  el  potencial  químico  de  cada  componente  k  resulta: 
pk  =  Pk°(T,P)  +  RT.In(Xk),  donde  pk°(T,P)  es  el  potencial  químico  de  la  sustancia  k  pura  y  el 
término  adicional  da  cuenta  del  efecto  de  la  mezcla. 

Para  una  mezcla  binarla  Ai_xBx  Ideal  entonces: 

g(T,P,x)  =  (1-x)pA(T,P,x)+xpB(T,P,x) 

g(T,P,x)  =  (1-x)pA°(T,P)  +  xpB°(T.P)  +  RT[(1-x).ln(1-x)+x.ln(x)] 

Omitiendo  las  variables  T  y  P,  podemos  escribir  g(x)  =  g°(x)  +  gm(x),  donde 

g°(x)  =  (1-x)pA°+xpB° 

representa  la  parte  energética  (entálplca)  del  sistema  que  corresponde  a  la  suma  pesada  de 
las  energías  libres  molares  de  sus  constituyentes  y 

gm(x)  =  RT[(1-x).ln(1-x)+x.ln(x)]  =  -T.sm(x) 
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describe  la  reducción  de  energía  libre  por  el  aumento  en  la  entropía  que  provoca  la  mezcla. 
sm  =  -R[(1-x).ln(1-x)+x.ln(x)]  =  -T.sm 

Notar  que  siendo  x  y  (1-x)  cantidades  inferiores  a  1 ,  sus  logaritmos  son  negativos  y  también 
lo  será  entonces  gm  que  toma  su  valor  mínimo  para  x=1/2. 

A  partir  de 

Pa(x)  =  Pa°  +  RT.In(l-x)  y 
Pb(x)  =  pB°  +  RT.In(x) 

podemos  estimar  también  dgldx  para  la  solución  Ideal,  usando  (3): 
dgldx  =  |iB°  -  ma°  +  RT.In[x/(1-x)] 

Se  ve  que  esta  derivada  diverge  en  los  extremos  (^-00  para  x=0  y  -^-+oo  para  x=1)  mostrando 
el  rápido  aumento  de  entropía  que  produce  la  dilución  de  Impurezas  en  un  elemento  puro. 
Asimismo  vemos  que  en  x=1/2,  vale  estrictamente  dgldx  =  pB°  -  Pa°- 
Finalmente  para  la  derivada  segunda: 

d2gldx2  =  RT/x(1-x)  >  0 

mostrando  que  la  mezcla  homogénea  es  estable  para  cualquier  composición.  Resultado  obvio 
ya  que  como  sabemos  un  sistema  Ideal  (sin  interacciones  entre  partículas)  no  presenta 
transiciones  de  fase. 

Sugerencia:  graficar  g(x)  para  una  mezcla  ideal  y  comprobar  que  la  misma  resulta  estable 
(gxx>0)  en  todo  el  rango  de  composiciones. 


Solución  regular 


Introducimos  entonces  el  modelo  de  solución  regular  que  intenta  dar  cuenta  de  manera 
sencilla  de  posibles  interacciones  entre  las  partículas  de  distinta  clase  presentes  en  la  mezcla, 
sin  alterar  el  término  entrópico.  Para  ello  proponemos: 

g(x)  =  gideal(x)  +  g'nter(x)  -  g°(x)  +  gm(x)  +  Ax(l-x) 

donde  incluimos  un  término  de  energía  de  interacción  proporcional  a  la  probabilidad  de 
encuentro  de  partículas  de  especies  diferentes:  xAxB  =  x(1-x)  que  será  positivo  o  negativo 
según  que  la  interacción  entre  ellas  sea  repulsiva  o  atractiva.  Tenemos  entonces: 

g(T,P,x)  =  (1-x)|iA°+xpB°+  RT[(1-x).ln(1-x)+x.ln(x)]  +  Ax(l-x). 
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Si  A<0  el  nuevo  término  no  hace  más  que  contribuir  a  la  estabilidad  de  la  mezcla 
homogénea  sumándose  al  efecto  entrópico,  ambos  negativos,  para  aumentar  la  curvatura  de  g. 
En  cambio,  si  A>0,  en  su  competencia  con  el  término  entrópico  podría  dar  lugar  a  la  aparición 
de  inestabilidades  en  la  forma  de  “jorobas”  en  la  curva  g(x).  En  efecto 
dg/dx  =  pB(x)  -  |ía(x)  =  pB°  -  Pa°  +  RT.In[x/(1-x)]  +  A(1-2x) 
y  entonces: 

d2gldx2  =  RT/x(1  -x)  -  2A  >  4RT  -  2A 

A  una  temperatura  dada,  4RT  -  2A  es  el  valor  mínimo  que  puede  tomar  gxx,  es  decir  que 
para  A<2RT,  será  gxx  >  0,  V  x  y  la  mezcla  homogénea  estable.  Pero  cuando  A  >  2RT  aparece 
un  intervalo  de  valores  de  x,  determinado  por  las  raíces  x'ayx'p  de  la  ecuación  x(1-x)  =  RT/2A, 
dentro  del  cual  gxx  <  0,  dando  lugar  a  inestabilidades  y  a  las  consecuentes  transiciones  de  fase 
tal  como  se  discutieron  en  el  punto  anterior.  Quedan  así  determinados,  la  temperatura  crítica  Tc 
=  A/2R  debajo  de  la  cual  puede  haber  separación  de  fases  y  la  curva  espinodal  gxx  =  0  dentro 
de  la  cual  los  estados  homogéneos  son  estrictamente  inestables.  Esta  curva  vendrá  dada  por 
las  raíces  x'a(T)  y  x'p(T)  de  la  ecuación: 

x2  -  x  +  T/4TC=  0 
x  =  (1/2).[1±V(1-T/Tc)j 
simétricas  respecto  de  x=1/2. 

Asimismo,  se  pueden  obtener  para  este  modelo  los  potenciales  químicos  individuales: 

pA(x)  =  pA°  +  RT.In(l-x)  +  Ax2, 
pB(x)  =  pb°  +  RT.In(x)  +  A(1-x)2, 

Sugerencia:  graficar  g(x)  para  una  mezcla  regular  y  analizar  la  aparición  de  inestabilidades  y  la 
separación  de  fases. 

Las  composiciones  xa  y  xp  de  las  fases  que  coexistirán  en  equilibrio  pueden  determinarse 
teniendo  en  cuenta  que  el  valor  común  de  la  tangente  a  la  curva  g(x)  en  esos  puntos  será  la 
pendiente  de  la  recta  envolvente  inferior:  |iB°  -  Pa° 

(5g/ax)a,3  =  |iB°  -  Ma°  +  RT.In[x/(1-x)]  +  A(1-2x)  =  pB°  -  Pa° 

entonces:  RT.In[x/(1-x)j  +  A(1-2x)  =  0 

o  bien:  ln[x/(1-x)]  =  2Tc(2x-1)/T. 

Esta  ecuación  se  satisface  trlvlalmente  en  x=1/2  donde  ambos  miembros  se  anulan  y  el 
estado  es  inestable  ya  que  gxx  =  4RT  -  2A  =  4R(T-TC)  <  0.  Para  T=TC,  gxx  =  0  y  la  estabilidad 
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debe  revisarse  en  órdenes  superiores.  Este  punto,  tiene  las  características  de  un  punto  crítico, 
donde  desaparece  la  distinción  entre  las  fases  que  coexisten,  determinando  una  composición 
crítica  xc=1/2  que  es  el  valor  al  que  van  a  tender  ambas  soluciones  x°  y  xp,  y  las  de  la  espinodal 
x'a  y  x'p,  cuando  T  — >  Tc. 

Esto  sugiere  definir  valores  reducidos  de  la  composición,  que  miden  el  apartamiento  relativo 
de  xa  y  xp  con  respecto  a  xc: 

X  =  (x-xc)/xc  =  (x-y2)/y2  =  2x-1. 

En  términos  de  x,  la  ecuación  anterior  resulta: 

ln[(1+x)/(1-X)]  =  2TcX/T 

ecuación  trascendente  que  puede  reescribirse: 

(1+X).exp(-TcX/T)  =  (1-X).exp(TcX/T) 

o  bien: 

X  =  tgh(Tcx/T). 

Esta  ecuación  trascendente,  similar  a  la  de  Curie-Welss  para  la  magnetización,  puede 
resolverse  gráficamente.  Como  sabemos,  por  debajo  de  T=TC,  dará  lugar  a  dos  soluciones 
simétricas  Xp  =  ~Xa,  distintas  de  la  trivial  x-0  (x=xc=1/2).  Estas  composiciones,  o  su  diferencia, 
juegan  el  papel  de  parámetro  de  orden  para  esta  transición,  con  un  comportamiento  crítico 
similar  al  predlcho  por  las  teorías  clásicas  o  de  campo  medio  (p=1/2).  En  efecto,  cuando  T->TC, 
X  ->XcY  podemos  aproximar: 

X  *  (TcX/T)  -  (T cX/T)3/3 

X2»  -3(T/Tc)3(1-Tc/T)  =  -36(1  +02) 

X~  ±(-30)1/2 


confirmando  el  valor  mencionado  del  exponente  crítico  p  =  %.. 
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Aplicaciones  del  modelo  de  solución  regular 


En  las  aplicaciones  a  mezclas  reales  fluidas ,  líquidas  o  gaseosas,  suele  definirse  un 
coeficiente  de  actividad  y  que  expresa  el  apartamiento  del  comportamiento  Ideal.  SI  en  una 
mezcla  Ideal:  pB(x)  =  Pb°  +  RT.In(x), 
para  el  caso  real  se  define: 

Pb(x)  =  pB°  +  RT.In(a) 

con  la  actividad  a  =  yx. 

En  el  modelo  de  solución  regular: 

pB(x)  =  |iB°  +  RT.In(x)  +  A(1-x)2. 

de  modo  que  por  comparación: 

RT.In(a/x)  =  A(1-x)2. 

Y  así: 

y  =  a/x=  exp[A(1-x)2/RT] 

y  expresiones  análogas  para  el  componente  A. 

Podemos  considerar  dos  casos  extremos: 

i)  x-»1  (sistema  B  casi  puro) 

y— >1;  a— »x 

comportamiento  Ideal  que  se  describe  empíricamente  por  la  ley  de  Raoult. 

ii)  x-^0  (B  muy  diluido  en  A) 

y— >  exp(A/RT)  =  cte  (<1  si  A>0;  >1  si  A<0) 
comportamiento  lineal  descripto  por  la  ley  empírica  de  Henry. 

En  su  aplicación  a  mezclas  sólidas  (solución  sólida  regular),  puede  relacionarse  A  con  la 
energía  de  interacción  de  pares  de  partículas  distintas,  e.  En  primer  lugar  debemos  tener  en 
cuenta  que  en  un  sólido  los  átomos  están  fijos  en  sitios  de  un  arreglo  periódico  (red  cristalina) 
cada  uno  de  ellos  rodeado  de  un  número  Z  de  vecinos  más  próximos.  Estos  podrán  ser  de  la 
misma  especie  o  no,  dando  lugar  a  un  cambio  en  la  energía  por  cada  par  de  átomos  distintos, 
que  puede  expresarse  como 

£  =  £ab_(£aa+£bbV2, 

por  comparación  con  la  situación  promedio  en  las  sustancias  A  y  B  puras. 
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Para  hallar  el  cambio  total  en  la  energía  debemos  determinar  el  número  de  pares  de  vecinos 
próximos  distintos,  NAB.  Considerando  un  mol  de  sustancia,  el  número  de  átomos  A  será: 

N0xa.  Cada  uno  tendrá  Z  primeros  vecinos  de  los  cuales,  suponiendo  una  distribución  al  azar, 
una  fracción  xB  será  de  la  especie  opuesta,  formando  pares  AB  en  número: 

Nab  =  N0xaxbZ  =  N0Zx(1-x) 

De  modo  que  la  contribución  energética  a  la  energía  de  Glbbs  por  mol  será: 
ginter(x)  =  NABe  =  N0Z£.x(1-x) 

lo  que  permite  Interpretar  la  constante  A  del  modelo,  como  N0Ze. 


Transiciones  orden-desorden 

En  las  soluciones  sólidas,  además  de  las  transiciones  de  fase  ligadas  a  la  composición 
como  las  ya  descrlptas,  pueden  aparecer  transiciones  desorden-orden  entre  una  solución 
sólida  al  azar  y  un  compuesto  donde  los  átomos  diferentes  ocupan  sitios  determinados  en  la 
estructura.  Esto  Implica  dar  un  paso  más  en  la  descripción  de  las  soluciones  sólidas  a  través  de 
g(T,P,x).  Al  término  de  sustancias  puras  (yuxtaposlclónde  A  y  B  en  las  proporciones  dadas) 
hemos  agregado  primero  un  término  entrópico  originado  por  la  simple  mezcla  de  sustancias 
distintas  y  luego  un  efecto  energético  por  el  cambio  en  la  interacción  entre  partículas  que 
Implica  esta  mezcla.  Pero  si  esta  interacción  es  suficientemente  intensa,  podrá  provocar 
cambios  en  la  distribución  de  partículas,  que  ya  no  será  al  azar,  modificando  la  entropía  y  la 
energía  total  de  interacción. 

Hay  casos  entonces  en  que,  para  una  composición  fija,  la  Interacción  atractiva  (e<0  )  entre 
átomos  de  distinta  clase  es  suficiente  para  forzar  ordenamientos  atómicos  que  favorecen  la 
formación  de  pares  de  primeros  vecinos  distintos  AB  por  sobre  el  ordenamiento  al  azar 
Implícito  en  la  solución  regular.  En  la  expresión  utilizada  para  g'nter(x)  hemos  supuesto  que  la 
probabilidad  PAB  de  que  un  sitio  vecino  a  un  átomo  A  esté  ocupado  por  un  átomo  B,  es 
justamente  x=xB.  Supondremos  entonces  que  pueden  producirse  cambios  en  dicha 
probabilidad  dando  lugar  a  una  contribución  a  la  energía  libre: 

ginter(x)  -  N0xAPABZe  =  N0Ze(1-x)  PAB 

Por  simplicidad  fijaremos  la  composición  en  x=1/2  y  entonces  la  parte  energética  de  la  energía 
libre  por  la  Interacción  será: 

ginter  =  N0Ze(1-x).PAB  =  A.Pab/2. 

Supongamos  además  que  la  estructura  cristalina  puede  descomponerse  en  dos  super-redes 
¡guales  ¡nterpenetradas,  1  y  2,  de  modo  que  para  un  punto  dado  en  una  de  ellas,  los  Z  vecinos 
próximos  se  encuentran  en  la  otra  y  viceversa.  La  red  cúbica  centrada  en  el  cuerpo  (bcc), 
compuesta  de  dos  redes  cúbicas  simples  (se),  constituye  un  buen  ejemplo.  Macroscópicamente 

el  grado  de  ordenamiento  producido  por  la  interacción  puede  caracterizarse  mediante  el 
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parámetro  de  orden  de  largo  alcance  L,  vinculado  a  las  probabilidades  de  que  los  átomos  de 
tipo  A  y  B  ocupen  correctamente  los  sitios  propios  en  la  estructura:  A  en  la  red  1  y  B  en  la  red  2 
o  viceversa: 

Pai  =Pb2  =  (1+L)/2  Pa2  =  Pbi  =  (1-L)/2 

Para  L=1  tendremos  orden  completo  y  entonces  cada  átomo  A  estará  rodeado  de  Z  átomos  B  y 
viceversa. 

Sugerencia:  analizar  la  situación  para  otros  valores  de  L  entre  -1  y  1. 

Finalmente  podemos  verifica  que 

PAB  =  Pai-Pb2  +  PA2-Pbi  -  (1+L2)/2 

y  entonces 


ginter(L)  =  A(1+L2)/4 


que  tiende  al  valor  de  solución  regular,  A/4,  para  L=0  como  se  espera. 

A  su  vez,  estos  ordenamientos  van  a  modificar  la  entropía  del  sistema  con  respecto  a  la  de 
la  mezcla  al  azar.  Mediante  métodos  de  la  Mecánica  Estadística  se  puede  establecer  que  la 
entropía  molar  será  en  este  caso: 


s(L)  =  -N0k[PA1.ln(PA1)  +  PBi.ln(PBi)  +  PA2.ln(PA2)  +  PB2.ln(PB2)]/2, 


s(L)  =  -R{(1  +L)ln[(1  +L)/2]+(1  -L)ln[(1  -L)/2]}/2, 

que  tiende  a  sm  (mezcla  Ideal)  para  L=0  como  se  espera.  Finalmente  tendremos: 
g(L)  =  g°  +  A.(1+L2)/4  +  RT.[(1+L).ln(1+L)+(1-L).ln(1-L)]/2. 

L  tomará  valores  que  minimicen  g,  es  decir,  que  satisfagan: 

dgldl  =  A.L/2  +  RT.In[(1+L)/(1-L)]/2  =  0 

Expresión  que  conduce  a  una  ecuación  trascendente  para  L: 

L  =  -  tgh(AL/2RT) 

del  tipo  de  las  que  ya  hemos  analizado  y  que  conduce  a  soluciones  no  triviales  (L^O)  para  T  < 
Tc  =  -Ze/2k  =  -A/2R  y  un  comportamiento  crítico  similar,  propio  de  teorías  de  campo  medio. 
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Anexo:  Entropía  de  mezcla 


Supongamos  una  mezcla  de  r  gases  Ideales  distintos  ocupando  un  volumen  V  a  la 
temperatura  T.  Si  hay  N¡  moles  de  cada  uno  de  ellos,  la  presión  resultante  será  P  =  NRT/V  con 
N  =  IN¡.  Queremos  determinar  la  entropía  de  mezcla,  es  decir,  aquella  parte  de  su  entropía  que 
puede  atribuirse  estrictamente  al  proceso  irreversible  de  haberlos  mezclado.  Con  ese  fin 
compararemos  la  entropía  S  del  estado  mezclado  con  la  de  un  estado  de  referencia  S°  donde 
los  r  gases  se  encuentren  a  la  misma  temperatura  y  presión,  pero  separados  en  recipientes  de 
volúmenes  V¡.  Para  que  la  presión  sea  la  misma  debe  cumplirse  V¡/V  =  N¡/N  =  x¡  y 
evidentemente  IV¡  =  V. 

Primeramente,  recordando  la  ecuación  fundamental  de  un  gas  Ideal,  podemos  escribir: 

S  =  Nct(T)  +  NR.In(V/N) 

donde  hemos  separado  la  dependencia  con  U  en  el  término  No(T),  teniendo  en  cuenta  que 
U=NCVT.  Asimismo  para  las  respectivas  energías  libres  resulta: 

F  =  Ncp(T)  -  NRT.In(V/N) 

G  =  Ny(T)  +  NRT.In(P) 

Ahora,  en  virtud  de  la  aditividad  de  la  entropía  podemos  escribir  para  el  estado  de 
referencia: 


S°  =  Z[N¡ct¡(T)  +  N¡R.ln(V¡/N¡)]. 


Para  el  estado  mezclado,  podemos  utilizar  la  misma  expresión,  teniendo  en  cuenta  que  ahora 
cada  gas  ocupa  todo  el  volumen  V,  como  si  los  otros  no  existieran: 


S  =  Z[N¡ct¡(T)  +  N¡R.ln(V/N¡)]. 


Resulta  así  una  entropía  de  mezcla: 

Sm  =  S  -  S°  =  ZN¡R.[ln(V/N¡)  -  ln(V¡/N¡)]  =  R.ZN¡.ln(V/V¡) 
que  puede  escribirse: 

Sm  =  NR.Z(N¡/N)ln(N/N¡) 
o  bien: 


Sm  =  -NR.Ix¡.ln(x¡) 
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Veamos  qué  resulta  para  la  energía  libre  de  Glbbs  G  =  U  +  PV  -  TS  =  Ip¡N¡  como 
consecuencia  de  esta  entropía  adicional.  Tanto  la  energía  como  el  volumen  son  aditivos  sobre 
los  sistemas  separados  para  dar  el  sistema  mezclado  a  la  misma  presión,  de  modo  que  la 
diferencia  entre  G  y  G°  sólo  vendrá  del  término  entrópico.  Podemos  escribir  entonces: 

G  =  G°  -  TSm  =  IN¡.[y¡(T)  +  RT.In(P)]  +  RT.IN¡.ln(x¡) 

G  =  IN¡.[y¡(T)  +  RT.In(P)  +  RT.In(x¡)] 

que,  introduciendo  las  presiones  parciales,  P¡  =  x¡P,  suele  escribirse: 

G  =  IN¡.[y¡(T)  +  RT.In(P¡)]. 

Finalmente  los  potenciales  químicos  de  cada  especie  serán: 

|i¡  =  y¡(T)  +  RT.In(P)  +  RT.In(x¡) 

donde  y¡(T)  +  RT.In(P)  representa  el  potencial  químico  de  la  sustancia  pura  que  llamaremos 
|i¡°(T,P)  y  entonces: 

W  =  Mí°(T,P)  +  RT.In(x¡) 

para  cada  componente  en  la  mezcla. 


111 


CAPÍTULO  14 
Equilibrio  químico 


En  este  capítulo  discutiremos  cómo  incluir  en  nuestra  descripción  la  posible  existencia  de 
reacciones  químicas  y  otras  transformaciones  internas  en  el  sistema. 


Reacciones  químicas 


Consideremos  ahora  cómo  extender  nuestro  tratamiento  termodinámico  a  sistemas  entre 
cuyos  constituyentes  pueden  ocurrir  reacciones  químicas.  Una  reacción  química  es  un  proceso 
en  el  que  algunos  constituyentes  desaparecen  dando  lugar  a  la  aparición  de  otros.  Como 
resultado,  los  diversos  números  de  moles  (variables  termodinámicas)  cambian,  aumentando 
unos  a  expensas  de  otros.  Las  relaciones  entre  los  números  de  moles  cambiantes  están 
regidas  por  las  leyes  básicas  de  la  química  y  se  describen  por  medio  de  ecuaciones  químicas 
tales  como: 

2H2  +  02  <-»  2H20 
02  <->  20 
CI5P  <-»  CI3P  +  Cl2 
2Fe203  +3T¡  <-»  4Fe  +  3Ti02 

El  significado  de  estas  ecuaciones,  por  ejemplo  la  primera,  es  que  los  cambios  en  los 
números  de  moles  de  H2,  02  y  H20  están  en  la  relación  -2:-1:2.  La  elección  de  signos  es 
arbitraria;  obviamente  tienen  que  ser  diferentes  a  ambos  lados  de  la  ecuación  para  indicar  que 
la  disminución  de  unos  ( reactivos )  implica  el  aumento  de  otros  ( productos ).  Los  números  de 
moles  de  cada  especie  química  o  componente  pueden  ser  cualesquiera  y  no  satisfarán  en 
general  una  relación  sencilla,  pero  sus  cambios  sí,  a  saber: 
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-dNh¡drogenc/2  -dNox¡geno/1  dNagUa/2  ÓC 


donde  dE,  es  un  factor  de  proporcionalidad  denominado  avance  de  la  reacción. 

Para  una  reacción  química  general,  que  involucre  a  las  especies  A1,A2,...,Ak  como  reactivos 
y  Ak+1,Ak+2,  ...,An  como  productos,  podemos  escribir: 


-  v-|Ai  -  v2A2  -  ...  -  vkAk  <-»  vk+iAk+i  +  vk+2Ak+2  +  ...  +  vnAn 
2r-v¡A¡  <-»Xpv¡A¡ 

donde  los  v¡  son  los  llamados  coeficientes  estequiométricos  y  los  subíndices  R  y  P  indican 
sumas  sobre  los  reactivos  y  productos  respectivamente.  Del  mismo  modo  que  en  el  ejemplo, 
esta  ecuación  significa  que: 

dN^v-i  =  dN2/v2  =  ...  =  dNn/vn  =  di; 

relacionando  las  variaciones  en  los  números  de  moles  de  las  especies  participantes  en  la 
reacción.  Es  conveniente  definir  coeficientes  estequiométricos  que  se  extiendan  a  todos  los 
constituyentes  del  sistema,  dándole  valor  nulo  a  los  de  aquellos  que  no  participen  de  la 
reacción  y  escribir  la  ecuación10  química  de  manera  compacta: 


£v¡A¡  <-»  0 


o  bien: 


dN¡  =  v¡.d£,  i  =  1,r 


La  identificación  de  productos  y  reactivos  y  con  ello  la  atribución  de  signos  a  los  coeficientes 
estequiométricos  es  arbitraria  ya  que  la  reacción  puede  ocurrir  en  ambas  direcciones.  Un 
criterio  bastante  usado  es  el  de  considerar  positiva  la  dirección  de  ocurrencia  en  que  la 
reacción  libera  calor  ( exotérmica )  y  negativa  la  dirección  contraria  ( endotérmica ).  Notemos  que 
el  signo  de  d^  define  el  sentido  de  avance  de  la  reacción. 

Las  consideraciones  presentes  pueden  extenderse  a  casos  en  que  en  el  sistema  ocurran 
otros  procesos  de  transformación  interna  como  ionización,  formación  de  pares,  reacciones 
nucleares  u  otras  reacciones  generalizadas  en  las  que  el  número  de  constituyentes  varíe  según 
una  dada  ecuación  “química”  característica. 


10  Si  bien  la  (14.1)  es  una  ecuación  simbólica,  ella  está  expresando  la  ecuación  de  conservación  de  la  masa.  En  efecto: 
dm  =  ZM¡dN¡  =  Xv¡M¡d5  =  0  implica  Xv¡M¡  =  0. 
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Termodinámica  del  equilibrio  químico 


Se  trata  ahora  de  revisar  las  condiciones  de  equilibrio  en  presencia  de  la  reacción  química. 
En  el  estado  de  equilibrio  las  variables  termodinámicas,  en  particular  los  N¡,  no  cambian  con  el 
tiempo.  Ello  no  implica  que  la  reacción  no  ocurra,  sino  mas  bien  que  está  ocurriendo  en  ambos 
sentidos,  con  igual  velocidad,  dando  lugar  a  concentraciones  constantes  (dN¡  =  0)  desde  el 
punto  de  vista  macroscópico.  En  un  sistema  a  temperatura  T  y  presión  P  constantes,  los 
estados  de  equilibrio  son  los  que  minimizan  la  energía  libre  de  Gibbs  G(T,P,N-i,N2,...,Nr). 
Consideremos  un  apartamiento  virtual  del  equilibrio  di;,  con  lo  que  los  N¡  variarán  según  2, 
dando  siempre  dG=0  que  es  el  requerimiento  para  el  equilibrio  en  las  condiciones  supuestas: 

dG  =  X|i¡dN¡  =  Zp¡v¡.dí;  =  0,  V  d^. 

Entonces  la  condición  de  equilibrio  químico  resulta: 

Z  v¡|i¡  =  0. 

que  es  una  ecuación  para  la  composición  del  sistema  ya  que  p¡  =  p¡  (T,P,Ni,N2,...,Nr)  y  los  v¡ 
son  constantes.  En  general  para  determinar  completamente  la  composición  será  necesario 
conocer  la  composición  inicial  {N¡0}  y  el  grado  de  avance  Para  cada  componente  tendremos 
en  el  estado  de  equilibrio:  N¡  =  N¡°  +  v£.  y  entonces 


p¡  =  p¡(T,P,^;comp.ini.)  y 


Z  v¡jLi¡  =  Z  v¡|i¡(T,P,i;;comp.ini.)  =  0 


ecuación  que  permite  determinar  i;  y  consecuentemente  la  composición  {N ¡}  para  las  dadas  T  y 

P. 

Al  relacionar  entre  sí  los  potenciales  químicos,  esta  ecuación  restringe  en  uno  los  grados  de 
libertad  del  sistema.  Si  el  sistema  consiste  de  varias  (cp)  fases  en  equilibrio,  como  todos  los  p¡ 
deben  ser  iguales  entre  las  fases,  la  condición  de  equilibrio  químico  se  cumplirá  en  todas  ellas. 
En  caso  de  que  en  el  sistema  ocurran  varias  (q)  reacciones,  cada  una  caracterizada  por  sus 
coeficientes  estequiométricos  {v¡j },  deberán  cumplirse  simultáneamente  ecuaciones  de  ese  tipo 
para  cada  reacción.  De  modo  que  el  número  de  grados  de  libertad  del  sistema  resulta: 
f  =  r+2-cp-q  y  la  regla  de  las  fases  de  Gibbs  puede  reescribirse: 


cp  <  r+2-q. 


Para  adecuar  nuestra  formulación  termodinámica  al  caso  presente  es  conveniente  definir 
una  nueva  magnitud,  intensiva,  la  afinidad  química, 

A  =  -Z  v¡p¡. 
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que,  según  la  expresión  para  dG,  resultará  la  variable  conjugada  del  grado  de  avance 
(extensiva)  en  la  representación  de  Glbbs: 

A  =  -(dGIdtyjp 

dG  =  -SdT  +  VdP  -Ad^ 

En  otras  representaciones  resultarán  expresiones  análogas,  por  ejemplo: 

dU  =  TdS  -  PdV  -Adq, 

dS  =  (1/T)dU  +  (P/T)dV  +  (AIT)d^. 

La  condición  de  equilibrio  adquiere  entonces  la  forma  A=0,  diferente  a  la  forma  habitual  que 
expresa  la  Igualdad  de  fuerzas  (presiones,  potenciales  químicos,  etc)  entre  las  partes.  SI 
consideramos  al  sistema  conceptualmente  descompuesto  en  dos  subsistemas  conteniendo 
respectivamente  los  reactivos  y  los  productos,  el  estado  de  equilibrio,  que  resulta  del  balance 
exacto  entre  las  fuerzas  que  tienden  a  pasar  partículas  de  uno  al  otro,  puede  describirse  por  la 
condición:  Zr-v¡|i¡  =  Zpv¡|i¡,  que  recuperaría  así  la  forma  habitual:  AR  =  AP. 

Si  bien  sólo  un  valor  determinado  de  la  afinidad,  A=0,  se  asocia  con  el  estado  de  equilibrio  a 
las  dadas  T  y  P,  podemos  extender  su  definición  a  otros  estados  vecinos  (T,P,x1,...,xM) 
conectados  con  el  de  equilibrio  por  la  reacción,  dN¡  =  v¡.d^,  y  que  serían  estados  de  equilibrio 
estables  si  suprimiéramos  la  reacción.  La  afinidad  resultante  para  esos  estados  será  no  nula 
indicando  el  apartamiento  del  equilibrio  y  la  consecuente  ocurrencia  de  la  reacción  en  la 
dirección  que  tienda  a  anular  A  (ver  figura  14.1). 


Estabilidad 

De  acuerdo  con  la  forma  general  de  los  criterios  de  estabilidad,  en  presencia  de  una 
reacción  química  deberá  cumplirse: 

G«  =  (52G/5^2)tp  =  -(cA/S^tp  >  0, 

además  de  los  criterios  habituales,  Cv>0  y  kt>0,  que  en  este  caso  deberán  tomarse  a  ^=cte.  En 
el  gráfico  de  la  figura  14.1  puede  hacerse  una  interpretación  directa  de  este  criterio. 

También  deberán  cumplirse  relaciones  de  Maxwell  adicionales,  que  establecen  propiedades 
para  las  derivadas  parciales  de  A: 
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Figura  14.1 


(3A/5P)T£  =  -(c>2G/0Pc^)t  =  -(oV/oc)TP, 

(oA/r3T)Pj  =  -(52G/5Ta^)P  =  (5S/5^)TP. 

Para  dar  cuenta  de  los  efectos  caloríficos  que  acompañan  a  una  reacción  suele  definirse  el 
calor  de  reacción  como  el  calor  entregado  al  sistema,  a  T  y  P  constantes,  por  grado  de  avance 
unitario: 

Qr  =  (dQ/d^)TP  =  T  (3S/^)tp  =  T  (oA/oT)P, 

donde  hemos  usado  una  de  la  relaciones  de  Maxwell.  Las  reacciones  que  liberan  calor  (Qr<0) 
se  denominan  exotérmicas  y  las  que  lo  absorben  (Qr>0)  endotérmicas.  Este  calor  está 
asociado  con  la  diferencia  en  entalpia11  entre  productos  y  reactivos  AH  y  cambia  de  signo  al 
cambiar  el  sentido  de  la  reacción.  Por  tanto  si  una  reacción  es  exotérmica  en  un  sentido  (d^>0, 
por  convención)  será  endotérmica  en  el  sentido  contrario  (d^<0)  y  viceversa. 

Es  interesante  analizar,  a  partir  del  criterio  de  estabilidad,  el  comportamiento  del  sistema 
cuando  se  intenta  modificar  su  temperatura  o  presión.  Analicemos  cómo  se  desplaza  el 
equilibrio  cuando  queremos  aumentar  la  temperatura  en  dT,  manteniendo  constantes  la  presión 
y  la  afinidad  (A=0): 


1 1  Podríamos  haber  definido  también:  Qr  =  (dH/d^)P  =  T  (dSI8QP  -  A  =  T  (5A/5T)P  pues  A=0,  usando  dH=TdS+VdP- 
Ad^. 
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di;  =  (8ydT)PA  dT 


(3^/3T)ap  =  -(dNoX)¡pl(dNot)jP  =  (Qr/T)/Gjj 

Por  estabilidad,  el  denominador  debe  ser  positivo  y  por  lo  tanto  el  signo  de  la  derivada  será 
el  mismo  que  el  del  calor  de  reacción.  Así,  las  reacciones  endotérmicas  avanzarán  (d^>0)  y  las 
exotérmicas  retrocederán  (d^<0),  tomando  siempre  calor  del  medio  y  disminuyendo  la  acción 
exterior  que  intentaba  aumentar  la  temperatura. 

Análogamente,  cuando  queremos  aumentar  la  presión,  manteniendo  constantes  la 
temperatura  y  la  afinidad  (A=0),  el  grado  de  avance  cambiará  en: 

di;  =  (3^/<9P)at  dP. 

(3^/3P)at=  -(3A/3P),t/(3A/3^)pt  =  -(3V/3^)PT/G^, 

y  entonces  el  signo  de  la  derivada  será  el  opuesto  al  del  cambio  de  volumen  por  grado  de 
avance.  De  manera  que  la  reacción  avanzará  en  el  sentido  en  que  disminuya  el  volumen  del 
sistema,  una  vez  más  reduciendo  los  efectos  de  la  acción  exterior  (compresión).  Estos  efectos 
pueden  ser  muy  Importantes  en  fase  gaseosa  para  reacciones  donde  cambian  los  números 
totales  de  moles  y  por  tanto  los  volúmenes. 


Ley  de  acción  de  masas 

Cuando  el  sistema  en  el  que  la  reacción  química  tiene  lugar,  puede  considerarse  como  una 
mezcla  Ideal,  tendremos: 


m  =  m°(T,P)  +  RT.In(Xi), 


Entonces,  la  condición  de  equilibrio  (4)  puede  expresarse: 
Zvi[|ii°(T,P)+RT.Inx¡]  =  Iviiü0  +  RTZIn(x¡,evi)  =  0 


o  bien: 

nx¡,evi  =  exp(-Iv¡pi°/RT)  =  K(T,P) 

que  es  la  llamada  ley  de  acción  de  masas  (Guldberg-Waage).  Hemos  introducido  la  constante 
de  equilibrio  K  que  define  una  relación  entre  las  concentraciones  de  equilibrio  para  cada  T  y  P. 
La  sumatorla  en  el  argumento  de  la  exponencial  puede  interpretarse  como  un  balance 
energético  de  la  reacción  considerando  los  componentes  puros  (figura  14.2): 

As  =  Zv¡|i¡°  =  Zpv¡|a¡°  -  Zr-v¡|í¡° 
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cuyo  signo  determinará  si  en  la  mezcla  prevalecen  los  reactivos  nRx¡  evl  >nPx¡evl  (para  Ae>0),  o 
los  productos  nPx¡ievl  >nRx¡  evl  (para  Ae<0). 

Invirtiendo  la  expresión  para  K(T,P): 

-EvíiV5  =  RT  InK 

y  la  afinidad  puede  entonces  reescribirse: 

A  =  -Z  v¡p¡  =  -Zvin0,  -  RTSIn(x¡vi)  =  RT.In(K/nx¡vi). 

Como  se  discutió  antes,  esta  definición  de  A  puede  considerarse  aplicable  no  sólo  al  estado  de 
equilibrio  sino  también  a  otros  estados  conectados  con  él.  Para  el  caso  de  equilibrio  (A=0)  la 
(15)  se  reduce  a  la  ley  de  acción  de  masas  que  es  la  que  define  el  valor  de  K(T,P)  cuando  se 
usan  las  concentraciones  x¡i6  de  equilibrio. 

A 
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Figura  14.2 

Veamos  cómo  estimar  el  calor  de  reacción  a  partir  de  esta  expresión  para  A: 

Qr  =  T(oA/oT)Pc  =  RT.In(K/nx¡vi)  +  RT2.(5lnK/5T)Pi=  =  A  +  RT2.(dlnK/r7T)P 
y  teniendo  en  cuenta  que  en  equilibrio  A=0,  resulta  la  conocida  relación  de  van’t  Hoff: 

(5lnK/5T)P  =  Qr/RT2. 


SI  además  el  sistema  donde  ocurre  la  reacción  puede  tratarse  como  una  mezcla  de  gases 
Ideales: 
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Hi°(T,P)  =  Y¡(T)  +  RT.InP 

lü  =  y¡(T)  +  RT.InP  +  RT.Inx¡  =  y¡(T)  +  RT.InP¡. 


Entonces  la  ley  de  acción  de  masas  resulta: 


nP¡vi  =  exp(-Iv¡Y¡/RT)  =  K’(T) 


o  bien 

nx¡vi  =  P'Zv¡  K'(T)  =  K(T,P). 

Como  en  el  caso  anterior,  podemos  escribir: 

A  =  -Zv¡Yi  -  RTZIn(P¡vi)  =  RT  ln(K’(T)  /nP¡vi ) 

Cinética 

En  el  estudio  de  la  termodinámica  de  procesos  irreversibles  las  reacciones  químicas 
proveen  buenos  ejemplos  de  sistemas  fuera  de  equilibrio.  Por  eso  es  conveniente  hacer  un 
análisis,  aunque  sea  sencillo,  de  la  cinética  de  las  mismas  y  de  su  relación  con  las  magnitudes 
termodinámicas  que  hemos  definido. 

Como  hemos  dicho,  el  equilibrio  puede  pensarse  como  resultado  de  un  balance  exacto 
entre  la  reacción  directa  R^P  que  ocurre  con  velocidad  vD  y  la  inversa  P^R  que  ocurre  con 
velocidad  V|.  La  velocidad  neta  v=vD-V|  se  anulará  en  equilibrio.  La  idea  más  simple  que 
podemos  hacernos  de  una  reacción  química  elemental  es  la  de  una  colisión  inelástica  en  la  que 
entran  los  reactivos  (productos)  y  salen  los  productos  (reactivos),  liberando  (consumiendo)  una 
cierta  cantidad  de  energía  (Qr).  Podemos  pues  escribir  las  velocidades  de  reacción  como 
proporcionales  a  las  probabilidades  de  encuentro,  es  decir: 

VD  =  kD-X1  V  .X2  ....Xk  V  =  koLlR  x¡ivl 

Vl  =  kl.xk+1vk+1.xk+2vk+2....xrvr=klnp  x¡vi 

donde  kD  y  k¡  son  las  llamadas  constantes  cinéticas.  Estas  dependerán  de  las  barreras  de 
potencial  que  deban  superarse  para  dar  lugar  a  las  respectivas  reacciones.  Entonces: 

v  =  d£,/dt  =  vD  (1-v,/vD)  =  vD  [1-(nx¡vi)/(  kD/k,)] 

Pero  en  el  equilibrio  (x¡=  x¡  e),  tanto  v  como  A  se  anulan  y  entonces  resulta: 

K  =  nx¡,evi  =  kD/k, 
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Así: 


v  —  vD  [i-(nx¡'7K)] 

y,  como,  al  invertir  la  expresión  para  la  afinidad: 
ln(nx¡v¡/K)  =  -A/RT, 

resulta: 

v  =  vD  [1-exp(-A/RT)] 

Esta  es  la  ecuación  cinética  que  relaciona  la  velocidad  de  reacción  neta,  con  la  afinidad 
(apartamiento  del  equilibrio).  Relación  que  también  puede  escribirse: 

A  =  RT  ln(vD/V|) 

Para  reacciones  no  muy  apartadas  del  equilibrio  A«RT  (2.5  kJ/mol  a  300K)  resulta  una 
relación  lineal: 

v  =  vd.(A/RT)  «  v0A/RT 

donde  v0  es  la  velocidad  de  reacción  directa  en  equilibrio  químico.  En  ese  caso,  entonces,  el 
coeficiente  cinético  (que  definiremos)  será:  L=v0/R. 

En  el  otro  extremo,  cuando  A/RT  oo,  muy  lejos  del  equilibrio: 

v  vDy  V|  «  0. 

Típicamente  es  fácil  alcanzar  valores  de  A  en  el  rango  10-100  kJ/mol,  apartándose  del  régimen 
lineal. 
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CAPÍTULO  15 

Efectos  de  superficie 


En  este  capítulo  extenderemos  la  aplicación  de  nuestro  formalismo  a  sistemas  en  los  que  no 
puede  despreciarse  el  efecto  de  las  superficies  límite  y  de  separación  entre  fases. 


Tensión  superficial 


Hasta  aquí  hemos  dejado  de  lado  los  efectos  debidos  a  la  existencia  de  una  superficie  de 
separación  entre  los  cuerpos.  Nuestra  definición  de  sistema  simple  requería  que  el  mismo  fuera 
suficientemente  grande  para  poder  despreciar  efectos  superficiales.  En  realidad,  las  fases  en 
contacto  están  separadas  por  una  delgada  capa  de  transición  pero  como  no  nos  interesamos 
por  su  estructura  podemos  considerarla  como  una  superficie  geométrica.  Al  aumentar  las 
dimensiones  de  un  cuerpo  (N)  los  efectos  de  superficie  crecen  mucho  más  lentamente  (N2)  que 
los  de  volumen  (N3),  por  lo  tanto  se  justifica  prescindir  de  aquellos  en  un  estudio  de  las 
propiedades  volumétricas.  Sin  embargo,  hay  una  serie  de  fenómenos  importantes  asociados 
con  las  propiedades  de  las  superficies  de  separación:  diferencia  de  presión  superficial, 
adsorción,  humectación,  capllarldad,  forma  de  las  fases,  nucleación  en  los  cambios  de  fase, 
entre  otros,  que  merecen  ser  tenidos  en  cuenta.  El  creciente  desarrollo  de  sistemas 
nanoestructurados:  materiales  nanocrlstallnos,  nanopartículas,  láminas  delgadas,  etc.  requiere 
de  una  descripción  que  tenga  en  cuenta  las  propiedades  superficiales. 

Microscópicamente,  las  fases  se  distinguen  por  el  tipo  e  intensidad  de  las  fuerzas  entre 
partículas  constituyentes:  cohesión.  En  la  superficie  de  separación  ocurre  entonces  una 
discontinuidad  en  tales  fuerzas  de  cohesión,  que  resultará  en  una  fuerza  efectiva  sobre  las 
partículas  superficiales.  Imaginemos  la  situación  en  el  seno  de  una  fase  líquida  y  en  la 
proximidad  de  la  superficie  de  separación  con  una  fase  gaseosa,  donde  las  fuerzas  de 
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interacción  son  despreciables.  Se  advierte  rápidamente  que  la  resultante  de  las  fuerzas 
atractivas  que  ejercen  las  partículas  vecinas  sobre  una  dada  partícula  en  la  fase  condensada, 
depende  de  su  distancia  a  la  superficie.  En  efecto,  en  un  punto  bien  en  el  interior  de  la  fase,  la 
resultante  es  prácticamente  nula,  en  cambio  en  puntos  próximos  a  la  superficie  hay  una  fuerza 
neta  dirigida  hacia  el  interior.  Esta  fuerza  resultante  Irá  disminuyendo  rápidamente  al  penetrar 
en  el  líquido.  En  consecuencia  para  llevar  una  partícula  a  la  superficie  hay  que  realizar  trabajo 
contra  estas  fuerzas.  Luego,  las  partículas  de  la  superficie  tienen  una  energía  adicional  o  bien 
la  superficie  tiene  energía  en  exceso:  energía  superficial. 

Las  propiedades  termodinámicas  de  una  superficie  de  separación  quedan  totalmente 
caracterizadas  por  una  magnitud:  la  tensión  superficial  y  (función  de  estado)  que  puede 
definirse  a  través  de  un  término  adicional  de  trabajo  en  la  expresión  para  la  energía  interna 
(ecuación  de  Glbbs).  Sea  A  el  área  de  la  superficie  de  separación  entre  las  fases.  En  un 
proceso  cuaslestátlco  en  que  ésta  se  modifica  en  dA  el  trabajo  invertido  es  proporcional  a  este 
cambio,  es  decir: 

dWsup  =  y.dA. 

Sobre  esta  base  se  puede  demostrar  que  sobre  la  unidad  de  longitud  del  contorno  que  limita 
una  porción  de  la  superficie  de  separación,  actúa  una  fuerza  equivalente  Igual  en  valor  absoluto 
a  y,  tangente  a  la  superficie  y  dirigida  según  la  normal  interior  al  contorno.  Esto  Implica  y>0. 
Esta  fuerza  tiende  a  minimizar  la  superficie  de  separación  entre  dos  fases,  para  un  dado 
volumen  de  ambas.  Caso  contrario  ésta  crecería  sin  límite,  entremezclándose  las  fases. 

Si  una  fase  Isótropa  se  introduce  en  otra,  toma  entonces  una  forma  esférica  si  pueden 
despreciarse  fuerzas  exteriores  volumétricas  como  la  gravltatorla.  Por  ejemplo,  gotas  o 
burbujas  pequeñas,  líquido  en  ausencia  de  gravedad,  gota  de  aceite  en  el  seno  de  una  mezcla 
de  alcohol  y  agua  que  tenga  su  misma  densidad.  Para  fases  anisótropas  (sólido  cristalino),  y 
será  función  de  la  orientación  de  las  caras  del  cristal  y  habrá  que  definir  los  correspondientes 
yhki  (hkl:  índices  de  Mlller  que  determinan  el  tipo  de  plano  cristalino).  Conocida  esta 
dependencia  podrá  determinarse  la  forma  de  equilibrio  de  esa  fase  como  aquella  que  minimiza 

la  energía  superficial  total:  Jy.dA.  Se  puede  demostrar  que  la  forma  de  equilibrio  está 

constituida  por  un  pequeño  número  de  partes  planas  (planos  cristalinos  de  bajos  índices  de 
Miller)  que  no  se  cortan  formando  ángulos  vivos  sino  que  están  unidos  por  superficies  curvadas 
(muchos  planos  pequeños  de  índices  de  Miller  superiores). 

La  tensión  superficial  entre  fases  de  una  dada  sustancia  es  función  únicamente  de  la 
temperatura.  Consideremos,  por  ejemplo,  la  superficie  de  separación  de  las  fases  líquido  y 
vapor  de  una  sustancia  pura.  Dado  que  y  caracteriza  a  fases  que  coexisten  en  equilibrio  mutuo, 
es  decir  con  valores  de  T  y  P  sobre  la  curva  de  coexistencia  P=P(T),  dependerá  de  una  única 
variable  independiente:  la  temperatura,  por  ejemplo.  Al  aumentar  la  temperatura  y  acercarnos 
al  punto  crítico,  donde  la  diferencia  entre  las  fases  desaparece,  y  debe  tender  a  cero.  Existen 
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algunas  fórmulas  propuestas  para  esta  dependencia12,  basadas  en  la  ley  de  estados 
correspondientes  f(P/Pc,V/Vc,T/Tc)  =  0.  Las  más  conocidas  son  las  siguientes: 
y  =  c.(Tc-T)/v2/3  Eotvos, 

y  =  d.Pc2/3Tc1/3(1-T/Tc)n ;  n~1.23  van  der  Waals. 

A  continuación  consideramos  algunos  casos  simples  en  los  que  el  término  adicional  y.dA 
puede  ser  relevante  para  determinar  las  propiedades  termodinámicas  del  sistema. 


Superficie  sin  espesor 


El  caso  más  sencillo  que  podemos  considerar  es  el  de  una  superficie  sin  espesor.  V  no  tiene 
sentido  en  este  caso  y  dejando  de  lado  N,  las  únicas  variables  a  considerar  serán  S  y  A: 

dll  =  T.dS  +  y.dA, 

o  bien: 

dS  =  (1/T)dU  -  (y/T).dA. 

Así,  en  las  representaciones  energética  y  entrópica,  la  variable  extensiva  A  tiene  a  y  y  a  -y/T 
como  Intensivas  conjugadas,  respectivamente. 

Las  energías  libres  pueden  definirse  de  la  manera  usual: 

F  =  G  =  U-TS, 
dG  =  -S.dT  +  y.dA, 


con  la  relación  de  Maxwell  correspondiente: 

(3S/cA)t  =  -(dy/3T)A. 

y  las  ecuaciones  de  estado: 

S  =  S(T,A)  =  A.s(T), 
y  =  y(T,A)  =  y(T). 

Aquí  hemos  usado  la  extenslvldad  de  S  y  la  intenslvldad  de  y  y  hemos  introducido  la  entropía 
específica  (por  unidad  de  área)  que  puede  entonces  escribirse: 

s  =  (8S/dA)j  =  -(5y/3T)A  =  -dy/dT 

12  Para  un  líquido  en  equilibrio  con  su  vapor,  teniendo  en  cuenta  que  [y]=erg/cm2=  [T1,3p2/3],  Landau  propone  y  = 
Tc1,3Pc2,3.f(T/Tc)=  TcVc  .g(T/Tc). 
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o  bien: 

S  =  -Ady/dT. 

Notemos  que  la  positividad  de  S  es  consistente  con  el  carácter  decreciente  de  y(T). 
Además,  utilizando  la  ecuación  de  Euler: 

U  =  TS  +  yA  =  (y  -  T.dy/dT)A, 

F  =  G  =  yA. 

Podemos  definir  también  una  capacidad  calorífica  a  área  constante: 
cA  =  T(5s/5T)a  =  -Td2y/dT2. 

Como  aplicación  podríamos  considerar  las  condiciones  de  equilibrio  en  un  sistema 
bidimensional  homogéneo  aislado.  Si  lo  suponemos  virtualmente  dividido  en  dos  subsistemas 
(1  y  2)  que  coexisten  en  equilibrio  es  fácil  encontrar  las  condiciones: 


Ti  -  T2  y  yi  —  y2. 


minimizando  la  energía  interna  para  los  dados  valores  de  S  y  A.  Arribamos  a  la  conclusión 
obvia  de  que  T  y  y  deben  ser  los  mismos  sobre  todo  el  sistema.  Podría  discutirse  en  este  marco 
la  posible  coexistencia  de  dos  fases  bidimensionales  en  equilibrio.  Para  que  ello  sea  posible 
debe  cumplirse:  y-i(T)  =  y2(T),  lo  cual  se  satisfará  a  una  temperatura  determinada  Tc. 


Sistema  tridimensional  con  superficie 


Consideremos  ahora  un  sistema  tridimensional  ordinario  compuesto  de  una  sustancia  pura 
con  una  superficie  límite.  La  ecuación  de  Gibbs  entonces  debe  escribirse: 

dU  =  TdS  -  PdV  +  pdN  +  ydA, 

y  la  correspondiente  ecuación  de  Euler: 

U  =  TS  -  PV  +  pN  +  yA. 

Se  suelen  separar  las  diferentes  magnitudes  termodinámicas  en  términos  aditivos  de 
volumen  y  de  superficie:  X  =  Xv  +  Xs.  Para  ello  debe  hacerse  alguna  hipótesis  sobre  la 
subdivisión  de  materia  entre  estas  dos  partes.  Habitualmente  se  supone  ( hipótesis  de  Gibbs) 
que  la  superficie  no  tiene  volumen  ni  partículas  es  decir:  Vs  =  Ns  =  0.  En  ese  caso  la  subdivisión 
es  inmediata  a  partir  de  que: 
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Fs  -  Gs  -  yA 


Ss  =  -  A.dy/dT, 


con  lo  que  resulta,  por  ejemplo: 

U  =  Uv  +  Us  =  Uv  +  (y  -  Tdy/dT)A, 

C  =  Cv  +  Cs  =  Cv  -  T.(d2y/dT2)A. 

Resultarán  ademas,  nuevas  relaciones  de  Maxwell.  Por  ejemplo,  en  la  representación  de 
Glbbs  a  N=cte,  como  dG  =  -SdT  +  VdP  +  ydA,  tendremos: 

(8SldA)jiP  =  -(oy/oT)PA 
(cV/3A)t,p  =  -(0y/5P)T,A- 

Esta  última  derivada  en  general  se  anula. 

Para  este  tipo  de  sistemas,  que  en  realidad  se  encuentran  abiertos  al  medio  (a  través  de  su 
superficie  límite)  y  entonces  en  equilibrio  con  él,  puede  ser  conveniente  trabajar  con  el  gran 
potencial  Q  convenclonalmente  definido  como  Q  =  U[T,p]  =  U  -  TS  -  pN.  De  acuerdo  con  lo 
visto: 

Q  s  =  U  s  -  TS  s  -  pN  s  =  Fs  =  yA 
pues  Ns  =  0.  Entonces: 

Q  =  Qv  +  Qs  =  -PV  +  yA 


mostrando  cierta  simetría  en  los  roles  de  y  y  -P  en  los  términos  de  superficie  y  volumen.  La 
utilidad  de  este  potencial  se  deriva  de  que  en  equilibrio  T  y  p  coincidirán  con  los  del  medio, 
pero  no  ocurrirá  necesariamente  lo  mismo  con  P,  como  veremos. 

Para  apreciar  la  Importancia  relativa  de  los  términos  de  superficie,  según  el  tamaño  del 
sistema,  podemos  mencionar  algunas  estimaciones  en  agua  a  373K.  La  contribución  superficial 
a  la  capacidad  calorífica  Cs/Cv  que  es  sólo  de  2.3.1 0"8  para  una  esfera  de  lOOmm  de  diámetro, 
puede  alcanzar  a  2.3.1 0"2  para  gotas  de  lOOnm  de  diámetro. 


Equilibrio  de  fases  en  presencia  de  interfase 


Consideramos  ahora  las  condiciones  de  equilibrio  para  un  sistema  aislado  compuesto  de 
dos  fases  coexlstentes,  teniendo  en  cuenta  la  capa  superficial  de  separación  entre  ambas: 
inferíase.  Suponemos  entonces  al  sistema  dividido  en  tres  subsistemas:  fase  1,  fase  2  y 
superficie  de  separación,  para  la  cual  nuevamente  Imponemos:  Vs  =  Ns  =  0.  La  condición  de 
equilibrio  será  la  de  U  mínimo  y  las  condiciones  de  cierre:  S,  V  y  N  constantes.  De  aquí 
resultan  las  relaciones  siguientes: 
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dSs  =  -  dS-i  -  dS2 
dV2  =  -  dV, 
dN2  =  -dN1 

En  un  apartamiento  virtual  del  estado  de  equilibrio: 
dll  =  dll-i  +  dl)2+  dlls  =  0. 

di)  =  -  P1dV1+  |x1dN1  +  T2dS2  -  P2dV2+  p2dN2  +  TsdSs  +  ydA 

dU  =  (Ti-Ts)dS1  +  (T2-Ts)dS2  -  (P,-  P2)dV.,+  (|i1-|i2)dN1  +  ydA. 

Como  dSi,  dS2  y  dN-i  son  independientes,  sus  coeficientes  deben  anularse  para  asegurar 
dl)=0.  Entonces  resultan  las  condiciones  de  equilibrio: 

Ti=T2=T8 

Pl  =  P2- 

Por  el  contrario  dV-,  y  dA  no  son  independientes  ya  que  los  cambios  de  volumen  y  superficie 
pueden  estar  relacionados.  Entonces13: 

-  (P1-P2)dV1  +  ydA  =  0 

obteniendo  así  la  nueva  condición  de  equilibrio: 

AP  =  (P  -| — P  2)  =  y(dA/dV1) 

En  general,  para  una  superficie  descripta  por  dos  radios  de  curvatura  (p3i  pb),  resulta  la 
ecuación  de  Laplace: 

(Pi-P2)  =  y(1/Pa+1/Pb) 

incluyendo  los  casos  particulares: 

plano  (pa  =  pb  =  oo)  AP  =  0 

cilindro  (pa  =  r;  pb  =  oo)  AP  =  y/r 

esfera  (pa  =  pb  =  r)  AP  =  2y/r 

Este  último  resultado  puede  obtenerse  directamente  notando  que  para  una  esfera  de  radio  r: 
dA=8jir.dr  y  dV-i=47ir2.dr 


Notar  que  esta  condición  es  equivalente  a  dü  =  -P1dV1-P2dV2+ydA  =  0 

126 


13 


Aplicaciones 


Como  aplicación  estimaremos  el  cambio  en  la  presión  8P  para  un  fluido  en  equilibrio  con  su 
vapor  a  través  de  una  superficie  curva.  En  ese  caso,  para  las  fases  que  coexisten: 

Pi(T,Pi)  =  p2(T,P2), 

donde  T  es  la  temperatura  común  de  equilibrio  y  Pi^P2.  Si  P0  es  la  presión  de  equilibrio  para 
interfase  plana  a  la  misma  T,  podemos  escribir: 

Hi(T,Pi)  -  w(T,Po)  =  p2(T,P2)  -  h2(T,P0), 

o  bien: 


5pi  =  5|í2 


Según  la  ecuación  de  Gibbs-Duhem,  a  T  constante:  dp=v.dP.  Entonces,  si  los  cambios  de 
presión  8P¡  =  P¡-P0  son  pequeños,  puede  aproximarse:  8p  *  v.SP  con  lo  cual: 


v-iSPí  =  v2SP2. 


Además: 


8Pi-8P2=  PrP2  =  AP 


Resolviendo  hallamos: 


SPi  =  AP.v2/(v2-Vi) 
SP2  =  AP.v1/(v2-v1). 


Por  ejemplo,  para  el  caso  de  gotas  de  líquido  en  vapor  (1=1;  2=g): 
V-]  <<  v2  =  RT/P2 «  RT/P0 
ya  que  consecuentemente:  SP2«SP1. 

Entonces: 

8P¡  «  AP,  8Pg  *  AP.PqV/RT, 

o  bien,  suponiendo  gotas  esféricas  de  radio  r: 
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Pi  *  Po+2y/r 


Pg  ~  Po[1+(2y/r)(V|/RT)]. 


Es  decir  que  la  presión  de  vapor  sobre  una  gota  es  mayor  que  sobre  la  superficie  plana  de 
un  líquido  y  esta  diferencia,  siguiendo  a  AP,  irá  aumentando  al  reducirse  el  tamaño  de  la  gota. 

En  el  caso  de  gotas  muy  pequeñas  5Pg/P0  ya  no  es  despreciable  y  la  aproximación  hecha 
(5p  *  v.SP)  no  es  válida,  ya  que  en  los  gases  v  depende  fuertemente  de  p.  Usando  para  el 
potencial  químico  la  expresión  p(T,P)  =  4>(T)  +  RT.InP  resulta: 


5|ig  -  |ig(T,Pg)  "  |ig(T,P0)  -  RT.In(P/Pg) 


Para  la  fase  líquida,  prácticamente  incompresible,  la  influencia  del  cambio  de  presión  es 
insignificante  y  podemos  continuar  aproximando  5|i|«  V|.8pi.  Luego: 

RT.In(P/Pg)  =  v,(PrPg)  «  v,(P,-P0)  -  AP.v, 

y  entonces: 

Pg  «  P0.exp(AP.v,/RT)  =  P0.exp[(2y/r)(v,/RT)]. 

Intercambiando  los  roles  de  gas  y  líquido,  burbujas  gaseosas  en  fase  líquida  (1  =g;  2=1), 
resultan  expresiones  semejantes  con  el  signo  cambiado. 
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